Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,727 G→T 9.8% intergenic (‑282/‑107) ylcI ← / → nohD DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7270GT 9.8% 1261.6 / 21.1 401intergenic (‑282/‑107)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (212/134);  new base T (2/52);  total (215/186)
Fisher's exact test for biased strand distribution p-value = 3.36e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 408 total aligned reads displayed.

CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580578‑580876
                                                                                                                                                       |                                                                                                                                                     
cGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATg                                                                                                                                                       >  2:7447286/1‑150 (MQ=255)
cGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATg                                                                                                                                                       >  1:6310385/1‑150 (MQ=255)
   tgtgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGtat                                                                                                                                                    <  2:4240775/150‑1 (MQ=255)
   tgtgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGtat                                                                                                                                                    <  2:5309957/150‑1 (MQ=255)
     tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                  >  2:8993316/1‑150 (MQ=255)
       aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                                <  2:11099430/150‑1 (MQ=255)
       aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                                >  2:4762896/1‑150 (MQ=255)
         aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              <  2:7813751/150‑1 (MQ=255)
         aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  2:8637015/1‑150 (MQ=255)
           gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            <  1:2446979/150‑1 (MQ=255)
           gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            >  2:8268952/1‑150 (MQ=255)
            tAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                           >  1:10704208/1‑150 (MQ=255)
                tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                       >  2:3538939/1‑150 (MQ=255)
                tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                       <  1:4230808/150‑1 (MQ=255)
                 gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                      <  2:6624225/150‑1 (MQ=255)
                    acttaataTCAATTTGTGGGTCCTATCCGGCGAACCGACAGGTTACGGGGCGGCGACCTCGCGGGTTAT‑CGCTATATATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATCTATTTATTTAAAATaccacc                                                                                                                                   >  2:10304619/3‑150 (MQ=14)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   >  2:9624317/1‑150 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcaccc                                                                                                                                   >  1:5983347/1‑150 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAacccc                                                                                                                                   >  2:7549685/1‑150 (MQ=255)
                     attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  <  2:6207279/150‑1 (MQ=255)
                       tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                >  2:6889454/1‑150 (MQ=255)
                         ttatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                              <  1:5636883/150‑1 (MQ=255)
                         ttatCATTGTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                              <  2:7113200/150‑1 (MQ=255)
                          tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             <  2:2984254/150‑1 (MQ=255)
                            tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                           >  2:3706663/1‑150 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         >  2:838023/1‑150 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:1659494/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         >  2:9478678/1‑150 (MQ=255)
                               ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        <  2:9327518/150‑1 (MQ=255)
                                 ttGTGGGTGGTTTCCGGCGATCCGAGTGGTTAGGGGGCGGCGTCCTCGCGGGTTTT‑GGCTATTTTTGTAATTTTTGGGGTTTTAGGTGTTTCCGT‑TCTTCTTCGTCGTTAGTTAATGTATTTATTTAAAATTCCCCCTGTaaagtatggt                                                                                                                      >  2:234559/1‑144 (MQ=14)
                                      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTCATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACgg                                                                                                                 >  2:6035065/1‑149 (MQ=255)
                                      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 >  2:1449820/1‑150 (MQ=39)
                                        tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACa                                                                                                               >  2:8204232/1‑150 (MQ=39)
                                         ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAg                                                                                                              >  1:7552320/1‑150 (MQ=39)
                                           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:3332247/1‑150 (MQ=39)
                                           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:345572/1‑150 (MQ=39)
                                           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:1090190/1‑150 (MQ=39)
                                            ttCGGGCGATCCGAATGTTTGGGGGGCGGCGACCTCGGGGTTTTT‑CGTTTTTTTTGAAATTTTTCCGTTTTATGTTTTTTCCGT‑TTTTCTTGGTCTTAATTTATTTTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTg                                                                                                           <  2:3987780/150‑1 (MQ=255)
                                              ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                                         >  2:8688254/1‑150 (MQ=39)
                                                  cGATCCGACTGGTTACGGGGCGGCGACCTCGCGGTTTTT‑CGCTATTTATGAAATTTTTCCGGTTTATGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                                     <  2:9779960/150‑1 (MQ=255)
                                                    aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  1:5576080/1‑150 (MQ=38)
                                                     tCCTACAGGTTCCGGGGCGGCGACCTGGCGGGTTTT‑CGCTATTTATGAAAATTTTCGGTTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCg                                                                                                  <  2:453037/150‑1 (MQ=14)
                                                     tCCGACTGGTTACGGGGGGGCGACCTCGCGGGTTTT‑CGATTTTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCg                                                                                                  <  2:3436307/150‑1 (MQ=14)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  1:9828688/1‑150 (MQ=39)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:8550768/1‑150 (MQ=39)
                                                           aGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            >  1:1184559/1‑150 (MQ=39)
                                                           aGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            >  1:7404492/1‑150 (MQ=39)
                                                             gTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTTTTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:4432541/150‑1 (MQ=255)
                                                                aCGGGGCGGCGACCTCGCGGGTTTT‑CGCGATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCGGAAAAGAAAGGAAACGACCGGTGCTGCAAGCGCGCTTTTTGGc                                                                                       >  1:2838619/1‑150 (MQ=14)
                                                                  ggggCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                                     >  2:6670552/1‑150 (MQ=255)
                                                                    ggcggcGACCTCGCGGTTTTT‑CGATATTTATTAAAATTTTCCGTTTTAAGGTTTTTCCGT‑TCTTCTTCGTCGTAACTTAATTTTTTTATTTAAAATACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctct                                                                                   <  2:10494346/150‑1 (MQ=11)
                                                                    ggcggcGACCTCGCGGGTTTT‑CGCTATTTATGAAAAATTTACGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctct                                                                                   >  1:3871986/1‑150 (MQ=255)
                                                                     gcggcgACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTg                                                                                  >  2:7054822/1‑150 (MQ=255)
                                                                        gcgACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAAAGTATTTATTTAAAATACCCCCTGAAAAGAAACGAACCGACCGGTGCGGAAAGCGAGCTTTTTGGCCGCTGTCg                                                                               >  1:9432777/1‑150 (MQ=14)
                                                                         cgACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTGATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                              >  2:10754197/1‑150 (MQ=255)
                                                                             cTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTcc                                                                          <  2:2652348/150‑1 (MQ=255)
                                                                               cgggggTTTTCCGCTATTTATTAAAATTTGCGGTTTTTTGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATTTTTTTATTAAAAAACCCCCCTTAAAAGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         <  1:7092706/147‑1 (MQ=11)
                                                                               cgcgGTTTTT‑CGCTATTTATCAAAATTTTCCGGTATAAGGTGTTACCGT‑TCTTCTTCGTCTTAACTTAATTTATTTATTTAAAATACCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCtt                                                                        <  2:10704208/150‑1 (MQ=14)
                                                                                       tt‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                >  2:10829692/1‑150 (MQ=255)
                                                                                           gCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                             >  1:4535915/1‑150 (MQ=255)
                                                                                            cTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTg                                                            >  2:5511457/1‑150 (MQ=255)
                                                                                             tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  1:9063860/1‑150 (MQ=255)
                                                                                              tttttTGAAATTTTTCGGTTTTTTGGTTTTTCCGT‑TATTCTTCGTCGTAATTTAATTTTTTTATTAAAAAACCCCCCTCAAAATAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTc                                                          <  1:1792629/149‑1 (MQ=11)
                                                                                                    gAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  2:3149312/1‑150 (MQ=255)
                                                                                                     aaatTTTCCCTGTTTCTGGTGTTCCCGT‑TCTTCTTCGTCGTAACTTAATTTATTTTTTTAAAATACCCCCCTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   <  1:8581313/150‑1 (MQ=14)
                                                                                                     aaaaTTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   >  2:1959560/1‑150 (MQ=255)
                                                                                                         ttttCCGTTTTCTGGT‑TTTGCGTGTCTTCTTCGTCTTAACTTATTTTTTTTATTAAAAAACCCCCCTTAAAAGAAAGCACACGACAGGTGCTGAAAGCGAGCCTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  1:4863546/150‑1 (MQ=11)
                                                                                                         ttttCCGGTTTAAGGTGTTTCCTT‑TCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  2:5622878/150‑1 (MQ=255)
                                                                                                         ttgtCCGTTTTATGTTGTGTCCTT‑TATTATTCTTCGTACCTTATTTTATTAATTTAAAACACCCCCTTAAAATAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGTAtgtg                                               <  1:9963749/147‑1 (MQ=14)
                                                                                                            tCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTGCGTGGAATGTGCaa                                            >  1:4733163/1‑150 (MQ=255)
                                                                                                               ggTTTATGGTTTTTCCTT‑TATTCTTCTTCTTAACTTAATTTATTTATTTAAAATACCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         <  1:5354674/150‑1 (MQ=14)
                                                                                                                gTTTAAGGTGTTTCCTT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                        <  2:1644237/150‑1 (MQ=255)
                                                                                                                gTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTTTCCGTGGACTGTGCAATgta                                        >  2:5712264/1‑148 (MQ=255)
                                                                                                                gTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                        >  1:2747898/1‑150 (MQ=255)
                                                                                                                 tttAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  2:9143953/1‑150 (MQ=255)
                                                                                                                 tttAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  2:1321563/1‑150 (MQ=255)
                                                                                                                     tggTGTTTCGGT‑TCTTCTTCGTCTTAACTTATTGTATTTATTAAAAATACCCCCTGAAAACAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   <  1:6926138/149‑1 (MQ=14)
                                                                                                                     aGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   >  1:6723736/1‑150 (MQ=255)
                                                                                                                       gtgtTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                 >  2:4985508/1‑150 (MQ=255)
                                                                                                                        tgtTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAAAGGAAGTcaaca                                >  1:432535/1‑150 (MQ=255)
                                                                                                                         gtTTCCGT‑ACTTCTTCGTCGTAAATAAATGTATTTATTAAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaacaa                               <  2:5858265/150‑1 (MQ=255)
                                                                                                                           ttCCGT‑TGTTGTTGGTCGTAACTTATTTTATTTATTTAAAATACCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             <  1:8268952/150‑1 (MQ=14)
                                                                                                                            tcggT‑TCTTCTTCGTCTACACTTAATTTATTTATTTAAAATACCCCCTCAAAACAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                            <  1:8486242/147‑1 (MQ=14)
                                                                                                                               cttactccttcttCGTATTTTATTTTATTATTTAAAAATACCCCCTCAAAAGAAACAAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagc                          <  1:1006280/143‑1 (MQ=14)
                                                                                                                                t‑tcttcttcGTCTTAACTTAATGTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcag                        <  1:8083023/150‑1 (MQ=255)
                                                                                                                                       ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  >  2:5112639/1‑150 (MQ=255)
                                                                                                                                        tcttcttAAATTAATTTTTTTATTAAAAATACCCCATTAAAACAAAGGAAACGACTGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                 <  1:10538420/144‑1 (MQ=14)
                                                                                                                                         cgtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                >  1:415918/1‑150 (MQ=255)
                                                                                                                                         cgtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                <  2:5364621/150‑1 (MQ=255)
                                                                                                                                           tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa              >  1:7802213/1‑150 (MQ=255)
                                                                                                                                            cgtAACCTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt             >  2:10235775/1‑150 (MQ=255)
                                                                                                                                                aCTTAATTTATTTATTTAAAATACCCCCTTAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTc         <  1:496180/150‑1 (MQ=255)
                                                                                                                                                aCTTAATTTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTc         <  1:11000368/150‑1 (MQ=255)
                                                                                                                                                aCTTAATGTATTCATTTAAAAACCCCCCTCAAACCAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTc         <  1:9624317/150‑1 (MQ=14)
                                                                                                                                                 cTTACTGTATTTATTTAAAATACCCCCTGAAACGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg        <  2:8034649/150‑1 (MQ=255)
                                                                                                                                                 cTTAATGTATTTATTTAAAATACCCCCTGAAACGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg        <  2:1060123/150‑1 (MQ=255)
                                                                                                                                                 cTTAATGTATTTATTTAAAAAACCCCCTGAAACCAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg        <  2:10563709/150‑1 (MQ=255)
                                                                                                                                                  ttAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCgg       >  1:8374996/1‑150 (MQ=255)
                                                                                                                                                  ttAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCgg       >  2:8936153/1‑150 (MQ=255)
                                                                                                                                                    aaTGTATTTATTTAAAATACCCCCTGAAAAGAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg     <  2:4998347/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGT‑TCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580578‑580876

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑