Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,727 G→T 12.2% 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,7270GT12.2% 556.9 / 18.5 187intergenic (‑282/‑107)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (102/56);  new base T (2/27);  total (104/83)
Fisher's exact test for biased strand distribution p-value = 3.24e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

TGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580581‑580876
                                                                                                                                                   |                                                                                                                                                     
  tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                  >  2:3155308/1‑150 (MQ=255)
   gACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattta                                                                                                                                                 >  1:1410647/1‑150 (MQ=255)
   gACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattta                                                                                                                                                 <  2:3852704/150‑1 (MQ=255)
   gACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattta                                                                                                                                                 <  1:1951514/150‑1 (MQ=255)
    aCATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                                <  1:3046144/150‑1 (MQ=255)
      aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  1:3655975/1‑150 (MQ=255)
      aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  1:3759657/1‑150 (MQ=255)
      aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  1:2967713/1‑150 (MQ=255)
        gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            <  2:2749727/150‑1 (MQ=255)
         tAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                           >  2:3295004/1‑150 (MQ=255)
          aGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                          <  2:1368315/150‑1 (MQ=255)
            atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                                        >  2:1396269/1‑150 (MQ=255)
            atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                                        >  2:4032256/1‑150 (MQ=255)
             tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                       <  2:1455386/150‑1 (MQ=255)
               atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     >  2:2038723/1‑150 (MQ=255)
                 aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   <  2:2782647/150‑1 (MQ=255)
                  attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  >  2:2103051/1‑150 (MQ=255)
                   ttattaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTg                                                                                                                                 >  2:2307827/1‑150 (MQ=255)
                   ttattaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTg                                                                                                                                 >  2:514690/1‑150 (MQ=255)
                    tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                >  1:1787737/1‑150 (MQ=255)
                    tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                <  1:736598/150‑1 (MQ=255)
                       tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             <  2:1454237/150‑1 (MQ=255)
                       tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             >  1:3017083/1‑150 (MQ=255)
                         tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                           >  1:402942/1‑150 (MQ=255)
                          ctttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaa                                                                                                                          <  1:3946939/148‑1 (MQ=255)
                          cATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaa                                                                                                                          >  2:968559/1‑150 (MQ=255)
                           aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:3260890/150‑1 (MQ=255)
                           aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:2132776/150‑1 (MQ=255)
                           aTTTTGGGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:3167524/150‑1 (MQ=255)
                            ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTGCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAAGTTAATGTATTTATTTAAATTACCACCTGAaaagaaag                                                                                                                        >  1:2340387/1‑150 (MQ=255)
                            ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        >  1:308953/1‑150 (MQ=255)
                            ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        >  1:309215/1‑150 (MQ=255)
                            ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTACTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        >  1:4092867/1‑150 (MQ=255)
                             tttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAgg                                                                                                                       >  2:4142978/1‑150 (MQ=255)
                                gtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                    <  1:1271712/150‑1 (MQ=40)
                                     tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACa                                                                                                               <  2:3411911/150‑1 (MQ=39)
                                       cTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAgg                                                                                                             >  2:309399/1‑150 (MQ=39)
                                        tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  1:3664346/1‑150 (MQ=39)
                                         ttCCGGCGATCCGACTGGTTCGGGGGGGGGGACGTCGGGGGTTTGCGCTTTTTATGAAATTTTTCGGG‑TTTACGGTGTTTCCGTTCTTCTTGTTCTTACCTTTTTGTTTTTATTTAAAATACCCCCTGAAACGAACGGACACGACAGGTg                                                                                                           <  2:909676/150‑1 (MQ=2)
                                              gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaa                                                                                                      >  1:685819/1‑150 (MQ=38)
                                              gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaa                                                                                                      >  1:2517607/1‑150 (MQ=38)
                                                 aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  2:1656263/1‑150 (MQ=38)
                                                 aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  2:3717870/1‑150 (MQ=38)
                                                 aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  1:3199792/1‑150 (MQ=38)
                                                    cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                                >  2:2853964/1‑150 (MQ=39)
                                                     gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:1198109/1‑150 (MQ=39)
                                                      aCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTCATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGATGAAAGAGAGCt                                                                                              >  1:2820601/1‑150 (MQ=255)
                                                        aGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            <  2:2699062/150‑1 (MQ=39)
                                                           ttACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                                         <  2:3376569/150‑1 (MQ=40)
                                                             aCGGGGCGGCGACCTCGCGGGTTTTCGCTAGTTATGAAAATTTTCCGG‑GTTAAGGTGTTTACGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCGTTTTTGc                                                                                       >  1:1059845/1‑150 (MQ=14)
                                                               ggggCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTGTTTGGCct                                                                                     >  2:3618626/1‑150 (MQ=255)
                                                                       gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATAACCCCAGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCACTGTCGtt                                                                             >  1:3708628/1‑150 (MQ=255)
                                                                           tCGGGGGTTTGCGATATTTATGAAATTTTTCGGT‑TTTATGGTGTTGCCGTTCTTCTGCGTCGTAACTTATTGTATTTATTTAAAATACCCCCTGAAAATAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         <  2:1857908/150‑1 (MQ=14)
                                                                           tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         >  1:3498039/1‑150 (MQ=255)
                                                                           tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         >  2:2845300/1‑150 (MQ=255)
                                                                             gcgGGTTTCCGCTATTTCTGAAAATTTTCCGG‑TTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTACTGTATTTATTTAAAATACCCCCTGAACAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       <  2:3350347/150‑1 (MQ=14)
                                                                               gggTTTTCGCTATTTATGAAAATTTTCCGG‑TTTATGGTGTTTCCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  2:56943/150‑1 (MQ=255)
                                                                               gggTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     >  2:125178/1‑150 (MQ=255)
                                                                               gggTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     >  2:4016233/1‑150 (MQ=255)
                                                                               gggTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     >  2:1962236/1‑150 (MQ=255)
                                                                                    ttCTCTATTTATTAAATTTTTCGGT‑TTTATGGTGTTGCCTTTCTTCTTGTCCTTAAATTAATTTATTAATTAAAAAAACCCCCTTAAAAGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                <  1:3385678/150‑1 (MQ=14)
                                                                                    ttCGCTATTAATGAAAATTATCCGG‑TTTAAGGTGTTACCGTTATACATCGTCGTAACTTAATGTATTTATTTAAAATACAACCTGAAAATAAATGAAACGAAAGGAGCTGAAAGAGAGCTATTTGGCCTCTGTCGTTTCCTGTCTATGtt                                                                >  2:1065504/1‑150 (MQ=255)
                                                                                    ttCGATTTTTATTAAATTTTTCGGG‑TTTCTGGTGTTACAGTTCTCCTCCGTCTTAATTTATTTTATTTATTTAAAATACCCCCTCAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                <  1:2098269/150‑1 (MQ=14)
                                                                                     tCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                               >  1:1652154/1‑150 (MQ=255)
                                                                                       gCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                             >  1:2263037/1‑150 (MQ=255)
                                                                                           tttATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTATTTTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                                                         <  2:1625688/150‑1 (MQ=255)
                                                                                            ttATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCg                                                        >  1:3473465/1‑150 (MQ=255)
                                                                                                ttaaaTTTTGCGTTTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTTTTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTgg                                                     <  1:3295634/148‑1 (MQ=14)
                                                                                                gAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  1:756866/1‑150 (MQ=255)
                                                                                                    aTTTTCCGG‑TTTATGGTGTTCCCTTTGTTCTTCGTCTTAAATTAATTTATTTATTTAAAACACCCCCTGAAAACAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgt                                                <  1:3004207/150‑1 (MQ=14)
                                                                                                     ttttCCTT‑TTTATGGTGTTTCCGTTCTTCTTCTTCTTAAATTAATTTTTTTATTTAAAATACCCCCTAAAAATAAATGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  1:2980863/150‑1 (MQ=14)
                                                                                                     ttttCCGT‑TTTAAGGTTTTTCCGTTTTTGTTCTTCGTAAATTAATTTATTTATTTAAAATACCCCCTCAAAAGAAAGTACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  1:2612347/150‑1 (MQ=14)
                                                                                                     tttcCCGT‑TTTATTGTTTTTCTTTTCTTCTTGGTATTAAATTAATTTATTTATTTAAAATACCCCCTTAAAATAAACCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                               <  1:1060503/150‑1 (MQ=14)
                                                                                                           gg‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         >  1:2203317/1‑150 (MQ=255)
                                                                                                              ttttcggTGTTTCCTTTCTTCTTCTTCTTAAATTATTTTTTTTTTTTAAAATCCCCCCTTAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:1247009/145‑1 (MQ=14)
                                                                                                              tttATGGTTTTTCCTTTCTTCTCCTTGTTATTTTATTTAATTTTTTAAAAACACCCCTTTAAATTAATGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:2794234/150‑1 (MQ=9)
                                                                                                              tttATGGTGTTACCGTTCTTCTTCGTCTAAAATTATTTTATTTATTTAAAATACCCCCTTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:968678/150‑1 (MQ=17)
                                                                                                              tttAAGGTTTTTCCGTTCTTCTTCGTCTTAAATTATTGTATTTATTTAAAATACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  2:2204555/150‑1 (MQ=14)
                                                                                                              tttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTCCGTGGAATGTGCAATGGaa                                       >  1:198826/1‑150 (MQ=255)
                                                                                                              tttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  1:917325/1‑150 (MQ=255)
                                                                                                               ttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAg                                      >  2:3487765/1‑150 (MQ=255)
                                                                                                                tatggTGTTCCGGTTCTTCGTCGTCTTAACTTATTTTTTTTATTTAAAAAACCCCCTTAAAATAAAGTACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt                                     <  1:1864066/147‑1 (MQ=14)
                                                                                                                 aaGGTGTTTCGGTTATTCTTGGTCGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTc                                    <  2:740308/150‑1 (MQ=17)
                                                                                                                        ttCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             <  2:2558903/150‑1 (MQ=255)
                                                                                                                         tCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                            >  1:533183/1‑150 (MQ=255)
                                                                                                                            gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         <  2:1284356/150‑1 (MQ=255)
                                                                                                                                 tcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctgg                    >  1:1473387/1‑150 (MQ=255)
                                                                                                                                  cttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggc                   >  2:542021/1‑150 (MQ=255)
                                                                                                                                   ttcttcTTAATTTAATTTATTTTTTTAAAAAACCCCCTCAAAATAATGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  <  1:2145253/150‑1 (MQ=14)
                                                                                                                                    tcttcgcAACTCAATGTATTAATTAAAAAACCCCCCTTAAACCAACGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                 <  1:3586968/147‑1 (MQ=14)
                                                                                                                                     cgtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAACAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                <  2:3708628/150‑1 (MQ=255)
                                                                                                                                       tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa              <  2:4041344/150‑1 (MQ=255)
                                                                                                                                       tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa              >  2:24346/1‑150 (MQ=255)
                                                                                                                                       tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTACAATGGAAGTCAACAAAAAGCAGCTGGCTGACa              >  2:1367127/1‑150 (MQ=255)
                                                                                                                                       tcgtAACTTAATGTATTTATTTAAAATACCCCATGAAAAGAAAGGAAACGACAGCTGCTGAAAGCCAGCTTTTTGTCCTCTCTCGTTTCCTGTCTCTTTTTTTGTCCGTGGACTGTTCAAAGGAAGACAACAAAAAGCAGCTCGCTGACa              >  2:2809945/1‑150 (MQ=255)
                                                                                                                                         gtAACTTAATGTATTTATTTAAAATACCCACTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtt            >  2:4083304/1‑150 (MQ=255)
                                                                                                                                            agttAATGTATTTATTAAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTc         <  1:567183/148‑1 (MQ=255)
                                                                                                                                            aCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTc         >  2:2913676/1‑150 (MQ=255)
                                                                                                                                              ttAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGGGGACTGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCgg       >  1:1928186/1‑150 (MQ=255)
                                                                                                                                              ttAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGATTTTTGGCATCTGTCGTTTCCTTTCTCTGTTGTTTTCCGTGGAATGTGCAAGGGAAGTCCACAAAAAGCAGCTGGCGGACATTTTCgg       >  1:2142188/1‑150 (MQ=14)
                                                                                                                                                   |                                                                                                                                                     
TGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580581‑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: ‑