Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,558,906 G→T 9.7% A362E (GCG→GAG)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,558,9060GT 9.7% 898.5 / 52.1 248A362E (GCG→GAG) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base G (113/111);  new base T (11/13);  total (124/124)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 268 total aligned reads displayed.

AACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAA  >  CP009273/1558758‑1559053
                                                                                                                                                    |                                                                                                                                                   
aaCACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgc                                                                                                                                                    >  1:3741068/1‑150 (MQ=255)
 acacCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                                                                                                                   >  2:2642042/1‑150 (MQ=255)
  cacCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcca                                                                                                                                                  >  2:1883677/1‑150 (MQ=255)
  cacCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcca                                                                                                                                                  >  2:4159947/1‑150 (MQ=255)
  cacCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcca                                                                                                                                                  <  1:2598725/150‑1 (MQ=255)
   acCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagctccag                                                                                                                                                 >  2:1967902/1‑150 (MQ=255)
   acCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagctccag                                                                                                                                                 >  1:16288/1‑150 (MQ=255)
     cGATGAGATACACCACGGGAGAGTCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagcg                                                                                                                                               <  1:3160694/150‑1 (MQ=255)
      gATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCcagcgccagcgc                                                                                                                                              >  2:2328699/1‑150 (MQ=255)
       aTGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGcgcg                                                                                                                                             <  2:3563356/150‑1 (MQ=255)
         gagaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGc                                                                                                                                           >  1:80726/1‑150 (MQ=255)
          agaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGcc                                                                                                                                          >  2:1842078/1‑150 (MQ=255)
          agaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGcc                                                                                                                                          <  1:927954/150‑1 (MQ=255)
            aTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCat                                                                                                                                        >  2:1035460/1‑150 (MQ=255)
            aTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCat                                                                                                                                        >  1:2271996/1‑150 (MQ=255)
                  cacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctg                                                                                                                                  <  1:2252564/150‑1 (MQ=255)
                   acGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctgg                                                                                                                                 <  1:3336637/150‑1 (MQ=255)
                    cGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggc                                                                                                                                >  2:2730302/1‑150 (MQ=255)
                      ggAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggctg                                                                                                                              <  1:1711868/150‑1 (MQ=255)
                               gCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                                                                                                                     >  2:78701/1‑150 (MQ=255)
                               gCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                                                                                                                     >  2:2367449/1‑150 (MQ=255)
                                 tttGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAc                                                                                                                   <  1:2560055/150‑1 (MQ=255)
                                   tGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc                                                                                                                 >  1:2317569/1‑150 (MQ=255)
                                   tGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc                                                                                                                 >  2:3010714/1‑150 (MQ=255)
                                   tGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc                                                                                                                 >  1:1383179/1‑150 (MQ=255)
                                      cGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGtt                                                                                                              >  2:375284/1‑150 (MQ=255)
                                         ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGa                                                                                                           >  2:2576940/1‑150 (MQ=255)
                                         ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGa                                                                                                           <  1:4169367/150‑1 (MQ=255)
                                         ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGa                                                                                                           <  1:1752824/150‑1 (MQ=255)
                                          cAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGAt                                                                                                          >  2:1437606/1‑150 (MQ=255)
                                          cAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGAt                                                                                                          <  2:1549799/150‑1 (MQ=255)
                                            ggTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaa                                                                                                        <  1:972367/150‑1 (MQ=255)
                                            ggTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaa                                                                                                        <  2:4540042/150‑1 (MQ=255)
                                                   gAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCt                                                                                                 >  1:1302534/1‑150 (MQ=255)
                                                   gAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCt                                                                                                 >  1:927089/1‑150 (MQ=255)
                                                    aGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTg                                                                                                <  2:4049112/150‑1 (MQ=255)
                                                     ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTgg                                                                                               <  1:2474236/150‑1 (MQ=255)
                                                      gTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGt                                                                                              <  1:280633/150‑1 (MQ=255)
                                                      gTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGt                                                                                              <  2:2039515/150‑1 (MQ=255)
                                                        aaTTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGttt                                                                                            >  1:738644/1‑150 (MQ=255)
                                                         aTTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTc                                                                                           >  1:3598487/1‑150 (MQ=255)
                                                          ttGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCt                                                                                          <  2:3709564/150‑1 (MQ=255)
                                                             tgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCt                                                                                       <  2:3122076/150‑1 (MQ=255)
                                                              gCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTc                                                                                      >  1:3906068/1‑150 (MQ=255)
                                                               cAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTcc                                                                                     <  1:4247728/150‑1 (MQ=255)
                                                               cAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTcc                                                                                     <  2:2925698/150‑1 (MQ=255)
                                                                aGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCt                                                                                    <  2:2118998/150‑1 (MQ=255)
                                                                  gTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGa                                                                                  >  1:3998615/1‑150 (MQ=255)
                                                                   ttATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGAc                                                                                 <  1:853781/150‑1 (MQ=255)
                                                                       tGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAgg                                                                             >  2:3427127/1‑150 (MQ=255)
                                                                           aTTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGtttt                                                                         >  2:63673/1‑150 (MQ=255)
                                                                           aTTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGtttt                                                                         >  1:1861327/1‑150 (MQ=255)
                                                                           aTTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGtttt                                                                         >  1:2677114/1‑150 (MQ=255)
                                                                           aTTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGtttt                                                                         <  1:4154764/150‑1 (MQ=255)
                                                                            tttGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTa                                                                        <  1:1115398/150‑1 (MQ=255)
                                                                               ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTCTCTGCTCCTGACGAGGTTTTAAtt                                                                     <  1:2447234/150‑1 (MQ=255)
                                                                                    aCCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAg                                                                <  2:2890120/150‑1 (MQ=255)
                                                                                     ccGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAgg                                                               <  2:1828343/150‑1 (MQ=255)
                                                                                      cGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGa                                                              >  2:657498/1‑150 (MQ=255)
                                                                                       ggTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGAt                                                             >  1:3819153/1‑150 (MQ=255)
                                                                                       ggTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGAt                                                             >  1:4454493/1‑150 (MQ=255)
                                                                                       ggTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGAt                                                             >  1:3622935/1‑150 (MQ=255)
                                                                                         tCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATcc                                                           <  1:1192339/150‑1 (MQ=255)
                                                                                          cATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCTTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATccc                                                          >  1:2727142/1‑150 (MQ=255)
                                                                                          cATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATccc                                                          <  2:1652918/150‑1 (MQ=255)
                                                                                          cATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATccc                                                          >  2:195707/1‑150 (MQ=255)
                                                                                              ggATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCa                                                      >  1:2756592/1‑150 (MQ=255)
                                                                                                aTCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAgg                                                    >  1:4144695/1‑150 (MQ=255)
                                                                                                 tCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGgc                                                   >  1:4489880/1‑150 (MQ=255)
                                                                                                   aaaTTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGgcgc                                                 <  1:1374380/150‑1 (MQ=255)
                                                                                                       tGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATACAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCa                                             <  2:4461355/150‑1 (MQ=255)
                                                                                                         gATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCCTc                                           >  1:3834874/1‑150 (MQ=255)
                                                                                                         gATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATc                                           >  2:422275/1‑150 (MQ=255)
                                                                                                          aTGAGTTGTTCAATATGCTGACGGGTTTTTTCGTGTTCCAGCTCCAGCGCGGCCAGAGGCGGGCTGATCTCTGCCACGCGTTCGATACAGGCGCTGGTTTCTGCTCCTGACGAGGGTTTTAGTTGCAGGATCCCCGCAGGCGCGCCATCa                                          >  2:3177823/1‑150 (MQ=255)
                                                                                                              gttgttCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCt                                      <  1:1985978/150‑1 (MQ=255)
                                                                                                                tgttCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGa                                    >  2:520776/1‑150 (MQ=255)
                                                                                                                 gttCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGAc                                   <  1:328817/150‑1 (MQ=255)
                                                                                                                 gttCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGAc                                   <  1:4300698/150‑1 (MQ=255)
                                                                                                                        atGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGACTGGt                            >  1:3068355/1‑150 (MQ=255)
                                                                                                                        atGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGt                            >  1:442214/1‑150 (MQ=255)
                                                                                                                             gACGGCTTTTTTCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTg                       <  2:226452/150‑1 (MQ=255)
                                                                                                                             gACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTg                       >  2:4307455/1‑150 (MQ=255)
                                                                                                                              aCGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGa                      >  2:1144022/1‑150 (MQ=255)
                                                                                                                               cGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGAc                     >  1:352377/1‑150 (MQ=255)
                                                                                                                               cGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGAc                     >  2:3314086/1‑150 (MQ=255)
                                                                                                                                ggCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGAcc                    >  1:428396/1‑150 (MQ=255)
                                                                                                                                 gCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCa                   <  2:3998615/150‑1 (MQ=255)
                                                                                                                                 gCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCa                   <  1:4093343/150‑1 (MQ=255)
                                                                                                                                    tttttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCt                >  1:1824552/1‑150 (MQ=255)
                                                                                                                                     ttgtCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCtt               <  2:3842112/147‑1 (MQ=255)
                                                                                                                                      tttCCTGTTCCAGCTCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCttt              <  2:2969164/150‑1 (MQ=255)
                                                                                                                                      tttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCttt              >  2:4263503/1‑150 (MQ=255)
                                                                                                                                      tttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCttt              <  1:269734/150‑1 (MQ=255)
                                                                                                                                      tttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCttt              <  1:3237261/150‑1 (MQ=255)
                                                                                                                                       ttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTg             <  1:2372395/150‑1 (MQ=255)
                                                                                                                                       ttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTg             >  2:2399049/1‑150 (MQ=255)
                                                                                                                                       ttCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTg             <  1:108555/150‑1 (MQ=255)
                                                                                                                                              tccagctccagcGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCAttt      >  2:2517242/1‑150 (MQ=255)
                                                                                                                                               ccagcgccagcgCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCAtttt     <  2:714640/150‑1 (MQ=255)
                                                                                                                                                  gcgccagcgcGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGaa  <  2:638399/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAA  >  CP009273/1558758‑1559053

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: