Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 576,192 T→A 10.1% intergenic (‑291/‑98) 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
*CP009273576,1920TA10.1% 629.1 / 21.0 237intergenic (‑291/‑98)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base T (131/79);  new base A (3/24);  total (134/103)
Fisher's exact test for biased strand distribution p-value = 3.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

ATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGT  >  CP009273/576043‑576340
                                                                                                                                                      |                                                                                                                                                    
aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATATATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTAGTCGTAACTTAATGTatttattt                                                                                                                                                      >  2:1036243/1‑150 (MQ=255)
 tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTAtttattta                                                                                                                                                     >  1:1439666/1‑150 (MQ=255)
  gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                                    <  1:149893/150‑1 (MQ=255)
  gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                                    <  1:150149/150‑1 (MQ=255)
    aGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                                  >  2:4924878/1‑150 (MQ=255)
        gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                              <  2:4146708/150‑1 (MQ=255)
         atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                             <  2:4835629/150‑1 (MQ=255)
         atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                             >  2:4683945/1‑150 (MQ=255)
         atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                             <  2:1393550/150‑1 (MQ=255)
         atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                             >  1:260117/1‑150 (MQ=255)
         atAATTATGATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGATTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                             <  2:2540469/150‑1 (MQ=255)
           aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                           <  2:1903463/150‑1 (MQ=255)
           aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                           >  1:589995/1‑150 (MQ=255)
            attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCTCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                          >  1:2393973/1‑150 (MQ=255)
            attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                          >  1:1078976/1‑150 (MQ=255)
            attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                          <  2:1963122/150‑1 (MQ=255)
             ttatgaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTg                                                                                                                                         <  2:5105285/150‑1 (MQ=255)
              tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                        <  1:2770575/150‑1 (MQ=255)
              tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                        <  1:3176576/150‑1 (MQ=255)
               attatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCTGaa                                                                                                                                       >  1:2465988/1‑150 (MQ=255)
                 tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaaa                                                                                                                                     <  1:4389744/150‑1 (MQ=255)
                  atCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaag                                                                                                                                    >  1:3871823/1‑150 (MQ=255)
                  atCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAaaag                                                                                                                                    >  1:3706930/1‑150 (MQ=255)
                   tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                                   >  2:366470/1‑150 (MQ=255)
                   tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                                   >  2:1955311/1‑150 (MQ=255)
                     aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCTCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                                 <  2:290253/150‑1 (MQ=255)
                     aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                                 >  1:3798842/1‑150 (MQ=255)
                      ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                                >  2:4240048/1‑150 (MQ=255)
                       tttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAgg                                                                                                                               <  1:4287166/150‑1 (MQ=255)
                           tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAc                                                                                                                           >  1:2628397/1‑150 (MQ=40)
                           tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAc                                                                                                                           <  2:2445063/150‑1 (MQ=40)
                           tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAc                                                                                                                           >  1:2244376/1‑150 (MQ=40)
                            gggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACg                                                                                                                          <  1:4660435/150‑1 (MQ=39)
                            gggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACg                                                                                                                          >  1:4768965/1‑150 (MQ=39)
                             ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                         >  1:4973250/1‑150 (MQ=39)
                             ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                         >  2:3283229/1‑150 (MQ=39)
                              gTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAc                                                                                                                        >  2:1789800/1‑150 (MQ=39)
                                 tttgccGGCGATCCGATTTGTTTGGGGGCGGCGACCTCGCGGTTTTTCGCTATTTATGAAAATTTTCGGTTTTATGGTGTTTCCGTTATTTTTGGTCTTAAATTAATGTATTTATTTAAAATACCCCCTTAAAAGAAAGCACACGACAgg                                                                                                                     <  2:4649078/146‑1 (MQ=14)
                                     ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTCATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                                                 >  1:2969098/1‑150 (MQ=255)
                                       ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                               >  1:221878/1‑150 (MQ=39)
                                       ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                               >  1:2743850/1‑150 (MQ=39)
                                         cGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                                             >  1:3689495/1‑150 (MQ=38)
                                          gATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAg                                                                                                            >  1:3846858/1‑150 (MQ=38)
                                          gATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAg                                                                                                            >  2:3591376/1‑150 (MQ=38)
                                            tCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCg                                                                                                          <  2:3562824/150‑1 (MQ=255)
                                              cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                                        >  1:1773254/1‑150 (MQ=39)
                                               gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                                       >  2:571724/1‑150 (MQ=39)
                                                   ggTTACGGGGCGGCGACCTCGCGGGTTTTCGCTCTTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTATTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAAGGAAACGACAGGTGCTGTAAGCGAGCtttt                                                                                                   >  1:1932361/1‑150 (MQ=255)
                                                    gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGTAAGCGAGCttttt                                                                                                  >  2:2332714/1‑150 (MQ=255)
                                                           ggcggcGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctct                                                                                           >  2:5168297/1‑150 (MQ=255)
                                                                 gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                                     >  2:312323/1‑150 (MQ=255)
                                                                 gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTACGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                                     <  2:4875068/150‑1 (MQ=255)
                                                                      cgcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCtt                                                                                >  2:4307520/1‑150 (MQ=255)
                                                                         ggtttttCTCTATTTATTAAATTTTTCGTTTTTTTGGTGTTTCGTTTCTTCTTCGACTAAAATTATTTTTTTTATTAAAAAACCCCCCACAAAACACCCCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                             <  1:51454/147‑1 (MQ=2)
                                                                         ggtttttCTCTATTTATTAAAATTTTCCGGTTTAAGGTGTTCCCGTACTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                             <  2:1194630/147‑1 (MQ=14)
                                                                          ggTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtctc                                                                            >  2:3878126/1‑150 (MQ=255)
                                                                          ggTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTtcta                                                                            >  1:5102064/1‑149 (MQ=255)
                                                                           tttttCTCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCTTAACTTATTGTATTTATTTAAAATACCCCCTGAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                           <  2:4768965/149‑1 (MQ=14)
                                                                           gTTTTCGCTATTTATTAAATTTTTCCGGTCTAAGGTGTTTCCGTTCTTCTTCGTCGTAAATTAATGTATTTATTTAAAATACCCCCTTAAACGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                           <  2:270039/150‑1 (MQ=14)
                                                                                ctcTATTTATTAAAATTTTCCGGTTTATGGTTTTTCCGTTCTTCTTCGTCTTAACTTATTTTTTTTATTTAAAATCCCCCATGAAAACAACCGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtttt                                                                      <  2:1597676/148‑1 (MQ=11)
                                                                                 gatATTTTTTAAATTTTTCCGGTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                                     <  2:3989417/148‑1 (MQ=255)
                                                                                   tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                                   >  2:4945399/1‑150 (MQ=255)
                                                                                     tttATTAAAATTTTCGGTTTTATGTTGTTTCCTTTCTATTTCTTCGAAACTTTTTTTTTTTATTAAAAATCCCCCCTCAAAACAACCCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                                                                 <  2:4535492/150‑1 (MQ=9)
                                                                                      ttttGAAAATTTTCGGGTTTCAGGTGTTTCCGTTCTTCGTCGTCGTAACTTACTTTATTTATTTAAAATACCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCg                                                                <  2:2884690/147‑1 (MQ=14)
                                                                                        aTGAAATTTTTCGGGTTTAAGGTGTTTCCGTTCTTTTTCTTCTTAACTTATTGTTTTTATTTAAAATACCCCCTGAAATGAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTg                                                              <  2:2715137/150‑1 (MQ=14)
                                                                                        aTGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTACGTg                                                              >  1:320143/1‑150 (MQ=255)
                                                                                         tGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTACGTgg                                                             >  1:2849451/1‑150 (MQ=255)
                                                                                          gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                            >  1:2565720/1‑150 (MQ=255)
                                                                                             aaTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTACGTGGAAtg                                                         >  2:5044819/1‑150 (MQ=255)
                                                                                               ttttCCGTTTTATGTTGTTCCTTTTCTTTTTGTTCTTAATTTATTTTTTTTTTTTAAAATACCCCATTAAACCAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                                       <  2:1354583/150‑1 (MQ=9)
                                                                                               ttttCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTCATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                                       <  2:2969098/150‑1 (MQ=255)
                                                                                                tttCCTGTTTATGTTTTTTCCTTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTTAAAACAACTGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGc                                                      <  2:744129/150‑1 (MQ=14)
                                                                                                    cGGTTTACGGTTTTTCCGTTATTCTTCGTCTTAACTTAATTTATATATTAAAAATACCCCCTTAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                                  <  2:3772011/150‑1 (MQ=14)
                                                                                                     ggTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTATCTCTGTTTTTGTACGTGGACTTTTCAATgg                                                 >  1:4184728/1‑150 (MQ=14)
                                                                                                      ttttATGGTGTTTCCTTTCTTCTTCGTCTTAACTTATTGTTTTTATTTAAAACCCCCCCTCAAAACAATGTACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                                <  1:4800774/149‑1 (MQ=14)
                                                                                                          ttggTGTTCCCTTTCATCTTCGTCTTAACTTATTTTATTTATTAAAAACCCCCCCTTAAAACAATGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTc                                            <  1:2938796/148‑1 (MQ=14)
                                                                                                          attgtgtTTCCGTTGTTCTTGTTATAAAATATTTTAATTTATAAAAAAAACCCCCATAAAACAAACCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTc                                            <  1:1036075/147‑1 (MQ=9)
                                                                                                          aaGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTCCGTGGTATGTGCAATGGAAGTc                                            >  2:3244158/1‑150 (MQ=255)
                                                                                                              tgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaaca                                        >  2:1324538/1‑150 (MQ=255)
                                                                                                               gtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTTTTTTTGTACGTGGAATGTGCAATGGAAGTcaacaa                                       >  2:4110063/1‑150 (MQ=255)
                                                                                                                     gtacttcttcttcttAACTTATTTTTTTTATTAAAAAACCCCCCTTAAACGAATGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                                 <  1:1165512/147‑1 (MQ=14)
                                                                                                                     gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                                 >  2:1021462/1‑150 (MQ=255)
                                                                                                                         ttatttgtcttAAATTATTTTTTTTATTAAAAACACCCCAACAAAATAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctg                             <  1:1625010/144‑1 (MQ=14)
                                                                                                                             tcgtctAAAATTAATTTATTTATTAAAAACCCCCCCACAAAACAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                         <  1:2706301/150‑1 (MQ=14)
                                                                                                                             tcgtcgTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                         >  2:5320320/1‑150 (MQ=255)
                                                                                                                                 cttaccTTAATGTATTAATTTAAAACACCCCCTGAAAAGAACCGAAACGACAGGTGCTGAAAGCGAGCTTTTGGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt                     <  2:3317674/145‑1 (MQ=14)
                                                                                                                                  gtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTGTCCGTGGAATGTGGAATGGAAGTCAACAAAAAGAAGCAGGCTGACAtt                    >  1:710037/1‑150 (MQ=255)
                                                                                                                                   tAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGTAAGTGCAAATGAAGTCAAAAAAAATAAGCTGGCTGTCAttt                   >  1:4547598/1‑150 (MQ=14)
                                                                                                                                   tAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTACGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCAGGCTGTCAttt                   >  1:279308/1‑150 (MQ=255)
                                                                                                                                    aaCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt                  >  2:2063927/1‑150 (MQ=255)
                                                                                                                                    aaCTTAATGTATTTATTAAAAATACCCCCAGAAAAGAAACGAAACCACAGCTGCTGAAAGCGAGCTTTATGCCCTCTGTCGTTTCCTTTCTCAGATTTTGTCCGTCGAAAGTGCAAAGGAAGTCAACAAAAACCAGCTGGCAGACaatat                  >  2:80618/1‑146 (MQ=14)
                                                                                                                                      cTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg                >  2:251777/1‑150 (MQ=255)
                                                                                                                                      cTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg                >  2:675080/1‑150 (MQ=255)
                                                                                                                                        tattttatttattAAAAAACCCCCAACAAAACAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt              <  1:1977147/145‑1 (MQ=14)
                                                                                                                                         aaTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg             >  2:4428728/1‑150 (MQ=255)
                                                                                                                                          aTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc            >  2:3750312/1‑150 (MQ=255)
                                                                                                                                          aTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc            >  1:3836183/1‑150 (MQ=255)
                                                                                                                                            gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCAGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa          >  2:2302089/1‑150 (MQ=255)
                                                                                                                                                 ttatttaAAATACCCCCTGAAAAGAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTAt      <  2:5130465/150‑1 (MQ=255)
                                                                                                                                                   atttaAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTATGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAATCAGCTGGCTGACATTTTCGGTGCGAGTATcc    >  1:3002630/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGT  >  CP009273/576043‑576340

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: