Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,413,542 A→T Y67N (TAT→AAT)  ydaC ← DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,5420AT98.9% 299.2 / NA 92Y67N (TAT→AAT) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (46/45);  minor base G (0/1);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTT  >  NC_000913/1413397‑1413684
                                                                                                                                                 |                                                                                                                                              
cGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTGTTAATTTGAATTCaa                                                                                                                                             >  1:1545516/1‑149 (MQ=255)
 gAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTCTGAAGTATTAATTTGAATTCAAt                                                                                                                                            >  2:817096/1‑149 (MQ=255)
 gAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGGATTCAAt                                                                                                                                            <  1:459992/149‑1 (MQ=255)
  aaGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATg                                                                                                                                           >  2:655513/1‑149 (MQ=255)
      aTCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGATTCTTCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAAt                                                                                                                                       >  1:1071200/1‑149 (MQ=255)
        cAGGCCGTATTCCACCGTCCAATCTTCCTTACAAATACCATGTGACAACTGCTATCGTCAGTCAGACGCTAATTGTAGACGGTGATTGACTCGCGATATTAATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATtc                                                                                                                                     <  2:136903/149‑1 (MQ=255)
          ggCCGCATTCAACCCTCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATtctc                                                                                                                                   <  1:188080/149‑1 (MQ=255)
            ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                 <  1:527959/149‑1 (MQ=255)
            ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                 <  2:809720/149‑1 (MQ=255)
            ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                 <  1:97366/149‑1 (MQ=255)
            ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAACCAGCTATACAGGTTTAATTTGAATTAAATGCAAGTCTCGa                                                                                                                                 >  1:323956/1‑149 (MQ=255)
                 ttCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGCCTAAAATTGACGCCCGATGTTGAGTATTAATAATCAGCTATGAAGTTTTAATTTGAGTTGAATGCAATTCTCGAGGACt                                                                                                                            >  1:1249383/1‑149 (MQ=255)
                 ttCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACt                                                                                                                            >  2:275885/1‑149 (MQ=255)
                 ttCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCCAGGACt                                                                                                                            >  1:965121/1‑149 (MQ=255)
                    aacgcaacACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                         <  2:548974/142‑1 (MQ=255)
                    aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                         >  1:348951/1‑149 (MQ=255)
                    aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                         >  1:790613/1‑149 (MQ=255)
                      ccATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGAGTGAAGt                                                                                                                       >  2:310939/1‑149 (MQ=255)
                      ccATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAACCGCCAATTTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCGCGAGAACGGAAGt                                                                                                                       >  2:480103/1‑149 (MQ=255)
                       cATCCTCTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGtt                                                                                                                      <  2:783621/149‑1 (MQ=255)
                       cATCCACTCTTCATTAAAGAGGACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAACTATGAAGTTTAAAGTTGAATTCAATGCAATTCTCGAGGACTGAAGtt                                                                                                                      >  2:325951/1‑149 (MQ=255)
                           cACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                  <  1:830322/149‑1 (MQ=255)
                            aCTCTTCATTAAAGAGTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCg                                                                                                                 <  2:360703/149‑1 (MQ=255)
                               ctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCaa                                                                                                              >  2:936405/1‑149 (MQ=255)
                                ttCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAt                                                                                                             >  2:786114/1‑149 (MQ=255)
                                ttCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAt                                                                                                             >  1:1349062/1‑149 (MQ=255)
                                       aaGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAAtt                                                                                                      >  1:1054734/1‑149 (MQ=255)
                                         gATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTa                                                                                                    <  2:649415/149‑1 (MQ=255)
                                          aTTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAACTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTAt                                                                                                   <  2:1148981/149‑1 (MQ=255)
                                               ttGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGtt                                                                                              >  2:1270466/1‑149 (MQ=255)
                                                tgtgCCATCGCATCACCTCTGACAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGttt                                                                                             <  2:950792/149‑1 (MQ=255)
                                                tgtgCCATCGCATCACCGCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGttt                                                                                             <  2:862646/149‑1 (MQ=255)
                                                 gtgCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGtttt                                                                                            >  1:1535791/1‑149 (MQ=255)
                                                  tgCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTa                                                                                           >  1:41211/1‑149 (MQ=255)
                                                  tgCCATCGCATCACCGCTGGCAGACGCCAATGTTAGACGGAAAGGGACTCCCGATGTTGCTTATGAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTa                                                                                           >  2:1145582/1‑149 (MQ=255)
                                                       tCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGCCGCCCGATGTTGATTATTAATAATCGGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCGATTACAATTTATACATTTTACttta                                                                                      >  1:1294765/1‑148 (MQ=255)
                                                         gCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTg                                                                                    <  2:259034/149‑1 (MQ=255)
                                                             cACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctct                                                                                <  2:1459437/149‑1 (MQ=255)
                                                              aCCTCTGGCAGGCGCCAATGTTAGACTTAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctc                                                                               <  2:263371/149‑1 (MQ=255)
                                                               cgtctGGCAGGTGGCAATGTTAGCCTGAAATTGACCCCCGTAGTGAGTTATGAATCATTAGCGATGAAGTTTTAACTTGAAATCAATGCAAGACACGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctct                                                                              <  2:712269/147‑1 (MQ=255)
                                                                ctctGGCAGGCGCCAATGTTAGACTTAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                             >  1:780393/1‑149 (MQ=255)
                                                                ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctctt                                                                             >  1:769045/1‑148 (MQ=255)
                                                                ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCAATGTTGATTATTAATAATCAGCTATGCAGTTTTAATTTGAATTCAATGCAATTCTCGAGGCCTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTCCTTTCTGCTCTCTg                                                                             >  1:163907/1‑149 (MQ=255)
                                                                  ctGGCAGGCGACAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTTTCGAGGAATGAAGTTTCTCTCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGa                                                                           >  2:425513/1‑149 (MQ=255)
                                                                    ggCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGaaa                                                                         <  1:786114/149‑1 (MQ=255)
                                                                       aGGCGCCAATGTTAGACGGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGc                                                                      <  2:505629/149‑1 (MQ=255)
                                                                           gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                  <  2:1054734/149‑1 (MQ=255)
                                                                            ccAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGATTTACTTTCTGCTCTCTGGAAACTCCTGct                                                                 >  2:825080/1‑149 (MQ=255)
                                                                                 gTTAGACTGAACTTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAACGCGATTCTCGAACCACGATGTGTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCttt                                                            <  1:1532875/149‑1 (MQ=255)
                                                                                 gTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGTCTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCttt                                                            <  1:665540/149‑1 (MQ=255)
                                                                                  atAGACAGAAATTGATGCCCGATGTTGATTATTAATAAGCAGCTATGAAGGTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCtttt                                                           <  2:1071200/148‑1 (MQ=255)
                                                                                   tAGACTGAAATTGACGCCCGATGTTGCTTCTTCCTACTCGTCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCttttt                                                          >  1:376443/1‑149 (MQ=255)
                                                                                   tAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCttttt                                                          <  1:1426162/149‑1 (MQ=255)
                                                                                    aGACTGAAATTGACGCCCGATGTTCATTATTAATGATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCtttttt                                                         >  1:785949/1‑149 (MQ=255)
                                                                                      aCTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAc                                                       >  1:473655/1‑149 (MQ=255)
                                                                                      aCTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAc                                                       >  1:629629/1‑149 (MQ=255)
                                                                                      aCTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTGAc                                                       <  2:725933/149‑1 (MQ=255)
                                                                                       cTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAcc                                                      <  1:1201654/149‑1 (MQ=255)
                                                                                        tGAAATTGACGGGGGATGTTGATTATTAGTAAGCAGCTATGAAGGTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCt                                                     >  2:1541686/1‑149 (MQ=255)
                                                                                        tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTActt                                                     >  2:130689/1‑147 (MQ=255)
                                                                                        tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAccc                                                     >  2:1452970/1‑148 (MQ=255)
                                                                                        tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCt                                                     >  1:1211790/1‑149 (MQ=255)
                                                                                        tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCt                                                     >  1:1445568/1‑149 (MQ=255)
                                                                                            aTTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGCAACCGCCTGCTTCTTTTTTACCTgacc                                                 >  2:1315400/1‑147 (MQ=255)
                                                                                                     cGATGTTGCTGCTTAATCGTCAGCTATGAGGGTTTAATTTGTATTCAATGCAATTCTCGAGGACTGCAGTTTCTCGCAATTAAACTTTATCCGTTATACTTGCTGCTTTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAtttttt                                        <  2:258379/149‑1 (MQ=255)
                                                                                                      tattttGTTTATTAATTAGCAGCTATAAAGTTTAAATTCGAATTAATTGCGATTCTCGAGCCCTGAAGTTTCTCGCAAGTAATATTTATCTGTTTTACTTTCTGCTCTCTGGAAACGCCTTCTTCTTTTTTACCTGAGAGCGTTTTTTc                                       <  2:362480/145‑1 (MQ=255)
                                                                                                      gatttagATTAGTAACAATCAGCTAAGAAGGTTTAATTTGCATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc                                       <  1:807474/143‑1 (MQ=255)
                                                                                                      gATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc                                       <  1:285390/149‑1 (MQ=255)
                                                                                                      gATGTTGATTATTAATAATCAGCTAGGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc                                       <  1:628380/149‑1 (MQ=255)
                                                                                                          ttGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCAt                                   >  1:942465/1‑149 (MQ=255)
                                                                                                            gATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGTCTGGAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACATGAGAGCATTTTTTCGCATTc                                 >  1:688131/1‑149 (MQ=255)
                                                                                                            gATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTc                                 <  2:632765/149‑1 (MQ=255)
                                                                                                              ttattaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTTGTATTCTg                               >  2:1420312/1‑149 (MQ=255)
                                                                                                                   cataatCAGCTATGAAGTTTTAATTGGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTc                          <  1:1491395/148‑1 (MQ=255)
                                                                                                                     taatCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGt                        >  2:421564/1‑149 (MQ=255)
                                                                                                                     aaatCAGCTAGGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGt                        <  1:349654/148‑1 (MQ=255)
                                                                                                                      aatCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTAATTTCGtt                       <  2:790613/149‑1 (MQ=255)
                                                                                                                          aGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGtt                   <  1:1015776/149‑1 (MQ=255)
                                                                                                                            cTATGAAGTTTTAACTTGATTTCACTGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCATTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTCCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTa                 <  2:46644/149‑1 (MQ=255)
                                                                                                                               tGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAt              <  2:965121/149‑1 (MQ=255)
                                                                                                                                   gTTTTAATTTGAATTCAATGCAATCGTCGAGGACTGAAGTTTCTCGCAATTAAGATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTg          <  1:1315400/149‑1 (MQ=255)
                                                                                                                                    ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGa         >  1:503317/1‑149 (MQ=255)
                                                                                                                                    ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGa         <  1:180014/149‑1 (MQ=255)
                                                                                                                                     tttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGaa        <  2:348951/149‑1 (MQ=255)
                                                                                                                                      ttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAACTTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTGTAGATTTTGAat       >  2:517887/1‑149 (MQ=255)
                                                                                                                                      ttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAat       >  2:1550516/1‑149 (MQ=255)
                                                                                                                                         aatagAATTCAATGCAATTATCGAGGAATGAAGTTTATCGCAATTCCAATTTATCAGTTTTACTTTCTGCTCTCCGCAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATc    <  1:325951/145‑1 (MQ=255)
                                                                                                                                         aTTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTCGTTTAGATTTTGAATATc    >  1:954612/1‑149 (MQ=255)
                                                                                                                                         aTTTGAATGCAATGCAAGTCTCGGAGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATc    <  2:1294765/149‑1 (MQ=255)
                                                                                                                                           ttGAATTCAATGCAATTCTCGGGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCCTTAGTTTAGATTTTGAATATCtt  <  1:816279/149‑1 (MQ=255)
                                                                                                                                           ttGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTTGTTTTTTTTTTTTGTTGGTTTAGATTTTGAATATCtt  >  2:6416/1‑149 (MQ=255)
                                                                                                                                           ttGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCtt  <  2:1462473/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTT  >  NC_000913/1413397‑1413684

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: