Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,287,829 C→A 64.6% 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,8290CA64.6% 203.3 / inf 274S185S (TCG→TCTccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (49/48);  new base A (72/105);  total (121/153)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 274 total aligned reads displayed.

CGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAG  >  CP009273/2287681‑2287972
                                                                                                                                                    |                                                                                                                                               
cGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCg                                                                                                                                                >  1:797037/1‑150 (MQ=255)
   gTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATAGCAAAAGCAAACGCCACAgaga                                                                                                                                             <  2:3199145/150‑1 (MQ=255)
    ttCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGaa                                                                                                                                            <  2:1083058/150‑1 (MQ=255)
      cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATCAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAc                                                                                                                                          <  2:1736294/150‑1 (MQ=255)
      cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAc                                                                                                                                          <  2:3987411/150‑1 (MQ=255)
      cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAAc                                                                                                                                          <  1:1966538/150‑1 (MQ=255)
         aaTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCa                                                                                                                                       >  1:3506535/1‑150 (MQ=255)
         aaTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCa                                                                                                                                       >  1:1235338/1‑150 (MQ=255)
           tAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACg                                                                                                                                     >  2:4029168/1‑150 (MQ=255)
           tAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACg                                                                                                                                     <  2:3066772/150‑1 (MQ=255)
             ggCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTa                                                                                                                                   <  1:3548768/150‑1 (MQ=255)
              gCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAc                                                                                                                                  <  2:1275296/150‑1 (MQ=255)
              gCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAc                                                                                                                                  <  2:23631/150‑1 (MQ=255)
                 cATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAcccc                                                                                                                               >  2:3509896/1‑150 (MQ=255)
                  aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGAGCTGTCCAGACGCCCGCTCAGCAAAGAAGCTATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCa                                                                                                                              <  2:4110575/150‑1 (MQ=255)
                   tGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCat                                                                                                                             >  1:3811315/1‑150 (MQ=255)
                     cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCatat                                                                                                                           <  1:2818983/150‑1 (MQ=255)
                      ggAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGACACCCACGTACCCCAtata                                                                                                                          >  2:2431206/1‑150 (MQ=255)
                      ggAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCAtata                                                                                                                          >  2:2430957/1‑150 (MQ=255)
                      ggAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCAtata                                                                                                                          >  1:1098222/1‑150 (MQ=255)
                         aCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAg                                                                                                                       <  1:2449401/150‑1 (MQ=255)
                            gAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCag                                                                                                                    >  1:1288740/1‑150 (MQ=255)
                             aGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   >  1:3625726/1‑150 (MQ=255)
                              gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCagag                                                                                                                  <  2:56885/150‑1 (MQ=255)
                              gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCagag                                                                                                                  <  2:3487020/150‑1 (MQ=255)
                               cACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAgg                                                                                                                 <  2:2412248/150‑1 (MQ=255)
                               cACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAgg                                                                                                                 <  2:2509213/150‑1 (MQ=255)
                                aCGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGc                                                                                                                <  1:4259318/150‑1 (MQ=255)
                                 cGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCg                                                                                                               >  1:3246120/1‑150 (MQ=255)
                                  gATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCgg                                                                                                              <  1:3165914/150‑1 (MQ=255)
                                  gATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCgg                                                                                                              >  2:2722221/1‑150 (MQ=255)
                                     gCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATg                                                                                                           <  2:2203578/150‑1 (MQ=255)
                                      ccGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGCTgg                                                                                                          >  2:4352863/1‑150 (MQ=255)
                                         aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAg                                                                                                       >  2:58108/1‑150 (MQ=255)
                                           cGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                                                                     <  2:3005891/150‑1 (MQ=255)
                                           cGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                                                                     <  1:2528363/150‑1 (MQ=255)
                                            gTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                    >  2:2159897/1‑150 (MQ=255)
                                            gTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                    >  1:4409766/1‑150 (MQ=255)
                                            gTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                    <  2:1708904/150‑1 (MQ=255)
                                                     aTCCACGCTGCCAGCGTCCACGGACGAGTATAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGa                                                                                           <  2:1400989/150‑1 (MQ=255)
                                                       ccACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATc                                                                                         <  1:1348726/150‑1 (MQ=255)
                                                          cGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTg                                                                                      <  1:3485115/150‑1 (MQ=255)
                                                           gCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGc                                                                                     <  1:701828/150‑1 (MQ=255)
                                                           gCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGc                                                                                     <  1:701727/150‑1 (MQ=255)
                                                             tGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCaa                                                                                   >  1:4170331/1‑150 (MQ=255)
                                                             tGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCaa                                                                                   >  1:2072072/1‑150 (MQ=255)
                                                               ccAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATa                                                                                 <  2:1619204/150‑1 (MQ=255)
                                                                cAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAg                                                                                <  1:4415191/150‑1 (MQ=255)
                                                                cAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAg                                                                                <  1:171307/150‑1 (MQ=255)
                                                                  gCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCg                                                                              <  2:1698781/150‑1 (MQ=255)
                                                                   cGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCgg                                                                             <  2:29810/150‑1 (MQ=255)
                                                                   cGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCgg                                                                             >  2:2453402/1‑150 (MQ=255)
                                                                   cGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATAAAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCAGCAATAGCgg                                                                             >  1:3543010/1‑150 (MQ=255)
                                                                       cACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTa                                                                         >  2:3681275/1‑150 (MQ=255)
                                                                               aGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcgcg                                                                 >  1:3046386/1‑150 (MQ=255)
                                                                                 tAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                               <  1:622789/150‑1 (MQ=255)
                                                                                 tAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                               >  2:2250846/1‑150 (MQ=255)
                                                                                  aaaaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAcc                                                              >  2:4520846/1‑150 (MQ=255)
                                                                                    aaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                            >  1:1141646/1‑150 (MQ=255)
                                                                                     aCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTc                                                           <  1:149707/150‑1 (MQ=255)
                                                                                        cgcATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAAt                                                        <  2:2998596/150‑1 (MQ=255)
                                                                                             aaGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGa                                                   >  2:3585447/1‑150 (MQ=255)
                                                                                             aaGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGa                                                   <  1:3763993/150‑1 (MQ=255)
                                                                                                 ggagTCCATACGCCCGCTCAGCAAATAAGCAATGGCAAAAGCAAACGCCACAGAGAAAACCAAGTACCCCATATAAATCAGATGCGGATGGTAGATCATCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGtt                                               <  2:1284699/147‑1 (MQ=255)
                                                                                                     tCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGACGTTCGGc                                           >  1:570461/1‑150 (MQ=255)
                                                                                                      ccAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCa                                          <  1:4008722/150‑1 (MQ=255)
                                                                                                      ccAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAACCCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCa                                          <  2:481932/150‑1 (MQ=255)
                                                                                                       cAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                         <  1:1586858/150‑1 (MQ=255)
                                                                                                        aGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAAc                                        >  2:4458885/1‑150 (MQ=255)
                                                                                                        aGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAAc                                        <  1:3088041/150‑1 (MQ=255)
                                                                                                        aGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAAc                                        >  1:480367/1‑150 (MQ=255)
                                                                                                         gACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACg                                       >  1:1653356/1‑150 (MQ=255)
                                                                                                           cGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTg                                     >  1:2626418/1‑150 (MQ=255)
                                                                                                            gCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCACCGTGc                                    >  2:4264404/1‑150 (MQ=255)
                                                                                                                   cAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                             >  2:738881/1‑150 (MQ=255)
                                                                                                                    aGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                            <  2:471672/150‑1 (MQ=255)
                                                                                                                       aaaGAAGCAATGGCAAAAGCAAACGCTACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACggg                         <  1:934767/150‑1 (MQ=255)
                                                                                                                       aaaGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACggg                         >  1:4319688/1‑150 (MQ=255)
                                                                                                                        aagaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt                        <  2:137179/150‑1 (MQ=255)
                                                                                                                         agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAATCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGtt                       >  1:3693348/1‑150 (MQ=255)
                                                                                                                           aagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTag                     >  2:2858798/1‑150 (MQ=255)
                                                                                                                            agCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTaga                    >  1:2952546/1‑150 (MQ=255)
                                                                                                                            agCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTaga                    <  2:1615757/150‑1 (MQ=255)
                                                                                                                             gCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTagag                   <  1:1270726/150‑1 (MQ=255)
                                                                                                                              cAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAgg                  <  1:453153/150‑1 (MQ=255)
                                                                                                                               aaTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                 <  1:3254626/150‑1 (MQ=255)
                                                                                                                               aaTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                 <  1:2798521/150‑1 (MQ=255)
                                                                                                                                aTGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTa                <  2:901304/150‑1 (MQ=255)
                                                                                                                                 tGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaa               <  1:2181078/150‑1 (MQ=255)
                                                                                                                                   gCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAag             <  1:1410398/150‑1 (MQ=255)
                                                                                                                                    cAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAaga            >  2:4080976/1‑150 (MQ=255)
                                                                                                                                    cAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAaga            >  2:2214202/1‑150 (MQ=255)
                                                                                                                                       aaGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGaa         <  2:188464/150‑1 (MQ=255)
                                                                                                                                         gCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATg       <  1:3145399/150‑1 (MQ=255)
                                                                                                                                         gCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATg       <  1:3145419/150‑1 (MQ=255)
                                                                                                                                          cAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGa      >  1:3355169/1‑150 (MQ=255)
                                                                                                                                            aaCGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAAc    <  1:2202033/150‑1 (MQ=255)
                                                                                                                                            aaCGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAAc    <  2:3057594/150‑1 (MQ=255)
                                                                                                                                              cGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAg  <  1:3841803/150‑1 (MQ=255)
                                                                                                                                              cGCCACAGAGAAACCCACGTAACCCATATAACCCAGAGGAGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGTGTTAAGATAGCGACCTTATATCGGGAAGTTCGGCAACGTGCGAGAGATCGGGTTAGAGGTTAAGAGAATGAACAg  <  2:4024530/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAG  >  CP009273/2287681‑2287972

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: