Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,612,820 T→A 7.6% S116C (AGT→TGT)  marC ← UPF0056 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,612,8200TA 7.6% 991.2 / 51.1 287S116C (AGT→TGT) marCUPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (131/133);  new base A (12/10);  total (144/143)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

GAAAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCA  >  CP009273/1612671‑1612969
                                                                                                                                                     |                                                                                                                                                     
gAAAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACt                                                                                                                                                       >  1:2414126/1‑150 (MQ=255)
   aaaGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTggg                                                                                                                                                    >  2:4414140/1‑150 (MQ=255)
   aaaGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTggg                                                                                                                                                    >  2:748864/1‑150 (MQ=255)
     aGATCAGCGGCGGCGCAACCATCAGTACCCACTCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGtt                                                                                                                                                  <  1:1343845/150‑1 (MQ=255)
     aGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCACTATTGGCACTGGGtt                                                                                                                                                  >  2:3564161/1‑150 (MQ=255)
      gATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTc                                                                                                                                                 <  2:2729038/150‑1 (MQ=255)
            cggcggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTc                                                                                                                                           >  2:317982/1‑150 (MQ=255)
             ggcggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTcc                                                                                                                                          <  2:2303580/150‑1 (MQ=255)
                 gcgcAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGct                                                                                                                                      <  1:3064730/150‑1 (MQ=255)
                    cAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCtt                                                                                                                                   <  1:3350382/150‑1 (MQ=255)
                     aaCCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTc                                                                                                                                  >  2:3102637/1‑150 (MQ=255)
                      aCCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTcc                                                                                                                                 <  1:1672032/150‑1 (MQ=255)
                        cATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGa                                                                                                                               >  2:945619/1‑150 (MQ=255)
                         aTCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGAc                                                                                                                              >  2:3569718/1‑150 (MQ=255)
                          tCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACt                                                                                                                             >  2:818659/1‑150 (MQ=255)
                            aGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTg                                                                                                                           >  1:1841073/1‑150 (MQ=255)
                                 ccAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCtttt                                                                                                                      <  1:3628579/150‑1 (MQ=255)
                                  cAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTc                                                                                                                     <  2:2952557/150‑1 (MQ=255)
                                   aaTCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCg                                                                                                                    >  2:1537637/1‑150 (MQ=255)
                                   aaTCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCg                                                                                                                    >  2:4024982/1‑150 (MQ=255)
                                    aTCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACAGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGc                                                                                                                   >  2:1018806/1‑150 (MQ=255)
                                      cAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCt                                                                                                                 >  1:717874/1‑150 (MQ=255)
                                        gCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCCTCGCTATGGTTCCCGGACCGGCAGTACGTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTcc                                                                                                               <  1:800504/150‑1 (MQ=255)
                                        gCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTcc                                                                                                               >  2:3512499/1‑150 (MQ=255)
                                          gAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCgg                                                                                                             >  2:3791680/1‑150 (MQ=255)
                                           aaGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                                                                                            >  1:2617156/1‑150 (MQ=255)
                                            aGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTg                                                                                                           <  1:3171415/150‑1 (MQ=255)
                                            aGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTg                                                                                                           >  1:1606009/1‑150 (MQ=255)
                                                 gAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatca                                                                                                      <  1:2759785/150‑1 (MQ=255)
                                                  aaCTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaa                                                                                                     >  1:2792978/1‑150 (MQ=255)
                                                       tGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCt                                                                                                >  1:3860441/1‑150 (MQ=255)
                                                          cGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                                                                                             <  1:1775185/150‑1 (MQ=255)
                                                          cGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                                                                                             <  1:1444831/150‑1 (MQ=255)
                                                            cACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTg                                                                                           >  1:1310557/1‑150 (MQ=255)
                                                                 tgGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCg                                                                                      <  2:4518910/150‑1 (MQ=255)
                                                                 tgGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCg                                                                                      >  2:2391559/1‑150 (MQ=255)
                                                                  ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCgg                                                                                     >  2:3995667/1‑150 (MQ=255)
                                                                   gACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGa                                                                                    >  2:3332174/1‑150 (MQ=255)
                                                                    aCGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaa                                                                                   <  2:1235090/150‑1 (MQ=255)
                                                                     cGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaaa                                                                                  >  1:2152034/1‑150 (MQ=255)
                                                                       ccGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                                                                                >  2:1508531/1‑150 (MQ=255)
                                                                        cGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAagag                                                                               >  2:2479031/1‑150 (MQ=255)
                                                                        cGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAagag                                                                               >  2:2217560/1‑150 (MQ=255)
                                                                          aGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCa                                                                             <  1:2344907/150‑1 (MQ=255)
                                                                           gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCAt                                                                            <  1:2710705/150‑1 (MQ=255)
                                                                            cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc                                                                           <  1:3270855/150‑1 (MQ=255)
                                                                            cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACAGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc                                                                           <  1:2663774/150‑1 (MQ=255)
                                                                              aataatCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCg                                                                         >  1:1207351/1‑150 (MQ=255)
                                                                                    catcGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAAcc                                                                   >  1:85861/1‑150 (MQ=255)
                                                                                       cGCTATGGTTCCCGGACCGGCAGTACCTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAt                                                                <  1:3417679/150‑1 (MQ=255)
                                                                                        gCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATa                                                               >  1:2874289/1‑150 (MQ=255)
                                                                                         cTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACAGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATaa                                                              >  2:2693739/1‑150 (MQ=255)
                                                                                          tATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATaaa                                                             >  1:2540828/1‑150 (MQ=255)
                                                                                            tGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCCCTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGGGAATCAATCGCTTTCTGTTGCGGACAGAGCATCCGAAAACCGATAAAcg                                                           >  2:651455/1‑150 (MQ=255)
                                                                                            tGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAcg                                                           <  2:4144434/150‑1 (MQ=255)
                                                                                            tGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAcg                                                           >  1:1792994/1‑150 (MQ=255)
                                                                                            tGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAcg                                                           <  1:4014093/150‑1 (MQ=255)
                                                                                             ggTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAcgc                                                          <  1:838572/150‑1 (MQ=255)
                                                                                               ttCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGcga                                                        >  1:4502425/1‑150 (MQ=255)
                                                                                                tCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGcgac                                                       >  1:2111728/1‑150 (MQ=255)
                                                                                                tCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACAGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGcgac                                                       <  2:3717619/150‑1 (MQ=255)
                                                                                                 cccGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGcgacg                                                      <  1:3220932/150‑1 (MQ=255)
                                                                                                   cGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGAt                                                    <  2:3077281/150‑1 (MQ=255)
                                                                                                    ggACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACAGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGAtt                                                   <  2:1697224/150‑1 (MQ=255)
                                                                                                           cAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgccgc                                            <  1:3751902/150‑1 (MQ=255)
                                                                                                           cAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgccgc                                            <  1:4094224/150‑1 (MQ=255)
                                                                                                           cAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgccgc                                            >  1:525560/1‑150 (MQ=255)
                                                                                                            aGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgccgcc                                           >  2:57498/1‑150 (MQ=255)
                                                                                                            aGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgccgcc                                           >  1:4164326/1‑150 (MQ=255)
                                                                                                            aGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgccgcc                                           >  2:3219269/1‑150 (MQ=255)
                                                                                                             gTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCt                                          <  2:2713751/150‑1 (MQ=255)
                                                                                                             gTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCt                                          >  2:2879115/1‑150 (MQ=255)
                                                                                                             gTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCt                                          <  2:3968709/150‑1 (MQ=255)
                                                                                                              tACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTg                                         <  2:1458943/150‑1 (MQ=255)
                                                                                                                cTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCa                                       >  2:1832649/1‑150 (MQ=255)
                                                                                                                     cATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCg                                  >  2:2551641/1‑150 (MQ=255)
                                                                                                                       ttGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCa                                >  1:1716162/1‑150 (MQ=255)
                                                                                                                       ttGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCa                                <  1:1766123/150‑1 (MQ=255)
                                                                                                                         gCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGa                              >  2:428081/1‑150 (MQ=255)
                                                                                                                           cAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGAcc                            <  2:1931211/150‑1 (MQ=255)
                                                                                                                               ggCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGa                        >  1:3133343/1‑150 (MQ=255)
                                                                                                                               ggCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGa                        >  1:2568543/1‑150 (MQ=255)
                                                                                                                               ggCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGa                        <  1:945619/150‑1 (MQ=255)
                                                                                                                                 cacaAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAAt                      <  1:3560908/150‑1 (MQ=255)
                                                                                                                                 cacaAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAAt                      <  2:1100407/150‑1 (MQ=255)
                                                                                                                                 cacaAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAAt                      <  1:1953075/150‑1 (MQ=255)
                                                                                                                                 cacaAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAAt                      >  1:2664228/1‑150 (MQ=255)
                                                                                                                                 cacaAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAAt                      >  1:1380718/1‑150 (MQ=255)
                                                                                                                                  acaAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATc                     <  1:1507506/150‑1 (MQ=255)
                                                                                                                                   caAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCg                    <  1:2603034/150‑1 (MQ=255)
                                                                                                                                      aGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGaaa                 >  1:4222942/1‑150 (MQ=255)
                                                                                                                                       ggCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAAt                <  1:3690628/150‑1 (MQ=255)
                                                                                                                                        gCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATg               <  2:19454/150‑1 (MQ=255)
                                                                                                                                         cAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGc              <  2:3762600/150‑1 (MQ=255)
                                                                                                                                         cAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCGCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGc              <  1:412519/150‑1 (MQ=255)
                                                                                                                                           atatTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCa            <  2:2414126/150‑1 (MQ=255)
                                                                                                                                                ggCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGt       >  2:4080255/1‑150 (MQ=255)
                                                                                                                                                ggCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGt       >  2:917656/1‑150 (MQ=255)
                                                                                                                                                  cACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt     <  2:4236113/150‑1 (MQ=255)
                                                                                                                                                     tGGGTCCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCa  >  1:1017241/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCA  >  CP009273/1612671‑1612969

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: