Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,287,829 C→A 52.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,8290CA52.6% 43.1 / inf 224S185S (TCG→TCTccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (41/65);  new base A (65/53);  total (106/118)
Fisher's exact test for biased strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 230 total aligned reads displayed.

CGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAA  >  CP009273/2287681‑2287978
                                                                                                                                                    |                                                                                                                                                     
 gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                                     >  1:3728570/1‑150 (MQ=255)
 gAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                                     >  2:2595986/1‑150 (MQ=255)
  aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACagag                                                                                                                                                    >  2:2541416/1‑150 (MQ=255)
   gTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAgaga                                                                                                                                                   >  1:3056460/1‑150 (MQ=255)
      cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAc                                                                                                                                                >  2:1536817/1‑150 (MQ=255)
       gTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAAcc                                                                                                                                               >  2:1612405/1‑150 (MQ=255)
         aaTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCa                                                                                                                                             <  1:3941819/150‑1 (MQ=255)
          aTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCAc                                                                                                                                            <  2:3490290/150‑1 (MQ=255)
            aGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGt                                                                                                                                          >  2:1869438/1‑150 (MQ=255)
                 cATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAACCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAcccc                                                                                                                                     <  2:1466272/150‑1 (MQ=255)
                  aTGCGTAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACTGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCa                                                                                                                                    <  1:3928712/150‑1 (MQ=255)
                  aTGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCa                                                                                                                                    <  2:922612/150‑1 (MQ=255)
                   tGCGGAACCGAGCACGATGCCGAGGGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACGCACGTACCCCat                                                                                                                                   <  1:3211967/150‑1 (MQ=255)
                    gCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCata                                                                                                                                  <  1:200893/150‑1 (MQ=255)
                     cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCatat                                                                                                                                 >  1:727244/1‑150 (MQ=255)
                        aaCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATaaa                                                                                                                              >  2:2675189/1‑150 (MQ=255)
                           cGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCa                                                                                                                           >  1:651951/1‑150 (MQ=255)
                            gAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCag                                                                                                                          >  2:592861/1‑150 (MQ=255)
                             aGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCACACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                         >  2:916430/1‑150 (MQ=255)
                             aGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                         <  2:1594308/150‑1 (MQ=255)
                                 cGATGCCGAGCTTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCAAGACGCCCGCTCAGAAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCg                                                                                                                     <  1:1775651/150‑1 (MQ=255)
                                 cGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCg                                                                                                                     >  1:1522409/1‑150 (MQ=255)
                                     gCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATg                                                                                                                 >  2:886860/1‑150 (MQ=255)
                                      ccGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATgg                                                                                                                >  2:3241030/1‑150 (MQ=255)
                                       cGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGa                                                                                                               >  2:1794888/1‑150 (MQ=255)
                                        gAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGaa                                                                                                              >  1:2529115/1‑150 (MQ=255)
                                           cGTCAGGAAGATCCACGCTGCCAGCTTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCTATGTCAAAAGCAAACGCCACCGAGAAACCCCCGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                                                                           <  2:3760247/150‑1 (MQ=255)
                                           cGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                                                                           >  2:2957142/1‑150 (MQ=255)
                                             tCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAGCCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCa                                                                                                         <  2:3637028/150‑1 (MQ=255)
                                               aGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGc                                                                                                       <  1:832142/150‑1 (MQ=255)
                                               aGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGc                                                                                                       <  2:771865/150‑1 (MQ=255)
                                                 gAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGccc                                                                                                     >  2:1968984/1‑150 (MQ=255)
                                                  aaGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGcccc                                                                                                    <  2:460332/150‑1 (MQ=255)
                                                      tCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGAt                                                                                                <  2:4117918/150‑1 (MQ=255)
                                                      tCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGAt                                                                                                >  2:176370/1‑150 (MQ=255)
                                                          cGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTg                                                                                            <  1:4051020/150‑1 (MQ=255)
                                                              gCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAAt                                                                                        <  1:1020535/150‑1 (MQ=255)
                                                                cAGCGTCCACTGACGAGTAATACGCGCATAAGTGCTGTCCAGACGCGCGCTGAGCTTATATGCAATGGAAAAAGCAATAGCGAAAGAGAATCACACGTAACCCATATAATGCTTATTCGTATGGAAGATGAGCCCCGGATCCTGCAAttt                                                                                      <  2:2197849/150‑3 (MQ=255)
                                                                cAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAg                                                                                      <  2:515253/150‑1 (MQ=255)
                                                                 aGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGc                                                                                     >  2:1251101/1‑150 (MQ=255)
                                                                       cACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTa                                                                               >  1:3170493/1‑150 (MQ=255)
                                                                          ggACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGa                                                                            >  2:694654/1‑150 (MQ=255)
                                                                           gACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                                           <  2:243172/150‑1 (MQ=255)
                                                                            aCGAGTAAAACGCGCATAAGTGCTGTCCAGCCGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATc                                                                          <  2:3395043/150‑1 (MQ=255)
                                                                            aCGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATc                                                                          <  2:1153956/150‑1 (MQ=255)
                                                                            aCGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATc                                                                          <  2:789532/150‑1 (MQ=255)
                                                                               aGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcgcg                                                                       >  1:2803336/1‑150 (MQ=255)
                                                                                 tAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                                     <  2:1001926/150‑1 (MQ=255)
                                                                                  aaaaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAcc                                                                    >  1:701171/1‑150 (MQ=255)
                                                                                    aaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                                  <  2:3593407/150‑1 (MQ=255)
                                                                                    aaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                                  >  1:2674172/1‑150 (MQ=255)
                                                                                    aaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                                  <  1:578453/150‑1 (MQ=255)
                                                                                    aaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                                  >  2:1708661/1‑150 (MQ=255)
                                                                                     aCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTc                                                                 <  1:2410575/150‑1 (MQ=255)
                                                                                      cgcgCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCa                                                                <  1:2218486/150‑1 (MQ=255)
                                                                                        cgcATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAAt                                                              <  1:2074833/150‑1 (MQ=255)
                                                                                          cATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCg                                                            >  2:2065704/1‑150 (MQ=255)
                                                                                           aTAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCgg                                                           >  1:4036242/1‑150 (MQ=255)
                                                                                            tAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCggg                                                          >  2:1557228/1‑150 (MQ=255)
                                                                                             aaGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGa                                                         >  2:468887/1‑150 (MQ=255)
                                                                                              aGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGaa                                                        <  1:3940503/150‑1 (MQ=255)
                                                                                              aGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGaa                                                        <  1:1705104/150‑1 (MQ=255)
                                                                                                 gCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGtt                                                     <  1:912752/150‑1 (MQ=255)
                                                                                                  cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                    <  1:1610563/150‑1 (MQ=255)
                                                                                                  cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTc                                                    >  1:623416/1‑150 (MQ=255)
                                                                                                   tGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCg                                                   <  2:1073535/150‑1 (MQ=255)
                                                                                                    gTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCgg                                                  >  1:3284772/1‑150 (MQ=255)
                                                                                                    gTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCgg                                                  >  1:3365985/1‑150 (MQ=255)
                                                                                                     tCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGc                                                 <  2:3518613/150‑1 (MQ=255)
                                                                                                      ccAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCa                                                <  1:1156409/150‑1 (MQ=255)
                                                                                                       cAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                               >  1:1859799/1‑150 (MQ=255)
                                                                                                        aGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAAc                                              <  1:3729488/150‑1 (MQ=255)
                                                                                                           cGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTg                                           >  2:2335273/1‑150 (MQ=255)
                                                                                                            gCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGc                                          <  1:509231/150‑1 (MQ=255)
                                                                                                              ccGCTCATCAAATAAGCAATGGCAAAATCAAACGCCACCGAGAAACCCACGTACCCAATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCTACCTTCAATCGGGAAGTTCGGCAACGTGCga                                        <  2:1304750/150‑1 (MQ=255)
                                                                                                                    aGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                                  <  2:2062522/150‑1 (MQ=255)
                                                                                                                     gCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACg                                 <  1:2938811/150‑1 (MQ=255)
                                                                                                                     gCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACg                                 >  1:4063611/1‑150 (MQ=255)
                                                                                                                       aaaGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACggg                               <  1:2385862/150‑1 (MQ=255)
                                                                                                                        aagaagCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt                              >  1:59170/1‑150 (MQ=255)
                                                                                                                         agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGtt                             <  1:2509645/150‑1 (MQ=255)
                                                                                                                         agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGtt                             <  2:3045571/150‑1 (MQ=255)
                                                                                                                          gaagCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAAGAGCGGGTTCAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTa                            >  2:431893/1‑150 (MQ=255)
                                                                                                                               aaTGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                       <  2:2898851/150‑1 (MQ=255)
                                                                                                                               aaTGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt                       <  1:1406014/150‑1 (MQ=255)
                                                                                                                                 tGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaa                     <  2:1609155/150‑1 (MQ=255)
                                                                                                                                 tGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaa                     >  2:793534/1‑150 (MQ=255)
                                                                                                                                 tGGCAAAAGCAAAAGCCAACGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCATGCAATAGCGGGTTCAGATCGCGACCTTAAATCGGGAAGTTAGGCAACGTGCGAGAGAACGGGTTAGAGGTaa                     >  2:367267/1‑150 (MQ=255)
                                                                                                                                   gCAAAAGCAAACGCCACCGAGAAACCCACGAACCCCAGATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTACAATCGGGAAGAACGGCAACGTGCGAGAGAACGGGATCGAGGTAAag                   >  2:3075872/1‑150 (MQ=255)
                                                                                                                                    cAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGACCGGGTTAGAGGTAAaga                  >  2:3807424/1‑150 (MQ=255)
                                                                                                                                     aaaaGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAagag                 >  2:2922560/1‑150 (MQ=255)
                                                                                                                                        aGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAAt              <  2:4108130/150‑1 (MQ=255)
                                                                                                                                        aGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAAt              <  2:3295040/150‑1 (MQ=255)
                                                                                                                                        aGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAAt              >  1:1684858/1‑150 (MQ=255)
                                                                                                                                                ccACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCa      <  1:3347414/150‑1 (MQ=255)
                                                                                                                                                ccACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCa      >  1:3443129/1‑150 (MQ=255)
                                                                                                                                                ccACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCa      <  1:1912027/150‑1 (MQ=255)
                                                                                                                                                 cACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCaa     <  1:671134/150‑1 (MQ=255)
                                                                                                                                                    cGAGAAACCCACGTACCCCATATCAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCaaaaa  <  2:4025056/150‑1 (MQ=255)
                                                                                                                                                    cGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCaaaaa  <  1:3823224/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAA  >  CP009273/2287681‑2287978

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