Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,459,603 A→T 100% intergenic (+24/+5) gspO → / ← bfr bifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,459,6030AT95.1% 953.3 / 23.9 299intergenic (+24/+5)gspO/bfrbifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (4/11);  new base T (141/143);  total (145/154)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 300 total aligned reads displayed.

CATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCAT  >  CP009273/3459456‑3459750
                                                                                                                                                   |                                                                                                                                                   
cATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACg                                                                                                                                                   <  2:4403149/150‑1 (MQ=255)
 aTTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCgg                                                                                                                                                  >  1:431415/1‑150 (MQ=255)
       tACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAAc                                                                                                                                            <  2:834682/150‑1 (MQ=255)
       tACGCCATCGTATCACAACTACTGCAACGCCGGACTAAAACCACTAACGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTAATAACCTGCAGGGTACTCCATCCCGACTCGGCTCAAc                                                                                                                                            >  2:3968016/1‑150 (MQ=255)
         cGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACct                                                                                                                                          >  2:3701067/1‑150 (MQ=255)
          gCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACctt                                                                                                                                         >  2:4476085/1‑150 (MQ=255)
          gCCATCGTTTCACAACATCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACctt                                                                                                                                         >  2:4748730/1‑150 (MQ=255)
           ccATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttc                                                                                                                                        >  1:1639422/1‑150 (MQ=255)
              tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCttcttc                                                                                                                                     <  2:5044826/150‑1 (MQ=255)
              tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCttcttc                                                                                                                                     <  2:1966939/150‑1 (MQ=255)
                gTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTcgc                                                                                                                                   >  1:4806594/1‑150 (MQ=255)
                 tttCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTcgcg                                                                                                                                  <  1:4096142/150‑1 (MQ=255)
                   tCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGa                                                                                                                                >  2:5130388/1‑150 (MQ=255)
                      caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCt                                                                                                                             <  2:3501971/150‑1 (MQ=255)
                      caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATGAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCt                                                                                                                             <  1:3069837/150‑1 (MQ=255)
                       aaCTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCTACTCGGCTCAACCTTCTTCGCGGATCtg                                                                                                                            <  2:4561832/150‑1 (MQ=255)
                       aaCTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCtg                                                                                                                            >  1:4020164/1‑150 (MQ=255)
                           tCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCt                                                                                                                        <  2:727443/150‑1 (MQ=255)
                             tGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                      >  2:4014077/1‑150 (MQ=255)
                             tGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                      <  1:4387088/150‑1 (MQ=255)
                              gCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGc                                                                                                                     <  2:1548996/150‑1 (MQ=255)
                                aaCGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                   >  1:1165574/1‑150 (MQ=255)
                                   gCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATa                                                                                                                >  2:2583441/1‑150 (MQ=255)
                                       gACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAttt                                                                                                            <  2:4667691/150‑1 (MQ=255)
                                         cTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTg                                                                                                          >  1:2493508/1‑150 (MQ=255)
                                             aaCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAgg                                                                                                      <  2:2767808/150‑1 (MQ=255)
                                              aCCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGc                                                                                                     >  1:793155/1‑150 (MQ=255)
                                                 aCTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCa                                                                                                  >  2:2922470/1‑150 (MQ=255)
                                                 aCTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCa                                                                                                  <  1:4182866/150‑1 (MQ=255)
                                                  cTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAt                                                                                                 >  1:3251791/1‑150 (MQ=255)
                                                     tCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtctt                                                                                              <  2:171787/150‑1 (MQ=255)
                                                       gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtcttct                                                                                            <  1:1626258/150‑1 (MQ=255)
                                                       gCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCAtcttct                                                                                            <  1:4059410/150‑1 (MQ=255)
                                                        cATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTg                                                                                           >  1:3390262/1‑150 (MQ=255)
                                                           tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAAt                                                                                        <  1:2007017/150‑1 (MQ=255)
                                                            cGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATc                                                                                       <  2:1169223/150‑1 (MQ=255)
                                                            cGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATc                                                                                       <  2:1169189/150‑1 (MQ=255)
                                                                 cGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATc                                                                                  >  1:539078/1‑150 (MQ=255)
                                                                  gTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCa                                                                                 <  1:1120363/150‑1 (MQ=255)
                                                                  gTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCa                                                                                 >  1:171887/1‑150 (MQ=255)
                                                                      cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGtt                                                                             >  2:765319/1‑150 (MQ=255)
                                                                      cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGtt                                                                             >  2:2329496/1‑150 (MQ=255)
                                                                      cTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGtt                                                                             <  1:2272347/150‑1 (MQ=255)
                                                                       ttGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTc                                                                            >  1:2895181/1‑150 (MQ=255)
                                                                          cgcTCGTCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGt                                                                         <  1:4449217/150‑1 (MQ=255)
                                                                          cgcTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGt                                                                         <  2:3922837/150‑1 (MQ=255)
                                                                          cgcTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGt                                                                         >  1:3933483/1‑150 (MQ=255)
                                                                           gcTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGtt                                                                        <  2:1467264/150‑1 (MQ=255)
                                                                            cTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGttt                                                                       <  2:4621923/150‑1 (MQ=255)
                                                                             tCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTc                                                                      <  1:488337/150‑1 (MQ=255)
                                                                             tCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTc                                                                      >  1:3654454/1‑150 (MQ=255)
                                                                               ggCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcca                                                                    <  2:4138188/150‑1 (MQ=255)
                                                                               ggCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcca                                                                    <  1:624953/150‑1 (MQ=255)
                                                                                 cagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagc                                                                  >  1:2589446/1‑150 (MQ=255)
                                                                                  agcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                                 >  1:3981646/1‑150 (MQ=255)
                                                                                      tGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTc                                                             <  1:3721642/150‑1 (MQ=255)
                                                                                       gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCg                                                            <  1:3427739/150‑1 (MQ=255)
                                                                                       gATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCg                                                            <  2:30442/150‑1 (MQ=255)
                                                                                        aTAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGa                                                           <  2:431348/150‑1 (MQ=255)
                                                                                         tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                          <  1:3455710/150‑1 (MQ=255)
                                                                                         tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                          >  2:2128025/1‑150 (MQ=255)
                                                                                             ccTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                      <  1:2355728/150‑1 (MQ=255)
                                                                                              cTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGc                                                     >  1:475797/1‑150 (MQ=255)
                                                                                               tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                    <  1:2210364/150‑1 (MQ=255)
                                                                                               tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGcc                                                    >  1:2980738/1‑150 (MQ=255)
                                                                                                    tATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCcttctt                                               >  1:4529760/1‑150 (MQ=255)
                                                                                                     aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCttcttc                                              <  1:1011089/150‑1 (MQ=255)
                                                                                                       cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcat                                            >  2:1879795/1‑150 (MQ=255)
                                                                                                       cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcat                                            >  1:2403917/1‑150 (MQ=255)
                                                                                                       cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcat                                            >  2:4065317/1‑150 (MQ=255)
                                                                                                          gCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCAc                                         >  1:556793/1‑150 (MQ=255)
                                                                                                           ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCAca                                        <  2:3252052/150‑2 (MQ=255)
                                                                                                              ggATCTCTTATTACCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                     <  1:873683/150‑1 (MQ=255)
                                                                                                              ggATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                     >  2:1026231/1‑150 (MQ=255)
                                                                                                              ggATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                     <  1:778405/150‑1 (MQ=255)
                                                                                                                       attaACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTat                            <  2:1518999/150‑1 (MQ=255)
                                                                                                                       attaACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTat                            >  1:2058069/1‑150 (MQ=255)
                                                                                                                         taACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                          <  1:3590515/150‑1 (MQ=255)
                                                                                                                          aaCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                         <  1:128041/150‑1 (MQ=255)
                                                                                                                          aaCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                         <  2:2116284/150‑1 (MQ=255)
                                                                                                                          aaCCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                         <  2:1124481/150‑1 (MQ=255)
                                                                                                                             cGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                      <  2:2167995/150‑1 (MQ=255)
                                                                                                                             cGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCAGTTTCCAGCCCGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                      >  2:3623517/1‑150 (MQ=255)
                                                                                                                               ggAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCatat                    >  1:4539930/1‑150 (MQ=255)
                                                                                                                               ggAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCatat                    <  1:2918608/150‑1 (MQ=255)
                                                                                                                               ggAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCatat                    <  1:2144314/150‑1 (MQ=255)
                                                                                                                                 aGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcg                  <  1:2753928/150‑1 (MQ=255)
                                                                                                                                   ggTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgcg                <  2:1229469/150‑1 (MQ=255)
                                                                                                                                   ggTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgcg                >  2:485460/1‑150 (MQ=255)
                                                                                                                                      tctcCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCt             <  2:1476018/150‑1 (MQ=255)
                                                                                                                                       ctccctccCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg            <  1:972292/150‑1 (MQ=255)
                                                                                                                                        tccctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATGTCGCGGCTGa           >  1:811702/1‑150 (MQ=255)
                                                                                                                                          cctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg         >  2:4362340/1‑150 (MQ=255)
                                                                                                                                          cctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg         >  1:350075/1‑150 (MQ=255)
                                                                                                                                            tcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTa       >  2:2032805/1‑150 (MQ=255)
                                                                                                                                             cccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCGTCATCACGCACAATTTCTCTCATCATATCGCGGCTGACGTaa      >  2:2553334/1‑150 (MQ=255)
                                                                                                                                               cGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc    <  1:4975154/150‑1 (MQ=255)
                                                                                                                                                gACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCa   <  1:3202546/150‑1 (MQ=255)
                                                                                                                                                gACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCa   >  1:4280735/1‑150 (MQ=255)
                                                                                                                                                 acacGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGATCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCAt  >  2:3914528/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCAT  >  CP009273/3459456‑3459750

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: