Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,727 G→T 7.1% 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 7.1% 2490.5 / 40.3 792intergenic (‑282/‑107)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (418/291);  new base T (5/78);  total (423/369)
Fisher's exact test for biased strand distribution p-value = 2.64e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 808 total aligned reads displayed.

CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580578‑580876
                                                                                                                                                      |                                                                                                                                                     
     tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                  >  2:11080757/1‑150 (MQ=255)
      gACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGtattta                                                                                                                                                 <  1:10981418/150‑1 (MQ=255)
           gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            >  1:7590048/1‑150 (MQ=255)
              gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                                         >  2:4148876/1‑150 (MQ=255)
               atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                                        >  1:10288611/1‑150 (MQ=255)
                 gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                      >  2:14656625/1‑150 (MQ=255)
                  atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     <  1:5618390/150‑1 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   >  1:4317829/1‑150 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   >  2:7031805/1‑150 (MQ=255)
                    aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   <  1:8102091/150‑1 (MQ=255)
                     attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  <  1:5795014/150‑1 (MQ=255)
                        attatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTAGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaa                                                                                                                               <  1:13522045/150‑1 (MQ=255)
                          tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             <  1:1523195/150‑1 (MQ=255)
                          tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             <  2:15249478/150‑1 (MQ=255)
                          tatCATTTTGTGGGCCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                             >  1:6147260/1‑150 (MQ=255)
                            tCATTTTGTGTGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCTGAaaaga                                                                                                                           >  2:1729144/1‑150 (MQ=255)
                             cATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaa                                                                                                                          <  2:9438978/150‑1 (MQ=255)
                             cATTTTGGGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaatgaa                                                                                                                          <  1:6323169/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:13158718/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:9709966/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:4666917/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:4053281/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:10798925/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:12374307/150‑1 (MQ=255)
                              aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:9146107/150‑1 (MQ=255)
                               ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        >  2:2108318/1‑150 (MQ=255)
                               ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTAATTAAAATACCCCCTGAaaagaaag                                                                                                                        >  1:2893805/1‑150 (MQ=255)
                                   gtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                    <  2:10299984/150‑1 (MQ=40)
                                      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 >  2:11027404/1‑150 (MQ=39)
                                      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 >  1:5027969/1‑150 (MQ=39)
                                        tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                          <  2:12649499/123‑1 (MQ=40)
                                          cTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAgg                                                                                                             >  1:14453454/1‑150 (MQ=39)
                                                ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                       >  2:777456/1‑150 (MQ=39)
                                                  cGATCCGACTGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGTTGTTT‑CCGTTATTCTTCGTCGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGaaa                                                                                                     <  2:4707146/150‑1 (MQ=11)
                                                   gATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAg                                                                                                    >  2:8892288/1‑150 (MQ=38)
                                                    aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGACACGACAGGTGCTGAAAGc                                                                                                   >  2:2937827/1‑150 (MQ=38)
                                                      ccGACAGTTTACGGGGCGGCGACCTCGCGGGTTTTCTCTATTTATGAAAATTTTCCGGTTTAAGTTGTTT‑CCGTTCTTCTTCGTCGTAAATTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGa                                                                                                 <  2:10065166/150‑1 (MQ=255)
                                                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         >  1:11189189/1‑124 (MQ=38)
                                                           aGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTCTTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            >  2:11807430/1‑150 (MQ=255)
                                                             gTTCCGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGTTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:5248883/150‑1 (MQ=255)
                                                             gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:5006279/150‑1 (MQ=40)
                                                              ttACGGGGCGGCGACCTCGCGGGTTTTCGCTCTTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                                         >  2:200019/1‑150 (MQ=255)
                                                                     gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTg                                                                                  >  1:13817870/1‑150 (MQ=255)
                                                                     gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTC‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTg                                                                                  <  2:6559564/150‑1 (MQ=255)
                                                                       ggcgACCTCGCGGTTTTTCGCTATTTATGAAATTTTACCGGTTTAAGTTGTTT‑CCGTTCTTCTTCGTCTTAAATTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTc                                                                                <  2:7988015/150‑1 (MQ=14)
                                                                       ggcgACCTCGCGGTTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTc                                                                                <  2:6063336/150‑1 (MQ=255)
                                                                         cgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                              >  2:8194241/1‑150 (MQ=255)
                                                                          gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACGGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                             >  2:3173155/1‑150 (MQ=255)
                                                                          gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                             >  2:3662331/1‑150 (MQ=255)
                                                                            ccTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAATGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTc                                                                           <  1:14042627/150‑1 (MQ=255)
                                                                             cTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTAGtttac                                                                          >  1:12654353/1‑148 (MQ=255)
                                                                              tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         >  1:3152135/1‑150 (MQ=255)
                                                                              tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         <  2:14940355/150‑1 (MQ=255)
                                                                              tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGCTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         >  2:12363372/1‑150 (MQ=255)
                                                                                  ggtttttCTCTATTTATCAAAATTGTCC‑GTTTAATGTTTTTCCCTTTCTTAGGGGTCTTAACTTATTTTTTTTATTTAAAAAACCCCCTTAAAACAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  1:2427478/147‑1 (MQ=2)
                                                                                  gggTTTTCGCTATTTATGAAAATTTGCCGGTTTAAGGGGTTA‑CCTTACTCCTCCGTCGTAAATTAATGTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  2:13632/150‑1 (MQ=14)
                                                                                  gggTTGTCTCTTTTTATGAAAATTTGCCGGTTTATGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAAGGTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  2:12085173/150‑1 (MQ=14)
                                                                                   ggTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACATAATGTATTTATTAAAAATACCCCCTGAAAACAAAGGAAACGACAGGAGCTGAAAGCGACCTTTATGGCCTCTGTCGTTACCATtcac                                                                    >  1:7957423/1‑148 (MQ=14)
                                                                                       ttctatATTTATGAAAATTTTCCTTTTTAAGGTGTTA‑CCTTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                <  2:13364216/145‑1 (MQ=14)
                                                                                       ttCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                >  2:9895502/1‑150 (MQ=255)
                                                                                         cGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGtttt                                                              >  1:1475524/1‑150 (MQ=255)
                                                                                            tATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         >  2:6762746/1‑88 (MQ=255)
                                                                                              tttATTAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTGTTCGTCGTGGTAACTTAATTTTTTTATTTAAAATACCCCCTCAAAACAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                                                         <  2:5713580/150‑1 (MQ=14)
                                                                                              tttATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                                                         >  2:5651719/1‑150 (MQ=255)
                                                                                                ttttaaatTTTTCCGTTTTATGGTGTTT‑CCTTTCTTTTTCGTCTTAAATTAATTTATTTATTTAAAAAACCCCCTCAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  1:10019773/146‑1 (MQ=14)
                                                                                                tATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       >  2:3947924/1‑150 (MQ=255)
                                                                                                tATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  2:3408274/150‑1 (MQ=255)
                                                                                                  tGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTgg                                                     <  2:1144062/150‑1 (MQ=255)
                                                                                                   gAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  1:1063458/1‑150 (MQ=255)
                                                                                                    aaaaTTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   >  1:7260978/1‑150 (MQ=255)
                                                                                                    aaaaTTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAAATTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   <  2:5908509/150‑1 (MQ=255)
                                                                                                     aaaTTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGTCCTCTGTCGTTTCCTTTCTCTGTTTTTGTAAGTGGAAt                                                  >  1:3783138/1‑150 (MQ=255)
                                                                                                         tttCCGGTTTATGGTGTTT‑CCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCATGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGc                                              <  2:14141229/150‑1 (MQ=255)
                                                                                                              ggTTTATGGTTTTT‑CCTTTCTTCTCCTTCTTAAATTAATTTATTTATTAAAAATACCCCCTCAAAATAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         <  1:12945892/150‑1 (MQ=14)
                                                                                                                tttATGGTTTTC‑CCTTTCTTCTTCGTCTTAACTTATTGTATTTATTTAAAATACCCCCTTAAAACAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:2236663/150‑1 (MQ=14)
                                                                                                                tttATGGTGTTT‑CGGTTCTTCTTGGTCGTAACTTAATTTATTTATTTAAAACACCCCCTTAAAATAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:6989508/150‑1 (MQ=17)
                                                                                                                tttATGGTGTTT‑CCTTTCTCCTTCGTCGTAACTTAATTTATTTATTTAAAATACCCCCTCAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:8500397/150‑1 (MQ=17)
                                                                                                                tttAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  1:15483432/1‑150 (MQ=255)
                                                                                                                tttAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       >  1:5197515/1‑150 (MQ=255)
                                                                                                                 ttAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAg                                      <  2:15475209/150‑1 (MQ=255)
                                                                                                                    tggTTTTT‑CTTTTGTTCTTCGTCTAAACTTAATTTATTTATTTAAAAAACCCCATGAAAATAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   <  1:6902640/149‑1 (MQ=14)
                                                                                                                     ggTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaa                                  >  2:2421948/1‑150 (MQ=255)
                                                                                                                      gtgtTC‑CGGTTCTTCTTCGTCGTAACTTATTTTATTTATTTAAAATACCCCCTCAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                 <  1:1700771/150‑1 (MQ=14)
                                                                                                                         ttt‑CCGTTCTTCTTCGTCTTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaa                              <  2:11240291/150‑1 (MQ=255)
                                                                                                                          tt‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             >  2:14539783/1‑150 (MQ=255)
                                                                                                                           ttacGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             <  2:13910566/147‑1 (MQ=255)
                                                                                                                           t‑CCGTTCTTCTTCTACTTAAATTAATGTATTTATTAAAAAAACCCCCTGAAAACAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                            <  1:2537473/150‑1 (MQ=14)
                                                                                                                           t‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGGAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaaa                            >  2:7341475/1‑150 (MQ=255)
                                                                                                                               gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGTCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         >  2:10471098/1‑150 (MQ=255)
                                                                                                                               gTTCTTCTTCGACGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         <  2:10761243/150‑1 (MQ=255)
                                                                                                                                 tcttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc                       >  1:6626606/1‑150 (MQ=255)
                                                                                                                                  cttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAGAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                      >  1:9132923/1‑150 (MQ=255)
                                                                                                                                    tcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctgg                    <  1:9891955/150‑1 (MQ=255)
                                                                                                                                       tcgtcgTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                 >  2:12855754/1‑150 (MQ=255)
                                                                                                                                          tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa              >  1:12044705/1‑150 (MQ=255)
                                                                                                                                                cTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg        >  1:3898044/1‑150 (MQ=255)
                                                                                                                                                   aaTTTATTTATTTAAAATACCCCCTTAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg     <  1:11502325/150‑1 (MQ=255)
                                                                                                                                                      gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  >  2:2602236/1‑150 (MQ=255)
                                                                                                                                                      gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  2:8564461/150‑1 (MQ=255)
                                                                                                                                                      gTATTTATTTAAAATACCCCCTCAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  1:3464879/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTT‑CCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  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: ‑