Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,727 G→T 5.9% 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 5.9% 1399.1 / 13.1 449intergenic (‑282/‑107)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (241/161);  new base T (2/45);  total (243/206)
Fisher's exact test for biased strand distribution p-value = 1.29e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 459 total aligned reads displayed.

CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580578‑580876
                                                                                                                                                      |                                                                                                                                                     
cGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                       <  2:3834888/150‑1 (MQ=255)
  ttGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGta                                                                                                                                                     >  2:4206988/1‑150 (MQ=255)
     tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                  >  2:7875042/1‑150 (MQ=255)
         aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  1:7892111/1‑150 (MQ=255)
             aGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                          <  1:8895034/150‑1 (MQ=255)
               atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                                        <  2:9374712/150‑1 (MQ=255)
               atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                                        >  2:1263705/1‑150 (MQ=255)
                tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                       >  2:9070378/1‑150 (MQ=255)
                 gatACTGATTCTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                      <  2:1109832/150‑1 (MQ=255)
                 gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                      >  1:2614891/1‑150 (MQ=255)
                  atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     <  1:9186619/150‑1 (MQ=255)
                  atAAATATTATCATTTTGTGGGTCCTTACCGGCGATCCGACAGGTTACTGGGCGGCGACCTCGCGGGTTTAAGCTATTTATGAAAAATTTCCGGTTTAA‑GGTGTTTCCGTTCTACTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     <  2:5817872/150‑1 (MQ=17)
                     attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  >  2:1612361/1‑150 (MQ=255)
                     attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  >  1:8017300/1‑150 (MQ=255)
                          tATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAA                                                                                                                             <  2:3413382‑M1/150‑1 (MQ=255)
                            tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGA                                                                                                                           >  2:7797023‑M1/1‑150 (MQ=255)
                             cATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAA                                                                                                                          >  2:7192400‑M1/1‑150 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAA                                                                                                                         <  2:6886305‑M1/150‑1 (MQ=255)
                                    tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAC                                                                                                                   <  1:3858178‑M1/150‑1 (MQ=255)
                                       gTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAC                                                                                                                >  1:2813624‑M1/1‑150 (MQ=255)
                                        tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACA                                                                                                               >  1:2022686‑M1/1‑150 (MQ=255)
                                         ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAG                                                                                                              >  1:4555414‑M1/1‑150 (MQ=255)
                                              gcGGCGATCCGACAGGTTACGGGGCGGCGTCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTTTTTTATTTACACCCTGATTTGAAAGGAAACGACAGGTGCt                                                                                                         >  2:8641317/2‑150 (MQ=14)
                                              ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTGGGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGTAATGAAAGGAAACGACAGGTGCt                                                                                                         >  2:2767271/1‑150 (MQ=255)
                                              ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                                         >  2:9179445/1‑150 (MQ=39)
                                               cGGCGATCGGACAGGTTCGGGGGCGGCGACCTCGGGGTTTTTCGCTTTTTATGAAATTTTTCCGGTTTCT‑GGTGTTTCCGTTTTTCTTGGTCTTAAATTATTTTATTTTTTTAAAATACCCCCTTAAAAGAAAGGAAACGACAGGTGCTg                                                                                                        <  2:4309289/150‑1 (MQ=14)
                                                ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                       >  1:34497/1‑150 (MQ=39)
                                                 gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaa                                                                                                      >  2:4152559/1‑150 (MQ=38)
                                                    aTCCGACTGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTTTTTCCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   <  2:6062622/150‑1 (MQ=255)
                                                    aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  1:8404208/1‑150 (MQ=38)
                                                      ccGACTGGTTACGGGGGGGCGACCTCGGGGTTTTTCGATATTTATGAAAATTTTCGGTTTTAT‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTTTTTATTTAAAATACCCCCGGAAAACAAAGGAAACGACAGGTGCTGAAAGCGa                                                                                                 <  2:5756344/150‑1 (MQ=11)
                                                       cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                                >  1:469558/1‑150 (MQ=39)
                                                       cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCGGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                                >  2:4995705/1‑150 (MQ=255)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:2909161/1‑150 (MQ=39)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  1:3564666/1‑150 (MQ=39)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  1:6776111/1‑150 (MQ=39)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTGTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:6686407/1‑150 (MQ=255)
                                                         aCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCt                                                                                              >  1:6880537/1‑150 (MQ=39)
                                                            ggTTAGGGGGCGGCGACCTCGCGGGTTTTCGCTTTTTCTGAAATTTTTCCGGTTTAT‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCtttt                                                                                           <  2:3613733/150‑1 (MQ=14)
                                                             gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:736296/150‑1 (MQ=40)
                                                              ttACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                                         >  1:5664897/1‑150 (MQ=40)
                                                                     gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctg                                                                                  <  2:7050764‑M2/150‑4 (MQ=255)
                                                                      cggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTGTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgt                                                                                 <  2:1413412‑M2/150‑5 (MQ=255)
                                                                        gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcg                                                                               <  2:1899782‑M2/150‑7 (MQ=255)
                                                                         cgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgt                                                                              <  1:1261820‑M2/150‑8 (MQ=255)
                                                                          gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCATTTTGGCCGctgtcgtt                                                                             <  1:6219179‑M2/150‑9 (MQ=255)
                                                                            ccTCGCGGTTTTTCTCTTTTTATTAAATTTTTCCGGTTTAT‑GGTGTTACCGTTTTTCTTCGTCTTAAATTATTGTTTTTATTTAAAATACCCCCTGAAAAGAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttc                                                                           >  2:23936‑M2/1‑140 (MQ=255)
                                                                            ccTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttc                                                                           >  2:2245803‑M2/1‑140 (MQ=255)
                                                                              tCGCGGGTTTTCGCTATTTATGAAAAATTTCCGGTTTAA‑GGTGTTTCCGTTCTTATTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttcct                                                                         <  1:5292288‑M2/150‑13 (MQ=255)
                                                                                gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttccttt                                                                       <  1:1062652‑M2/150‑15 (MQ=255)
                                                                                gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttccttt                                                                       <  1:193644‑M2/150‑15 (MQ=255)
                                                                                  gggTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttcctttct                                                                     >  2:6389380‑M2/1‑134 (MQ=255)
                                                                                     ttttCTTTATTTATTAAATTTTGCGGTTTTTT‑GTTGTTGCGTTTCTTCTTCGTCTAAACTTTATTTATTTATTAAAAAACCCCCCTCAAACCAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTg                                                                  <  1:1368141/150‑1 (MQ=9)
                                                                                       ttCGCTTTTTATTAAAATTTTCCTTTTTAT‑GTTGTTTCCTTTCTTCTTCTTCTTAAATTAATTTATTTATTAAAAAAACCCCCTCAAAATAATGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                <  1:517045/150‑1 (MQ=14)
                                                                                        tctttttttATTAAATTTTTCCTTTTTTT‑GGTGTTCCTGTTCTTCTCCTTTGTAATTTATTTTTTTTATTTAAAATACCCCCTTAAAACAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                               <  1:6721351/144‑1 (MQ=2)
                                                                                        tCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                               >  2:6281805/1‑150 (MQ=255)
                                                                                          gCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                             >  1:7602628/1‑150 (MQ=255)
                                                                                            tATTTTTTAAAATTTT‑CGGTTTAATGGTGTTTCCTTTCTTCTTCGTCTTAACTTATTTTTTTTATTAAAAAAACCCCATCAAAAGAATGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           <  1:9706256/150‑1 (MQ=11)
                                                                                              tttATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCATGAAAAGAAAGGAAACGACAGGTGATGAAAGCGAGCTTTTTGGCCTCTGTCGATACCTATCTATGTTTTTTTcc                                                         >  2:590270/1‑150 (MQ=14)
                                                                                               ttATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCg                                                        >  1:7301884/1‑150 (MQ=255)
                                                                                                tATTAAATTTTTCGGTTTTTT‑GGTTTTTCGTTTCTTCTTCGTCTTAAATTAATTTATTTATTAAAAAAACCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  1:1521095/150‑1 (MQ=14)
                                                                                                tATTAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  2:5344128/150‑1 (MQ=255)
                                                                                                tATGAAATTTTTCCGGTTTAT‑GGTTTTTCCTTTCTTGTTGGTGTTAACTTATTGTATTTATTTAAAATACCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  2:2572356/150‑1 (MQ=14)
                                                                                                  tGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTATAATACCCCCTGAAATGAAAGGAAACGACAGGTGCTGTAAGCGTGGTTTTTGGCATCTGTCGTTTCCTTTCTCTGTGTTTgtgggtgg                                                     >  2:6161829/1‑150 (MQ=14)
                                                                                                   gAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  2:3207982/1‑150 (MQ=255)
                                                                                                   gAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  2:791830/1‑150 (MQ=255)
                                                                                                    aaaaTTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATATAAAATACCCCCTTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACATTATATGTTTTTGTACGTGGaa                                                   >  1:2141272/1‑150 (MQ=14)
                                                                                                     aatttttCCGGTTTTT‑GGTGTTGCCGTTCTTCTTGGTGGTAACTTAATGTTTTTATTTAAAATACCCCCTCAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAt                                                  <  1:9451187/147‑1 (MQ=14)
                                                                                                     aaaTTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAt                                                  <  2:7281846/150‑1 (MQ=255)
                                                                                                        ttttCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  2:9515581/150‑1 (MQ=255)
                                                                                                         tttCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCGGTGGAATGTGc                                              >  2:4981922/1‑150 (MQ=255)
                                                                                                         tttCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTAATAAAAATAACACCTGAAAAGAAAGGAAACTAAAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTTTCCGTGGAAagtta                                              >  1:8619757/1‑148 (MQ=14)
                                                                                                           tCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                            >  2:4705932/1‑150 (MQ=255)
                                                                                                           tCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                            >  2:3670313/1‑150 (MQ=255)
                                                                                                             cGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                          <  2:5234971/150‑1 (MQ=255)
                                                                                                             cGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                          >  2:8879526/1‑150 (MQ=255)
                                                                                                              ggTTTAA‑GTTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         <  2:490209/150‑1 (MQ=255)
                                                                                                                tttAA‑GGTTTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  2:453502/150‑1 (MQ=255)
                                                                                                                tttAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  2:1382442/150‑1 (MQ=255)
                                                                                                                tttAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  2:6595954/1‑150 (MQ=255)
                                                                                                                   tatggTGTTTCCGTTCTTCTTCGTCGTAAATTAATGTATTTATTTAAAATACCCCCTTAAAACAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt                                     <  1:5526070/147‑1 (MQ=255)
                                                                                                                    a‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   >  2:8273321/1‑150 (MQ=255)
                                                                                                                       gtttttCCGTTCTTCTCCGTCTTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                 <  1:5990371/147‑1 (MQ=255)
                                                                                                                       gtgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                 >  1:1590539/1‑150 (MQ=255)
                                                                                                                       gtgtTTCCGTACTTCTTCGTCGAAACTTAATTTATTTATTTAAAATACCCCCTGAAAACAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                 <  1:6450249/150‑1 (MQ=255)
                                                                                                                        tgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTCCGTGGAATGTGCAATTGAATTcaaca                                >  1:5849029/1‑150 (MQ=255)
                                                                                                                          tttCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCGGTTTTTGTCCGTGGACTGTGCAATGGAAGTCAACaaa                              >  1:419100/1‑150 (MQ=255)
                                                                                                                              cGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagc                          >  2:1019640/1‑150 (MQ=255)
                                                                                                                              cGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACATTCTCAGTTTTTTTCCGTGGAAAGTGCAATGGAAGTCAACAAAAagc                          >  1:9008229/1‑150 (MQ=14)
                                                                                                                               gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         >  1:2712607/1‑150 (MQ=255)
                                                                                                                               gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         >  1:7117494/1‑150 (MQ=255)
                                                                                                                                    tcttcGTCGTAACTTAATTTATTTATTTAAAACACCCCCTGAAAAGAATGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctgg                    <  1:290727/150‑1 (MQ=14)
                                                                                                                                     cttcGTCTTAAATTATTGTTTTTTTTAAAAAAACCCCCTTAAAACAAAGCAAACGACTGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggc                   <  1:3160989/150‑1 (MQ=14)
                                                                                                                                      ttcGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  <  1:1413412/150‑1 (MQ=255)
                                                                                                                                      ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGACTGTGCAATGGAAGTCAACAAAAAGCAgctggct                  >  2:6761387/1‑150 (MQ=255)
                                                                                                                                      ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  >  2:5983898/1‑150 (MQ=255)
                                                                                                                                        cgtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                >  2:4012713/1‑150 (MQ=255)
                                                                                                                                                      gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  2:5399570/150‑1 (MQ=255)
                                                                                                                                                      gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  >  2:8062852/1‑150 (MQ=255)
                                                                                                                                                      gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  >  2:3509505/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAA‑GGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  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: ‑