Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,287,829 C→A 44.5% S185S (TCG→TCT ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,287,8290CA44.5% 204.4 / inf 272S185S (TCG→TCTccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (74/77);  new base A (43/78);  total (117/155)
Fisher's exact test for biased strand distribution p-value = 2.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 275 total aligned reads displayed.

CGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAA  >  CP009273/2287681‑2287976
                                                                                                                                                    |                                                                                                                                                   
 gtgTTGGTAATCGGCCCATGCGTAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCTGCGTCCACGGACGAGTAAAACGCGCATATGTGCTGTCCAGACGCCCGCTAAGAAAAGAAGCATTGGCAAAAGCAAACGCCACCga                                                                                                                                                   <  2:432173/148‑1 (MQ=255)
 gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCga                                                                                                                                                   <  2:2748841/150‑1 (MQ=255)
  aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACagag                                                                                                                                                  >  2:3542388/1‑150 (MQ=255)
    ttCGTAATAGGCCCATGCGGAACCTAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACTAGTAAAACGCGAATAAGTGCTGTCCAGACGCCCGCTCAGAAAAGAAGCAATGGAAAAAGCAAACGCCACCGAGaa                                                                                                                                                <  2:3282587/150‑1 (MQ=255)
     tCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGaaa                                                                                                                                               >  1:604411/1‑150 (MQ=255)
     tCGTAATAGGCCAAGGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGaaa                                                                                                                                               <  1:1692094/150‑1 (MQ=255)
       gTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAcc                                                                                                                                             >  1:2863280/1‑150 (MQ=255)
       gTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAcc                                                                                                                                             >  1:2766342/1‑150 (MQ=255)
         aaTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCa                                                                                                                                           >  2:3876168/1‑150 (MQ=255)
          aTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCAc                                                                                                                                          >  2:1161389/1‑150 (MQ=255)
           tAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACg                                                                                                                                         <  1:314329/150‑1 (MQ=255)
            tggCCCATGCGGAACCGAGCACGATGCCGAGCTTCTGTAAGATCCACGCGGCCAGCGTCCACGGACGAGTAAAACGCGAATATGTGCTGTCCAGACGCCCGCTCAGAAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGt                                                                                                                                        <  1:1276221/149‑1 (MQ=255)
             ggCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTa                                                                                                                                       <  2:2072591/150‑1 (MQ=255)
              gCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCATTGGCAAAAGCAAACGCCACCGAGAAACCCACGTAc                                                                                                                                      <  2:2482964/150‑1 (MQ=255)
                     cgtaaCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCAAGCGTCAACGGACGATTAAACCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCatat                                                                                                                               <  2:2919879/147‑1 (MQ=255)
                     cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCatat                                                                                                                               <  2:4222759/150‑1 (MQ=255)
                     cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCatat                                                                                                                               >  1:1086901/1‑150 (MQ=255)
                        aaCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAATGCAAACGCCACCGAGAAACCCACGTACCCCATATaaa                                                                                                                            <  2:2152473/150‑1 (MQ=255)
                        aaCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATaaa                                                                                                                            >  2:605376/1‑150 (MQ=255)
                           cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCa                                                                                                                         >  1:1941634/1‑150 (MQ=255)
                             aGCACGTTGCCGAGCGTCTGGAATATCCACGCTGCCAGCGTCCACGGACTTTTAAAAGGCGAATATGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGAAATGGAAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                       <  1:470966/150‑1 (MQ=255)
                              gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCagag                                                                                                                      >  2:2645945/1‑150 (MQ=255)
                                 cGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGAAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCg                                                                                                                   <  1:672376/150‑1 (MQ=255)
                                  gATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCgg                                                                                                                  >  2:3184827/1‑150 (MQ=255)
                                   aTGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGa                                                                                                                 <  2:4110029/150‑1 (MQ=255)
                                   aTGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGa                                                                                                                 <  2:2221268/150‑1 (MQ=255)
                                    tGCCGTGCTTCTGTAAGATCCACGCTGCCAGGGACCACGTATGATTAAAACGCTATTTTGTGCTGTCCAGACGCCCGATCACAAAATACGAATTGGAAAAAGAAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGAt                                                                                                                <  1:14269/150‑1 (MQ=255)
                                    tGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGAt                                                                                                                <  2:1919347/150‑1 (MQ=255)
                                        gAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGaa                                                                                                            >  1:3360148/1‑150 (MQ=255)
                                        gAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGaa                                                                                                            <  1:4270562/150‑1 (MQ=255)
                                        gAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGAAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGaa                                                                                                            <  1:2084445/150‑1 (MQ=255)
                                         aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAg                                                                                                           >  1:3002144/1‑150 (MQ=255)
                                         aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAg                                                                                                           >  2:4082420/1‑150 (MQ=255)
                                           cGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                                                                         <  1:2174578/150‑1 (MQ=255)
                                            gTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCACACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                        <  2:1236089/150‑1 (MQ=255)
                                            gTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                        <  2:2809898/150‑1 (MQ=255)
                                             tCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCa                                                                                                       <  2:848894/150‑1 (MQ=255)
                                             tCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCa                                                                                                       >  2:92011/1‑150 (MQ=255)
                                              cAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAg                                                                                                      <  2:2518044/150‑1 (MQ=255)
                                              cAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACACACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAg                                                                                                      >  2:527309/1‑150 (MQ=255)
                                              cAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAg                                                                                                      <  1:303715/150‑1 (MQ=255)
                                                ggAAGATCCACGCGGCCAGCGTCCACGGACGAGTAACCCGCGCATATGGGCGGTCCAGCCGCCCGCTCAGCAAAGAAGCAAGGGCAAAAGCACACGCCACAGAGAAACCCACGTACCCCATATAAAGCCGAGGCGGATGGAAGATCAGcc                                                                                                    <  1:3340684/150‑1 (MQ=255)
                                                   aGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCg                                                                                                 <  1:1989757/150‑1 (MQ=255)
                                                         aCGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCt                                                                                           <  1:3559243/150‑1 (MQ=255)
                                                            cTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCa                                                                                        <  2:1392907/150‑1 (MQ=255)
                                                             tGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCac                                                                                       >  2:491875/1‑149 (MQ=255)
                                                              gCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAAt                                                                                      <  2:1352409/150‑1 (MQ=255)
                                                                cAGCGTCCACGGACGAGTAACACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAg                                                                                    <  2:3035404/150‑1 (MQ=255)
                                                                cAGCGTCCACGGACGAGTAAAACGCGCATATGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAg                                                                                    <  1:4063516/150‑1 (MQ=255)
                                                                 aGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGc                                                                                   >  2:3720007/1‑150 (MQ=255)
                                                                 aGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGc                                                                                   >  1:2188461/1‑150 (MQ=255)
                                                                  gCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCg                                                                                  <  1:1775174/150‑1 (MQ=255)
                                                                  gCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCg                                                                                  <  1:3964227/150‑1 (MQ=255)
                                                                   cGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCgg                                                                                 >  1:1910586/1‑150 (MQ=255)
                                                                        tcGGTCGAGTAATTCGCGCATAAGTGCGGTCCAGACGCCCGCTCAGCAATGAAGCATTGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATAGAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTaa                                                                            <  1:2374258/149‑1 (MQ=255)
                                                                        aCGGACGAGAAAACCGCGCATACGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCACACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTaa                                                                            <  1:3659519/150‑1 (MQ=255)
                                                                         cGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAg                                                                           >  1:2554762/1‑150 (MQ=255)
                                                                           gACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                                         <  1:3210114/150‑1 (MQ=255)
                                                                                   aaccgcgCATCAGTGCGGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAcca                                                                 <  2:778101/147‑2 (MQ=255)
                                                                                           aTAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCgg                                                         <  2:1848509/150‑1 (MQ=255)
                                                                                             aaGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGa                                                       <  2:3806621/150‑1 (MQ=255)
                                                                                                  cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                  <  2:1425294/150‑1 (MQ=255)
                                                                                                   tGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCg                                                 >  2:1592934/1‑150 (MQ=255)
                                                                                                   tGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCg                                                 >  1:3499713/1‑150 (MQ=255)
                                                                                                   tGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCg                                                 >  2:1978285/1‑150 (MQ=255)
                                                                                                        aGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAAc                                            >  2:749/1‑150 (MQ=255)
                                                                                                          aCGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGt                                          >  2:2664064/1‑150 (MQ=255)
                                                                                                           cGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTg                                         <  1:4121537/150‑1 (MQ=255)
                                                                                                             cccGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCg                                       >  2:3379650/1‑150 (MQ=255)
                                                                                                              ccGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCga                                      <  2:1755177/150‑1 (MQ=255)
                                                                                                               cGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgag                                     >  2:3581576/1‑150 (MQ=255)
                                                                                                               cGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgag                                     >  2:4259221/1‑150 (MQ=255)
                                                                                                                 cTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                                   <  2:101793/150‑1 (MQ=255)
                                                                                                                 cTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                                   <  1:1724203/150‑1 (MQ=255)
                                                                                                                   cAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                                 >  1:63814/1‑150 (MQ=255)
                                                                                                                   cAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                                 >  1:1617746/1‑150 (MQ=255)
                                                                                                                    aGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                                >  1:24430/1‑150 (MQ=255)
                                                                                                                     gCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACg                               <  2:625248/150‑1 (MQ=255)
                                                                                                                      cAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACgg                              >  2:580919/1‑150 (MQ=255)
                                                                                                                      cAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACgg                              <  1:2421079/150‑1 (MQ=255)
                                                                                                                        aagaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt                            >  1:2434068/1‑150 (MQ=255)
                                                                                                                         agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGtt                           <  1:2532828/150‑1 (MQ=255)
                                                                                                                            agCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTaga                        <  2:1076531/150‑1 (MQ=255)
                                                                                                                            agCAATGGCAAAAGCAAACGCCACCGAGAAACACACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTAAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTaga                        >  1:3924454/1‑150 (MQ=255)
                                                                                                                              cAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAgg                      >  1:2671102/1‑150 (MQ=255)
                                                                                                                               aaTGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                     >  1:931786/1‑150 (MQ=255)
                                                                                                                                aTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTa                    <  2:1518313/150‑1 (MQ=255)
                                                                                                                                     aaaaGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAagag               >  1:772514/1‑150 (MQ=255)
                                                                                                                                      aaaGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGTGAGAGAACGGGTTAGAGGTAAAgaga              >  2:1478965/1‑150 (MQ=255)
                                                                                                                                      aaaGCAAACGCCACAGAGAAACCCACGTACCCCATAGAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAgaga              <  2:3456167/150‑1 (MQ=255)
                                                                                                                                      aaaGAAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAgaga              <  1:4218137/150‑1 (MQ=255)
                                                                                                                                       aaGCAAACACCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGaa             >  2:3113256/1‑150 (MQ=255)
                                                                                                                                        aGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTagag                       >  1:3636822/1‑139 (MQ=255)
                                                                                                                                        aGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAAt            <  1:2783379/150‑1 (MQ=255)
                                                                                                                                          cAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGa          <  2:2928189/150‑1 (MQ=255)
                                                                                                                                          cAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGa          <  2:2926492/150‑1 (MQ=255)
                                                                                                                                              cGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAg      >  2:1891976/1‑150 (MQ=255)
                                                                                                                                               gCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGc     >  1:1099581/1‑150 (MQ=255)
                                                                                                                                               gCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGc     >  1:3843149/1‑150 (MQ=255)
                                                                                                                                               gCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGc     >  1:3287014/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAA  >  CP009273/2287681‑2287976

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