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.

CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580578‑580876
                                                                                                                                                       |                                                                                                                                                     
cGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                       >  2:7447286/1‑150 (MQ=255)
cGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                       >  1:6310385/1‑150 (MQ=255)
   tgtgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtat                                                                                                                                                    <  2:4240775/150‑1 (MQ=255)
   tgtgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtat                                                                                                                                                    <  2:5309957/150‑1 (MQ=255)
     tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                  >  2:8993316/1‑150 (MQ=255)
       aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                                <  2:11099430/150‑1 (MQ=255)
       aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                                >  2:4762896/1‑150 (MQ=255)
         aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              <  2:7813751/150‑1 (MQ=255)
         aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  2:8637015/1‑150 (MQ=255)
           gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            <  1:2446979/150‑1 (MQ=255)
           gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            >  2:8268952/1‑150 (MQ=255)
            tAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                           >  1:10704208/1‑150 (MQ=255)
                tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                       >  2:3538939/1‑150 (MQ=255)
                tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                       <  1:4230808/150‑1 (MQ=255)
                 gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                      <  2:6624225/150‑1 (MQ=255)
                    acttaataTCAATTTGTGGGTCCTATCCGGCGAACCGACAGGTTACGGGGCGGCGACCTCGCGGGTTATCGCTATATATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATCTATTTATTTAAAATaccacc                                                                                                                                   >  2:10304619/3‑150 (MQ=14)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   >  2:9624317/1‑150 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcaccc                                                                                                                                   >  1:5983347/1‑150 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAacccc                                                                                                                                   >  2:7549685/1‑150 (MQ=255)
                     attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  <  2:6207279/150‑1 (MQ=255)
                       tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                >  2:6889454/1‑150 (MQ=255)
                         ttatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                              <  1:5636883/150‑1 (MQ=255)
                         ttatCATTGTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                              <  2:7113200/150‑1 (MQ=255)
                          tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             <  2:2984254/150‑1 (MQ=255)
                            tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                           >  2:3706663/1‑150 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         >  2:838023/1‑150 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:1659494/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         >  2:9478678/1‑150 (MQ=255)
                               ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        <  2:9327518/150‑1 (MQ=255)
                                 ttGTGGGTGGTTTCCGGCGATCCGAGTGGTTAGGGGGCGGCGTCCTCGCGGGTTTTGGCTATTTTTGTAATTTTTGGGGTTTTAGGTGTTTCCGTTCTTCTTCGTCGTTAGTTAATGTATTTATTTAAAATTCCCCCTGTaaagtatggt                                                                                                                      >  2:234559/1‑144 (MQ=14)
                                      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTCATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACgg                                                                                                                 >  2:6035065/1‑149 (MQ=255)
                                      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 >  2:1449820/1‑150 (MQ=39)
                                        tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACa                                                                                                               >  2:8204232/1‑150 (MQ=39)
                                         ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAg                                                                                                              >  1:7552320/1‑150 (MQ=39)
                                           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:3332247/1‑150 (MQ=39)
                                           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:345572/1‑150 (MQ=39)
                                           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:1090190/1‑150 (MQ=39)
                                            ttCGGGCGATCCGAATGTTTGGGGGGCGGCGACCTCGGGGTTTTTCGTTTTTTTTGAAATTTTTCCGTTTTATGTTTTTTCCGTTTTTCTTGGTCTTAATTTATTTTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTg                                                                                                           <  2:3987780/150‑1 (MQ=255)
                                              ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                                         >  2:8688254/1‑150 (MQ=39)
                                                  cGATCCGACTGGTTACGGGGCGGCGACCTCGCGGTTTTTCGCTATTTATGAAATTTTTCCGGTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                                     <  2:9779960/150‑1 (MQ=255)
                                                    aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  1:5576080/1‑150 (MQ=38)
                                                     tCCTACAGGTTCCGGGGCGGCGACCTGGCGGGTTTTCGCTATTTATGAAAATTTTCGGTTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCg                                                                                                  <  2:453037/150‑1 (MQ=14)
                                                     tCCGACTGGTTACGGGGGGGCGACCTCGCGGGTTTTCGATTTTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCg                                                                                                  <  2:3436307/150‑1 (MQ=14)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  1:9828688/1‑150 (MQ=39)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:8550768/1‑150 (MQ=39)
                                                           aGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            >  1:1184559/1‑150 (MQ=39)
                                                           aGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            >  1:7404492/1‑150 (MQ=39)
                                                             gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTTTTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:4432541/150‑1 (MQ=255)
                                                                aCGGGGCGGCGACCTCGCGGGTTTTCGCGATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCGGAAAAGAAAGGAAACGACCGGTGCTGCAAGCGCGCTTTTTGGc                                                                                       >  1:2838619/1‑150 (MQ=14)
                                                                  ggggCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                                     >  2:6670552/1‑150 (MQ=255)
                                                                    ggcggcGACCTCGCGGTTTTTCGATATTTATTAAAATTTTCCGTTTTAAGGTTTTTCCGTTCTTCTTCGTCGTAACTTAATTTTTTTATTTAAAATACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctct                                                                                   <  2:10494346/150‑1 (MQ=11)
                                                                    ggcggcGACCTCGCGGGTTTTCGCTATTTATGAAAAATTTACGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctct                                                                                   >  1:3871986/1‑150 (MQ=255)
                                                                     gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTg                                                                                  >  2:7054822/1‑150 (MQ=255)
                                                                        gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAAGTATTTATTTAAAATACCCCCTGAAAAGAAACGAACCGACCGGTGCGGAAAGCGAGCTTTTTGGCCGCTGTCg                                                                               >  1:9432777/1‑150 (MQ=14)
                                                                         cgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTGATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                              >  2:10754197/1‑150 (MQ=255)
                                                                             cTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTcc                                                                          <  2:2652348/150‑1 (MQ=255)
                                                                               cgggggTTTTCCGCTATTTATTAAAATTTGCGGTTTTTTGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATTTTTTTATTAAAAAACCCCCCTTAAAAGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         <  1:7092706/147‑1 (MQ=11)
                                                                               cgcgGTTTTTCGCTATTTATCAAAATTTTCCGGTATAAGGTGTTACCGTTCTTCTTCGTCTTAACTTAATTTATTTATTTAAAATACCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCtt                                                                        <  2:10704208/150‑1 (MQ=14)
                                                                                       ttCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                >  2:10829692/1‑150 (MQ=255)
                                                                                           gCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                             >  1:4535915/1‑150 (MQ=255)
                                                                                            cTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTg                                                            >  2:5511457/1‑150 (MQ=255)
                                                                                             tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  1:9063860/1‑150 (MQ=255)
                                                                                              tttttTGAAATTTTTCGGTTTTTTGGTTTTTCCGTTATTCTTCGTCGTAATTTAATTTTTTTATTAAAAAACCCCCCTCAAAATAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTc                                                          <  1:1792629/149‑1 (MQ=11)
                                                                                                    gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  2:3149312/1‑150 (MQ=255)
                                                                                                     aaatTTTCCCTGTTTCTGGTGTTCCCGTTCTTCTTCGTCGTAACTTAATTTATTTTTTTAAAATACCCCCCTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   <  1:8581313/150‑1 (MQ=14)
                                                                                                     aaaaTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   >  2:1959560/1‑150 (MQ=255)
                                                                                                         ttttCCGTTTTCTGGTTTTGCGTGTCTTCTTCGTCTTAACTTATTTTTTTTATTAAAAAACCCCCCTTAAAAGAAAGCACACGACAGGTGCTGAAAGCGAGCCTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  1:4863546/150‑1 (MQ=11)
                                                                                                         ttttCCGGTTTAAGGTGTTTCCTTTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  2:5622878/150‑1 (MQ=255)
                                                                                                         ttgtCCGTTTTATGTTGTGTCCTTTATTATTCTTCGTACCTTATTTTATTAATTTAAAACACCCCCTTAAAATAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGTAtgtg                                               <  1:9963749/147‑1 (MQ=14)
                                                                                                            tCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTGCGTGGAATGTGCaa                                            >  1:4733163/1‑150 (MQ=255)
                                                                                                               ggTTTATGGTTTTTCCTTTATTCTTCTTCTTAACTTAATTTATTTATTTAAAATACCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         <  1:5354674/150‑1 (MQ=14)
                                                                                                                gTTTAAGGTGTTTCCTTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                        <  2:1644237/150‑1 (MQ=255)
                                                                                                                gTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTTTCCGTGGACTGTGCAATgta                                        >  2:5712264/1‑148 (MQ=255)
                                                                                                                gTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                        >  1:2747898/1‑150 (MQ=255)
                                                                                                                 tttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  2:9143953/1‑150 (MQ=255)
                                                                                                                 tttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  2:1321563/1‑150 (MQ=255)
                                                                                                                     tggTGTTTCGGTTCTTCTTCGTCTTAACTTATTGTATTTATTAAAAATACCCCCTGAAAACAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   <  1:6926138/149‑1 (MQ=14)
                                                                                                                     aGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   >  1:6723736/1‑150 (MQ=255)
                                                                                                                       gtgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                 >  2:4985508/1‑150 (MQ=255)
                                                                                                                        tgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAAAGGAAGTcaaca                                >  1:432535/1‑150 (MQ=255)
                                                                                                                         gtTTCCGTACTTCTTCGTCGTAAATAAATGTATTTATTAAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaacaa                               <  2:5858265/150‑1 (MQ=255)
                                                                                                                           ttCCGTTGTTGTTGGTCGTAACTTATTTTATTTATTTAAAATACCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             <  1:8268952/150‑1 (MQ=14)
                                                                                                                            tcggTTCTTCTTCGTCTACACTTAATTTATTTATTTAAAATACCCCCTCAAAACAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                            <  1:8486242/147‑1 (MQ=14)
                                                                                                                               cttactccttcttCGTATTTTATTTTATTATTTAAAAATACCCCCTCAAAAGAAACAAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagc                          <  1:1006280/143‑1 (MQ=14)
                                                                                                                                ttcttcttcGTCTTAACTTAATGTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcag                        <  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)
                                                                                                                                                       |                                                                                                                                                     
CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  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: