Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,070,058 C→A 17.4% A90A (GCC→GCA yeeS → CP4‑44 prophage; putative DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,070,0580CA17.4% 675.4 / 105.3 240A90A (GCC→GCAyeeSCP4‑44 prophage; putative DNA repair protein
Reads supporting (aligned to +/- strand):  ref base C (125/72);  new base A (25/17);  total (151/89)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCAGGTTTT  >  CP009273/2069913‑2070207
                                                                                                                                                 |                                                                                                                                                     
cGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTg                                                                                                                                                   <  2:1283480/150‑1 (MQ=255)
 gggACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAAtgctgctgc                                                                                                                                                  >  2:4886295/1‑150 (MQ=255)
 gggACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGc                                                                                                                                                  >  2:2589654/1‑150 (MQ=255)
      tGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTgg                                                                                                                                             >  2:951611/1‑150 (MQ=255)
        gAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTa                                                                                                                                           >  2:3430713/1‑150 (MQ=255)
          gCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACt                                                                                                                                         <  1:698534/150‑1 (MQ=255)
          gCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACt                                                                                                                                         <  1:3323823/150‑1 (MQ=255)
           cGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTg                                                                                                                                        >  2:4517394/1‑150 (MQ=255)
            gTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTgg                                                                                                                                       >  2:2866188/1‑150 (MQ=255)
              gaagaaGTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGca                                                                                                                                     <  2:5192254/150‑1 (MQ=255)
                  aaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGcacaca                                                                                                                                 >  2:3560166/1‑150 (MQ=255)
                  aaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGcacaca                                                                                                                                 >  1:1924863/1‑150 (MQ=255)
                  aaTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGcacaca                                                                                                                                 >  1:4494060/1‑150 (MQ=255)
                   cttCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACaa                                                                                                                                <  2:961036/149‑1 (MQ=255)
                   aTTCCGGTTTCGGTATCAGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTTATTAAACTCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGAACACaa                                                                                                                                <  2:3583169/150‑1 (MQ=255)
                   aTTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACaa                                                                                                                                >  1:3435374/1‑150 (MQ=255)
                     tCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGAGTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATc                                                                                                                              <  1:3593525/150‑1 (MQ=255)
                      ccGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCa                                                                                                                             <  1:2709870/150‑1 (MQ=255)
                         ggTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCAccc                                                                                                                          <  1:1955116/150‑1 (MQ=255)
                           tGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGt                                                                                                                        >  2:1202462/1‑150 (MQ=255)
                            gCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTc                                                                                                                       >  1:4165661/1‑150 (MQ=255)
                            gCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTc                                                                                                                       >  1:3096972/1‑150 (MQ=255)
                             cTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTcc                                                                                                                      <  1:2437280/150‑1 (MQ=255)
                                tATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACCATCACCCGTCCGGt                                                                                                                   >  2:2166632/1‑150 (MQ=255)
                                tATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGAGTAAACGCGCCCTGTCCCACAATGCAGCTGCCGTGGTACTGGCCCACAATCACCCGTCCGGt                                                                                                                   >  2:3910184/1‑150 (MQ=255)
                                 aTCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCGACACCGGCACCATCAACCGCACGGACGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTg                                                                                                                  <  2:4397619/150‑1 (MQ=255)
                                 aTCTGAACAACCAGAATCAGCTGAGTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGCACTGGCACACAATCACCCGTCCGGTg                                                                                                                  <  2:4220039/150‑1 (MQ=255)
                                    tGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAg                                                                                                               >  1:19035/1‑150 (MQ=255)
                                    tGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAg                                                                                                               >  1:19458/1‑150 (MQ=255)
                                        caacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTcac                                                                                                           >  2:4241585/1‑150 (MQ=255)
                                        caacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTcac                                                                                                           >  1:3344703/1‑150 (MQ=255)
                                         aacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTcaca                                                                                                          >  1:3442524/1‑150 (MQ=255)
                                         aacCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTcaca                                                                                                          >  2:767921/1‑150 (MQ=255)
                                             aGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCCTCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCa                                                                                                      <  2:680579/150‑1 (MQ=255)
                                             aGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCa                                                                                                      >  2:5038431/1‑150 (MQ=255)
                                             aGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCAAACCCa                                                                                                      >  2:5038233/1‑150 (MQ=255)
                                                aTCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTa                                                                                                   <  1:2701694/150‑1 (MQ=255)
                                                 tCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTaa                                                                                                  <  2:3975039/150‑1 (MQ=255)
                                                 tCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTaa                                                                                                  <  1:2866188/150‑1 (MQ=255)
                                                 tCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTaa                                                                                                  <  1:768907/150‑1 (MQ=255)
                                                  cAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGACGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTCCTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAg                                                                                                 >  2:2803641/1‑150 (MQ=255)
                                                        aTTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGAc                                                                                           >  1:2137900/1‑150 (MQ=255)
                                                          tGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCg                                                                                         <  2:1740655/150‑1 (MQ=255)
                                                           gCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCgg                                                                                        <  1:1192295/150‑1 (MQ=255)
                                                            ccGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGc                                                                                       >  1:4980676/1‑150 (MQ=255)
                                                            ccGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGc                                                                                       >  1:2088953/1‑150 (MQ=255)
                                                              ggTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGAAGCTGCCGTGGTACTGTCACACAATCAACCGTCCGGTGAAGTCACACCCAGTAATGCAGACCGGCtt                                                                                     >  2:4928453/1‑150 (MQ=255)
                                                                tGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTAt                                                                                   <  2:1413072/150‑1 (MQ=255)
                                                                  aaaCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCa                                                                                 <  2:3791313/150‑1 (MQ=255)
                                                                  aaaCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCa                                                                                 >  2:427123/1‑150 (MQ=255)
                                                                    aCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCAcc                                                                               >  2:1699680/1‑150 (MQ=255)
                                                                    aCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCAcc                                                                               >  2:1699693/1‑150 (MQ=255)
                                                                    aCCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCAcc                                                                               <  2:4041770/150‑1 (MQ=255)
                                                                      ccTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGa                                                                             >  1:4762992/1‑150 (MQ=255)
                                                                        tctTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAAc                                                                           >  1:1713688/1‑150 (MQ=255)
                                                                        tctTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAAc                                                                           <  1:3915845/150‑1 (MQ=255)
                                                                          ttCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGt                                                                         >  2:3319894/1‑150 (MQ=255)
                                                                           tCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTc                                                                        >  2:1976312/1‑150 (MQ=255)
                                                                             aCCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTg                                                                      >  2:3454328/1‑150 (MQ=255)
                                                                               cGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGt                                                                    >  2:1060347/1‑150 (MQ=255)
                                                                                ggCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTa                                                                   >  1:4436101/1‑150 (MQ=255)
                                                                                  cACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACa                                                                 >  1:4685756/1‑150 (MQ=255)
                                                                                   aCCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAg                                                                >  1:3218923/1‑150 (MQ=255)
                                                                                      aTCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCa                                                             >  2:4102271/1‑150 (MQ=255)
                                                                                      aTCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCa                                                             <  1:3113689/150‑1 (MQ=255)
                                                                                        cAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACt                                                           <  1:3719343/150‑1 (MQ=255)
                                                                                           ccGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTggg                                                        <  2:2934300/150‑1 (MQ=255)
                                                                                             gCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGcc                                                      >  2:3311445/1‑150 (MQ=255)
                                                                                                 ggAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggt                                                  <  1:2083696/150‑1 (MQ=255)
                                                                                                 ggAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggt                                                  >  1:2906537/1‑150 (MQ=255)
                                                                                                 ggAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggt                                                  >  2:3397539/1‑150 (MQ=255)
                                                                                                 ggAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggt                                                  <  2:2380756/150‑1 (MQ=255)
                                                                                                  gAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggtg                                                 >  2:1494905/1‑150 (MQ=255)
                                                                                                   aaGTCCATCACCGGGAAGTGAATAAACGCGCCCTGTACCACAATGCCGCTGACTTTGTACTGGCACACAATCACCCGTCCTGTGAAGTCACACCCAGAAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCtggagt                                                >  2:3474907/1‑147 (MQ=255)
                                                                                                    aGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGa                                               >  2:1942065/1‑150 (MQ=255)
                                                                                                     gTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGat                                              <  2:1806473/150‑1 (MQ=255)
                                                                                                        cATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATc                                           <  1:3840693/150‑1 (MQ=255)
                                                                                                           ccccGGGAAGTGATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCgg                                        <  2:645600/150‑1 (MQ=255)
                                                                                                               gggAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGc                                    >  1:1062930/1‑150 (MQ=255)
                                                                                                               gggAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGc                                    >  1:1062152/1‑150 (MQ=255)
                                                                                                                      gATTAAACGCGCCCTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCa                             >  2:3015930/1‑150 (MQ=255)
                                                                                                                        ttAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATc                           >  1:4820716/1‑150 (MQ=255)
                                                                                                                         tAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCt                          >  2:2448871/1‑150 (MQ=255)
                                                                                                                         tAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACCGGGCCTGGTGGATATCCGGGTGCCGGACCATCt                          <  1:2535023/150‑1 (MQ=255)
                                                                                                                           aaCGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGa                        <  1:1574280/150‑1 (MQ=255)
                                                                                                                           aaCGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGa                        <  2:4165661/150‑1 (MQ=255)
                                                                                                                           aaCGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGa                        >  1:3220839/1‑150 (MQ=255)
                                                                                                                             cgcgCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATa                      <  2:364590/150‑1 (MQ=255)
                                                                                                                             cgcgCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATa                      <  1:2292545/150‑1 (MQ=255)
                                                                                                                              gcgcCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCTCCGAACGGCTGGTCCAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGTTAg                     >  2:3134250/1‑150 (MQ=255)
                                                                                                                              gcgcCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAg                     <  2:4451203/150‑1 (MQ=255)
                                                                                                                                 cccTGTACCACAATGCAGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCg                  >  2:2250543/1‑150 (MQ=255)
                                                                                                                                  ccTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCgg                 <  2:3664508/150‑1 (MQ=255)
                                                                                                                                   cTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGt                <  1:4102271/150‑1 (MQ=255)
                                                                                                                                    tGTACCACACTGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTg               <  1:1660480/150‑1 (MQ=255)
                                                                                                                                    tGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTg               <  1:4268783/150‑1 (MQ=255)
                                                                                                                                         cacaATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAAc          >  2:485542/1‑150 (MQ=255)
                                                                                                                                           caATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCa        >  1:1616646/1‑150 (MQ=255)
                                                                                                                                             aTGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCAgg      <  1:1463484/150‑1 (MQ=255)
                                                                                                                                                 cGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCAGGtttt  <  1:2002563/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCATCCCCGGGAAGTGATTAAACGCGCCCTGTACCACAATGCCGCTGCCGTGGTACTGGCACACAATCACCCGTCCGGTGAAGTCACACCCAGTAAGGCAGACCGGCTTATCACCGAACGTCTGGTACAGGCACTGGGCCTGGTGGATATCCGGGTGCCGGACCATCTGATAGTCGGTGGCAACCAGGTTTT  >  CP009273/2069913‑2070207

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: