Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 560,498 A→T 9.4% F313I (TTC→ATC)  intD ← DLP12 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273560,4980AT 9.4% 751.9 / 55.1 245F313I (TTC→ATC) intDDLP12 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base A (102/120);  new base T (12/11);  total (114/131)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

Only 100 of 256 total aligned reads displayed.

CCTCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCTTGCCTA  >  CP009273/560349‑560647
                                                                                                                                                     |                                                                                                                                                     
ccTCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGaa                                                                                                                                                       <  1:1302889/150‑1 (MQ=255)
 cTCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGaaa                                                                                                                                                      >  2:19154/1‑150 (MQ=255)
  tCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATAc                                                                                                                                                     <  1:4207516/150‑1 (MQ=255)
  tCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATAc                                                                                                                                                     <  1:844736/150‑1 (MQ=255)
  tCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAAc                                                                                                                                                     <  1:2952363/150‑1 (MQ=255)
  tCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAAc                                                                                                                                                     <  1:1813050/150‑1 (MQ=255)
    gCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACgg                                                                                                                                                   >  2:3018686/1‑150 (MQ=255)
    gCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACgg                                                                                                                                                   >  2:2965475/1‑150 (MQ=255)
        gCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAAt                                                                                                                                               >  2:1103325/1‑150 (MQ=255)
           ctGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCtt                                                                                                                                            >  1:4392148/1‑150 (MQ=255)
            tGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTc                                                                                                                                           >  1:75078/1‑150 (MQ=255)
            tGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTc                                                                                                                                           >  2:2554852/1‑150 (MQ=255)
            tGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTc                                                                                                                                           >  1:2069863/1‑150 (MQ=255)
             gTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCa                                                                                                                                          <  1:2141804/150‑1 (MQ=255)
              tCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATACGGAAATCTTCaa                                                                                                                                         <  1:3435504/150‑1 (MQ=255)
               cAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAAt                                                                                                                                        <  1:4467054/150‑1 (MQ=255)
                aaaTGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAAtt                                                                                                                                       >  1:3848209/1‑150 (MQ=255)
                 aaTGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTc                                                                                                                                      <  2:4540036/150‑1 (MQ=255)
                  aTGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTcc                                                                                                                                     <  1:1682150/150‑1 (MQ=255)
                  aTGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTcc                                                                                                                                     <  2:3185052/150‑1 (MQ=255)
                    gATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATACGGAAATCTTCAATTCCTg                                                                                                                                   <  2:3777030/150‑1 (MQ=255)
                    gATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATACGGAAATCTTCAATTCCTg                                                                                                                                   >  2:293280/1‑150 (MQ=255)
                        aGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCacg                                                                                                                               >  1:1088351/1‑150 (MQ=255)
                         ggCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCacga                                                                                                                              >  2:1422630/1‑150 (MQ=255)
                         ggCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCacga                                                                                                                              <  2:62199/150‑1 (MQ=255)
                           cgcAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCacgacg                                                                                                                            >  2:645901/1‑150 (MQ=255)
                            gcAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCAcgacga                                                                                                                           >  1:493870/1‑150 (MQ=255)
                             cAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACgacgac                                                                                                                          <  2:2418416/150‑1 (MQ=255)
                             cAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACgacgac                                                                                                                          <  1:3227006/150‑1 (MQ=255)
                             cAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACgacgac                                                                                                                          <  2:239867/150‑1 (MQ=255)
                              aaGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACTGAAATCTTCAATTCCTGCACGACGACa                                                                                                                         >  2:4211763/1‑150 (MQ=255)
                              aaGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACa                                                                                                                         >  2:3106775/1‑150 (MQ=255)
                               aGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACaa                                                                                                                        <  2:390477/150‑1 (MQ=255)
                                 gTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGc                                                                                                                      >  2:3723610/1‑150 (MQ=255)
                                  tGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCt                                                                                                                     >  2:3914160/1‑150 (MQ=255)
                                   gAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTg                                                                                                                    <  1:2170845/150‑1 (MQ=255)
                                   gAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTg                                                                                                                    >  1:571441/1‑150 (MQ=255)
                                      cATACCTACGAACCAGTTCTATGGACTCACATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATa                                                                                                                 >  2:4400051/1‑150 (MQ=255)
                                           cTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCat                                                                                                            >  2:1044610/1‑150 (MQ=255)
                                                aaCCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGt                                                                                                       >  1:2516365/1‑150 (MQ=255)
                                                aaCCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGt                                                                                                       >  1:1311521/1‑150 (MQ=255)
                                                   cATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCtt                                                                                                    >  1:2183827/1‑150 (MQ=255)
                                                   cATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCtt                                                                                                    >  1:1307923/1‑150 (MQ=255)
                                                   cATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCtt                                                                                                    <  1:2329967/150‑1 (MQ=255)
                                                   cATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCtt                                                                                                    <  1:236041/150‑1 (MQ=255)
                                                   cATTTCTATGGACTCACATCCGCCCATTTCCTGAAGCACTGATAATGGTACGCCTGACTGAATCAGCCAGCTTGCCAAGGTGTTTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGAACGACGACAAGCAGATAGCCATAATGTcag                                                                                                    >  2:669918/1‑148 (MQ=255)
                                                    aTTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTg                                                                                                   <  2:2164346/150‑1 (MQ=255)
                                                    aTTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCAGGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTg                                                                                                   <  2:1813001/150‑1 (MQ=255)
                                                      ttCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCt                                                                                                 <  2:3641016/150‑1 (MQ=255)
                                                         tATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTc                                                                                              <  1:1009967/150‑1 (MQ=255)
                                                         tATGGACTCCCATCCGCCCATGTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTc                                                                                              <  2:2612011/150‑1 (MQ=255)
                                                          aTGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCg                                                                                             <  2:1618473/150‑1 (MQ=255)
                                                            ggACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGAt                                                                                           <  2:504684/150‑1 (MQ=255)
                                                                  ccATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCAt                                                                                     >  1:3130686/1‑150 (MQ=255)
                                                                     tCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCtt                                                                                  <  2:1177174/150‑1 (MQ=255)
                                                                       cGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTcc                                                                                >  1:1750757/1‑150 (MQ=255)
                                                                       cGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTcc                                                                                >  1:3169704/1‑150 (MQ=255)
                                                                        gCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCt                                                                               <  1:4012458/150‑1 (MQ=255)
                                                                            aTTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGAcc                                                                           >  2:4085005/1‑150 (MQ=255)
                                                                            aTTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCGTCGTGAcc                                                                           >  2:3065441/1‑150 (MQ=255)
                                                                              ttCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGc                                                                         >  2:2940317/1‑150 (MQ=255)
                                                                                  tGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGc                                                                     <  2:1352790/150‑1 (MQ=255)
                                                                                      gCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTg                                                                 <  1:947605/150‑1 (MQ=255)
                                                                                       cACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGa                                                                <  2:1544799/150‑1 (MQ=255)
                                                                                       cACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGa                                                                <  2:1543450/150‑1 (MQ=255)
                                                                                       cACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGa                                                                <  2:516670/150‑1 (MQ=255)
                                                                                        aCTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGAt                                                               >  2:3383716/1‑150 (MQ=255)
                                                                                              aaTGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCa                                                         >  2:1967662/1‑150 (MQ=255)
                                                                                                  ggACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTg                                                     <  1:4511342/150‑1 (MQ=255)
                                                                                                  ggACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTg                                                     <  1:3711058/150‑1 (MQ=255)
                                                                                                    aCGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCt                                                   <  1:4045479/150‑1 (MQ=255)
                                                                                                      gCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCg                                                 <  1:1680755/150‑1 (MQ=255)
                                                                                                         tgactgaATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCtt                                              <  2:3169704/150‑1 (MQ=255)
                                                                                                         tgactgaATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCtt                                              >  1:3606770/1‑150 (MQ=255)
                                                                                                            ctgaATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAgc                                           >  1:3823921/1‑150 (MQ=255)
                                                                                                             tgaATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAgcc                                          >  1:2229126/1‑150 (MQ=255)
                                                                                                             tgaATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGATACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAgcc                                          >  1:2229358/1‑150 (MQ=255)
                                                                                                              gaATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAgccg                                         >  1:1617363/1‑150 (MQ=255)
                                                                                                                   agccagcTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTg                                    <  1:2877504/150‑1 (MQ=255)
                                                                                                                     ccagcTTGCCCAGGTGTGTCTGAGGTCATGGATACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGt                                  <  2:689394/150‑1 (MQ=255)
                                                                                                                     ccagcTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGt                                  <  1:986094/150‑1 (MQ=255)
                                                                                                                      cagcTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGt                                                             >  1:4412151/1‑122 (MQ=255)
                                                                                                                      cagcTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGt                                                             <  2:4412329/122‑1 (MQ=255)
                                                                                                                       agcTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTAt                                <  2:2098545/150‑1 (MQ=255)
                                                                                                                           tGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTAc                            >  1:4326650/1‑150 (MQ=255)
                                                                                                                              ccAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAAGTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTAGCTTAGCCGCCTTGGTATGTACaaa                         <  1:4434626/150‑1 (MQ=255)
                                                                                                                                  gtgtgtCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACAcc                     >  2:2327232/1‑150 (MQ=255)
                                                                                                                                     tgtCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCAt                  >  2:4314737/1‑150 (MQ=255)
                                                                                                                                       tCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCAttt                <  2:2858339/150‑1 (MQ=255)
                                                                                                                                       tCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCAttt                >  1:623770/1‑150 (MQ=255)
                                                                                                                                       tCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCAttt                <  1:2867649/150‑1 (MQ=255)
                                                                                                                                         tGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTtgt              >  2:4498812/1‑150 (MQ=255)
                                                                                                                                            ggTCATGGATACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGAt           <  2:1685302/150‑1 (MQ=255)
                                                                                                                                            ggTCATGGATACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGAt           <  2:2257540/150‑1 (MQ=255)
                                                                                                                                                aTGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCtt       <  2:1261047/150‑1 (MQ=255)
                                                                                                                                                aTGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCtt       <  2:493870/150‑1 (MQ=255)
                                                                                                                                                  ggAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCTTGc     <  2:1307923/150‑1 (MQ=255)
                                                                                                                                                   gAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCTTGcc    <  2:2638600/150‑1 (MQ=255)
                                                                                                                                                    aaaCGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCTTGCCt   >  1:1857334/1‑150 (MQ=255)
                                                                                                                                                     aaCGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCTTGCCTa  >  2:3821869/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTCGCATGCTCTGTCAAATGATTAGGCGCAAGGTGAGCATACCTACGAACCATTTCTATGGACTCCCATCCGCCCATTTCCTGAAGCACTGATAATGGGACGCCTGACTGAATCAGCCAGCTTGCCCAGGTGTGTCTGAGGTCATGGAAACGGAAATCTTCAATTCCTGCACGACGACAAGCTGATAGCCATGATGTCTTGCTGTCGATGCGCATCTTCCTGACCGCAGGCGTTGATGTTCCATCTGCTCGCTTAGCCGCCTTGGTATGTACAAACACCCATTTGTGATGCTTGCCTA  >  CP009273/560349‑560647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: