Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,717 G→T 6.0% intergenic (‑272/‑117) 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,7170GT 6.0% 1004.2 / 7.8 306intergenic (‑272/‑117)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (167/110);  new base T (2/26);  total (169/137)
Fisher's exact test for biased strand distribution p-value = 4.79e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

CGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAT  >  NC_000913/580568‑580865
                                                                                                                                                      |                                                                                                                                                    
         aCGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAt                                                                                                                                             <  2:896582/150‑1 (MQ=255)
         aCGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAt                                                                                                                                             <  2:598819/150‑1 (MQ=255)
         aCGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAt                                                                                                                                             <  2:3454605/150‑1 (MQ=255)
         aCGGTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTCTGAAAATTTTCCGGTTTAAGGTGTTTCCGGTCTTCTTCGTCGTAACTTAAt                                                                                                                                             >  1:3085299/1‑150 (MQ=255)
           gTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGt                                                                                                                                           <  1:1069149/150‑1 (MQ=255)
               tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                       >  2:5781762/1‑150 (MQ=255)
                 aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAt                                                                                                                                             <  1:3207782/142‑1 (MQ=255)
                 aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                     >  1:2939094/1‑150 (MQ=255)
                   aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                   >  1:5617710/1‑150 (MQ=255)
                    tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTAtttattta                                                                                                                                  <  1:396979/150‑1 (MQ=255)
                     gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTCTTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                 >  2:4401249/1‑150 (MQ=255)
                      tAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                >  2:683274/1‑150 (MQ=255)
                        gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                              >  2:293502/1‑150 (MQ=255)
                         atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                             >  1:6146131/1‑150 (MQ=255)
                          tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                            >  2:6295286/1‑150 (MQ=255)
                          tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                            <  1:2885666/150‑1 (MQ=255)
                           gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                           >  1:792950/1‑150 (MQ=255)
                            atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                          <  1:4897284/150‑1 (MQ=255)
                            atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                          <  2:5986251/150‑1 (MQ=255)
                             tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcccc                                                                                                                         <  2:2750926/150‑1 (MQ=255)
                              aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcaccc                                                                                                                        >  2:1286710/1‑150 (MQ=255)
                               attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                       >  2:1962380/1‑150 (MQ=255)
                                ttattaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTg                                                                                                                      >  1:4134396/1‑150 (MQ=255)
                                  attatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCTGaa                                                                                                                    >  1:5265747/1‑150 (MQ=255)
                                    tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                  <  2:1290200/150‑1 (MQ=255)
                                      tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                >  1:1391959/1‑150 (MQ=255)
                                      tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                >  1:5806974/1‑150 (MQ=255)
                                       cATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaa                                                                                                               <  2:4771670/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  2:2338301/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  2:2685903/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  2:4920946/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  2:5369907/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  1:103090/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  1:1765493/150‑1 (MQ=255)
                                        aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                              <  1:4850795/150‑1 (MQ=255)
                                         ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGGTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                             >  2:1471642/1‑150 (MQ=255)
                                              tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAc                                                                                                        <  1:5369867/150‑1 (MQ=40)
                                              tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAc                                                                                                        >  1:667020/1‑150 (MQ=40)
                                                ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                      >  1:2754737/1‑150 (MQ=39)
                                                ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTGTAAGGTGTTTCCGTTCTTCTTCGTCGTGACTTAGTGTGTTTAGGTAATTTACCCCCTGTCGTGAAAGGAGACGa                                                                                                      >  1:3807342/1‑150 (MQ=11)
                                                 gTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAc                                                                                                     >  2:163463/1‑150 (MQ=39)
                                                        ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                              <  2:6064556/150‑1 (MQ=39)
                                                            ctaTCCGACTGGTTAGGGGGGGGCGACCTCGCGGGTTTT‑CGATATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                          <  2:5668257/148‑1 (MQ=255)
                                                            cGATCCCACAGGTTAGGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATTAAATTTTTCGGGTTTATGGTGTTGCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                          <  2:5889719/150‑1 (MQ=11)
                                                              aTCCGACTGGTTAGGGGGCGGCGACCTCGGGGGTTTT‑CGCTATTTATTAAATTTTTCCGGTTTACGGTGTTTCCGTTCTTCTTGGTCGTAACTTAATGTTTTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGc                                                                                        <  2:1913179/150‑1 (MQ=11)
                                                                ccGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGa                                                                                      <  2:3805375/150‑1 (MQ=39)
                                                                 cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                     >  2:4496408/1‑150 (MQ=39)
                                                                 cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                     >  1:2956087/1‑150 (MQ=39)
                                                                  gACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                    >  1:4222351/1‑150 (MQ=39)
                                                                       gTTACGGGGCGGCGACCTCGCGGTTTTT‑CTCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                               <  2:6514791/150‑1 (MQ=255)
                                                                       gTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACCACAGGTGCTTAAAGCGAGCttttt                                                                               >  1:4644214/1‑150 (MQ=255)
                                                                        ttAGGGGGCGGCGACCTCGGGGTTTTT‑CGCTATTTATGAAATTTTTCCGTTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTTTTGTTTTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                              <  2:2499108/150‑1 (MQ=11)
                                                                        ttACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                              >  1:2229297/1‑150 (MQ=40)
                                                                         tagggggCGGCGACCTCGGGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTACGGTGTTTCCGTTCTTCTTCGTCGTAAATTAATTTATTTATTTAAAATACCCCCTGAAAAGAACGTAAACGACAGGTGCTGAAAGCGAGCTTTTTgg                                                                             <  2:1821112/147‑1 (MQ=14)
                                                                          aCGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGc                                                                            <  2:459225/150‑1 (MQ=255)
                                                                           cGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGcc                                                                           >  1:1479961/1‑150 (MQ=255)
                                                                            ggggCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                          >  2:2673074/1‑150 (MQ=255)
                                                                            ggggCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAAATTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                          >  1:3888661/1‑150 (MQ=255)
                                                                             gggCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAATGAAACGACAGGTGCTGAAAGCGATCTTTTTGGCctc                                                                         >  1:2112391/1‑150 (MQ=255)
                                                                                       cctgggggTTTTCCTCTATTTATTAAATTTTTCCGGTTTCTGGTGTTTCGGTTCTATTTCGCCTTAACTTAATTTATTTATTAAAAATACCCCCTCAAAACAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTc                                                                <  1:6010298/145‑1 (MQ=255)
                                                                                        tCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                              >  1:6243727/1‑150 (MQ=255)
                                                                                        gcgcgGGTTTT‑CGGTTTTTCTGAAACTTTTCCGGTTTCAGGTGTGTCCGTTCTTGTTCGTCGTAACTTAATGTATTTATTTAACATACCCCCTGAACAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                              <  2:3755248/149‑1 (MQ=14)
                                                                                            gggTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                          >  1:2478963/1‑150 (MQ=255)
                                                                                               tttt‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTGTTTGGCCTCTGTCGTTTCCTTTCTCTg                                                       >  2:255687/1‑150 (MQ=255)
                                                                                                ttt‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTCCGTCGTAACTTAATGTATTTTTTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGt                                                      <  2:755906/150‑1 (MQ=255)
                                                                                                  t‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                    <  2:6568891/150‑1 (MQ=255)
                                                                                                   tcgttttttATTAATTTTTCCCGTTTTATGGTGTTTCGTTTCTTCTTGGTCTTAAATTATTTTTTTTATTTAAAACCCCCCCTCAAAATAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                    <  1:5036759/144‑1 (MQ=2)
                                                                                                     gCTATTTTTTAAAATTTTCCGGTTTATGGTTTTTCCGTTCTTATTCGTAGTAACTTATTTTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                  <  2:6240802/150‑1 (MQ=14)
                                                                                                      cTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTg                                                 >  1:5852384/1‑150 (MQ=255)
                                                                                                       tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                           >  1:1969732/1‑75 (MQ=255)
                                                                                                       tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                >  2:5497978/1‑150 (MQ=255)
                                                                                                       tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                >  2:848291/1‑150 (MQ=255)
                                                                                                           tATTAAAATTTGCTGTTTTTTGGTGTTTCCTTTCTTCTTCGTCTTAAATTAATTTTTTTATTAAAAATACCCCCTCAAAACAATGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                            <  1:6492536/150‑1 (MQ=14)
                                                                                                           tATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCGGAACAGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                            <  2:4134332/150‑1 (MQ=255)
                                                                                                              gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                         >  1:4746831/1‑150 (MQ=255)
                                                                                                              gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTAATTTCTCTGTTTTTGTCCGTGGa                                         >  1:4086126/1‑150 (MQ=255)
                                                                                                              gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAAGTATTTATATAAAATACCCCCTGAAAAGAAATGAAACGACAGGTGCTGAAAGCGAGCTTTTTTGACTCTGTCGTATAATTTATCTGTTTTTGTCCTTGGa                                         >  1:3315631/1‑150 (MQ=14)
                                                                                                               aaaaTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                        >  1:4393810/1‑150 (MQ=255)
                                                                                                                   ttttCCGGTTTACGGTGTTTCGGTTCTTCTTCTTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                    <  2:3962942/150‑1 (MQ=255)
                                                                                                                   ttttCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                    >  2:3915896/1‑150 (MQ=255)
                                                                                                                       cgggtttttggtttttCCGTTCTTCTTCGTCTTAAATTATTTTTTTTTTTTAAAATACCCCCTTAAATTAATGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAAt                                <  1:3697427/148‑1 (MQ=14)
                                                                                                                        ggttttATGGTGTTTCCGTTCTTCTTCGTCTTAACTTATTTTATTTATTTAAAATACCCCCTTAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                               <  2:3073960/147‑1 (MQ=17)
                                                                                                                        cgttttATGGTTTTTCCTTTCTTCTTCGTCTTAATTTATTTTATTTATTTAAAATACCCCCTGAAACTAAATTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                               <  2:1367383/147‑1 (MQ=14)
                                                                                                                        cGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                               >  1:1255389/1‑150 (MQ=255)
                                                                                                                           ttttAGGTTTTTCGTTTCTTCTTCTTCGTAAATATATGTTTTTATTTAAAAAACCCCCTTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                            <  1:771116/150‑1 (MQ=14)
                                                                                                                           tttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                            >  2:115848/1‑150 (MQ=255)
                                                                                                                            ttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCAGTCGTTTCCTTTCTCTGATTTTGTACGTGGAATGTGCAATGGAAg                           >  2:4009299/1‑150 (MQ=255)
                                                                                                                              aaGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTc                         >  2:5063360/1‑150 (MQ=255)
                                                                                                                                ggTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaa                       <  2:2164128/150‑1 (MQ=255)
                                                                                                                                 gtgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                      >  2:3328696/1‑150 (MQ=255)
                                                                                                                                     ttCCGTTGTTCTTCGTCGTAACTTATTGTTTTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                  <  1:942287/150‑1 (MQ=255)
                                                                                                                                      tCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                 >  2:180032/1‑150 (MQ=255)
                                                                                                                                      tCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                 >  1:3525581/1‑150 (MQ=255)
                                                                                                                                         gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca              >  1:2139206/1‑150 (MQ=255)
                                                                                                                                            cttcttcGTCTTAAATTATTTTATTTATTTAAAATACCCCCTTAAAACAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct           <  1:132622/150‑1 (MQ=14)
                                                                                                                                            cttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct           >  2:5412326/1‑150 (MQ=255)
                                                                                                                                            cttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct           >  2:4286080/1‑150 (MQ=255)
                                                                                                                                             ttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctg          >  1:2664294/1‑150 (MQ=255)
                                                                                                                                                ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct       >  2:5102328/1‑150 (MQ=255)
                                                                                                                                                    tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa   >  1:5216498/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTT‑CGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAT  >  NC_000913/580568‑580865

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: ‑