Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,841,044 C→A 87.2% A519A (GCG→GCT hycC ← hydrogenase 3, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,841,0440CA87.2% 764.5 / 110.6 282A519A (GCG→GCThycChydrogenase 3, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base A (115/130);  total (135/147)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 283 total aligned reads displayed.

TTGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTAC  >  CP009273/2840898‑2841193
                                                                                                                                                  |                                                                                                                                                     
  gCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCAt                                                                                                                                                  <  2:3577502/150‑1 (MQ=255)
     tAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAat                                                                                                                                               >  2:336859/1‑150 (MQ=255)
      aGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAatg                                                                                                                                              <  1:2889183/150‑1 (MQ=255)
        tttCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAatgat                                                                                                                                            <  1:4358375/150‑1 (MQ=255)
           cAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGaa                                                                                                                                         <  1:3987167/150‑1 (MQ=255)
           cAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGaa                                                                                                                                         >  1:2296954/1‑150 (MQ=255)
            aGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAAt                                                                                                                                        >  1:2340461/1‑150 (MQ=255)
              cACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGACCGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATgg                                                                                                                                      >  2:4084256/1‑150 (MQ=255)
               aCCGGCGCACACGCCTGTTTCACCGGCTTGGCAAAACCGTGAGGGGTACTCACCATTGATTTTTCGTGGTCGTATCCGCACACCCAGGCCGCTCCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATTgc                                                                                                                                     <  2:1258977/150‑1 (MQ=255)
               aCCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgc                                                                                                                                     >  2:4531540/1‑150 (MQ=255)
                ccGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgca                                                                                                                                    <  1:270886/150‑1 (MQ=255)
                ccGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgca                                                                                                                                    <  1:771952/150‑1 (MQ=255)
                 cGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCCTTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgcag                                                                                                                                   <  2:594115/150‑1 (MQ=255)
                 cGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGgcag                                                                                                                                   >  2:3870586/1‑150 (MQ=255)
                 cGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAAGAGCCATAATGATGAATGgcag                                                                                                                                   <  2:2661472/150‑1 (MQ=255)
                    cgcAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGcagcag                                                                                                                                >  2:931251/1‑150 (MQ=255)
                    cgcAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGcagcag                                                                                                                                >  2:4334567/1‑150 (MQ=255)
                     gcAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCagcagc                                                                                                                               >  1:2876169/1‑150 (MQ=255)
                           gCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGTACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGggcag                                                                                                                         >  1:4522546/1‑150 (MQ=255)
                           gCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGggcag                                                                                                                         >  2:1240051/1‑150 (MQ=255)
                           gCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGggcag                                                                                                                         <  2:2998849/150‑1 (MQ=255)
                            ccTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGggcagg                                                                                                                        <  1:104228/150‑1 (MQ=255)
                            ccTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGggcagg                                                                                                                        <  1:304132/150‑1 (MQ=255)
                              tGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGggcaggca                                                                                                                      >  2:404418/1‑150 (MQ=255)
                                 ttCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATc                                                                                                                   <  2:257461/150‑1 (MQ=255)
                                 ttCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATc                                                                                                                   >  1:1799392/1‑150 (MQ=255)
                                  tCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATca                                                                                                                  <  2:1863146/150‑1 (MQ=255)
                                       ggCATGGCAAAAGCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTATCCGCACACCCAGGCCGCACCGCGGGTACGCGATGGCAATCGATCGCCTTTGCAAAGTGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCagcagc                                                                                                             <  1:1002051/150‑1 (MQ=255)
                                        gCATGGCAAAACCGTGAGCGGTAATCACCAGGGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAgcagca                                                                                                            <  2:4409622/150‑1 (MQ=255)
                                         cATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCaa                                                                                                           <  1:54254/150‑1 (MQ=255)
                                         cATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCaa                                                                                                           <  1:4134093/150‑1 (MQ=255)
                                            ggCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGt                                                                                                        <  2:1338174/150‑1 (MQ=255)
                                             gCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTg                                                                                                       <  1:3910249/150‑1 (MQ=255)
                                              cAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGa                                                                                                      <  1:2292870/150‑1 (MQ=255)
                                               aaaaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAAGGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGat                                                                                                     <  2:147098/150‑1 (MQ=255)
                                                 aaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGatca                                                                                                   <  1:3126022/150‑1 (MQ=255)
                                                 aaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGatca                                                                                                   <  1:1500085/150‑1 (MQ=255)
                                                    cGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACTCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCg                                                                                                >  1:1158024/1‑150 (MQ=255)
                                                    cGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCg                                                                                                <  2:3477577/150‑1 (MQ=255)
                                                     gTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCgg                                                                                               >  2:2409929/1‑150 (MQ=255)
                                                       gAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGtt                                                                                             <  1:3771372/150‑1 (MQ=255)
                                                       gAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGtt                                                                                             <  1:1742398/150‑1 (MQ=255)
                                                          cGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTgag                                                                                          <  1:1727977/150‑1 (MQ=255)
                                                          cGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTgag                                                                                          >  2:1988817/1‑150 (MQ=255)
                                                           ggTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTgaga                                                                                         >  2:1583621/1‑150 (MQ=255)
                                                           ggTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTgaga                                                                                         >  2:2747136/1‑150 (MQ=255)
                                                            gTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaa                                                                                        <  1:2162120/150‑1 (MQ=255)
                                                            gTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaa                                                                                        <  1:4513103/150‑1 (MQ=255)
                                                             tAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaaa                                                                                       >  1:872283/1‑150 (MQ=255)
                                                              aaTCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAAc                                                                                      >  1:3335971/1‑150 (MQ=255)
                                                              aaTCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAAc                                                                                      >  1:3694922/1‑150 (MQ=255)
                                                                 cACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggt                                                                                   >  2:4121496/1‑150 (MQ=255)
                                                                 cACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggt                                                                                   <  2:367662/150‑1 (MQ=255)
                                                                  aCCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtg                                                                                  >  1:3631146/1‑150 (MQ=255)
                                                                   ccATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtgg                                                                                 <  2:3253393/150‑1 (MQ=255)
                                                                    cATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtggt                                                                                <  2:1771018/150‑1 (MQ=255)
                                                                    cATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACATGATCATCGGTTGAGAAACggtggt                                                                                >  2:2156172/1‑150 (MQ=255)
                                                                    cATTGATTTTTCGTGGTCGGAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtggt                                                                                <  2:727007/150‑1 (MQ=255)
                                                                      ttGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGgtgt                                                                              >  1:1858037/1‑150 (MQ=255)
                                                                        gATTTTTCGTGGTCGTAACCGCACCCGCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGACTGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTa                                                                            <  1:4315880/150‑1 (MQ=255)
                                                                        gATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCATACGATCGCCTTTGCAAATAGCCATAATGTTGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTa                                                                            <  2:2227595/150‑1 (MQ=255)
                                                                         aTTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAg                                                                           <  2:170052/150‑1 (MQ=255)
                                                                            tttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCAGCGGTTGAGAAACGGTGGTGTTAGCAg                                                                        >  2:1278374/1‑150 (MQ=255)
                                                                                gTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGTATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTc                                                                    <  1:1129118/150‑1 (MQ=255)
                                                                                 tGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTcc                                                                   >  1:1604529/1‑150 (MQ=255)
                                                                                  ggTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCa                                                                  <  1:1975888/150‑1 (MQ=255)
                                                                                  ggTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCa                                                                  >  2:4198439/1‑150 (MQ=255)
                                                                                     cGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCg                                                               <  1:155242/150‑1 (MQ=255)
                                                                                         aCCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCag                                                           >  1:738363/1‑150 (MQ=255)
                                                                                           cGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCagag                                                         >  1:608178/1‑150 (MQ=255)
                                                                                             cacaCCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGt                                                       <  2:3181313/150‑1 (MQ=255)
                                                                                             cacaCCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGt                                                       >  2:352028/1‑150 (MQ=255)
                                                                                              acacCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTa                                                      <  2:2599067/150‑1 (MQ=255)
                                                                                                acCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAca                                                    <  2:3071061/150‑1 (MQ=255)
                                                                                                    aGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcag                                                >  1:3548684/1‑150 (MQ=255)
                                                                                                    aGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcag                                                <  1:1242233/150‑1 (MQ=255)
                                                                                                        cGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCaga                                            >  2:4498820/1‑150 (MQ=255)
                                                                                                         gCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                           >  2:4155991/1‑150 (MQ=255)
                                                                                                         gCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                           >  2:1063431/1‑150 (MQ=255)
                                                                                                           aCCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagagag                                         <  1:2611041/150‑1 (MQ=255)
                                                                                                            ccGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGc                                        <  1:3501067/150‑1 (MQ=255)
                                                                                                                gggAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCg                                    >  1:49100/1‑150 (MQ=255)
                                                                                                                  gAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGc                                  >  2:2316535/1‑150 (MQ=255)
                                                                                                                  gAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGc                                  >  2:1447866/1‑150 (MQ=255)
                                                                                                                     cgcgATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGt                               <  2:3902472/150‑1 (MQ=255)
                                                                                                                         aTGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGcc                           >  2:1378146/1‑150 (MQ=255)
                                                                                                                             cAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACgg                       >  2:2488804/1‑150 (MQ=255)
                                                                                                                                  gATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcag                  <  1:719201/150‑1 (MQ=255)
                                                                                                                                   aTCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcagc                 <  2:3267495/150‑1 (MQ=255)
                                                                                                                                    tCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcagca                <  2:1799392/150‑1 (MQ=255)
                                                                                                                                       ccTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAaca             <  2:1765501/150‑1 (MQ=255)
                                                                                                                                       ccTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAaca             >  1:2415025/1‑150 (MQ=255)
                                                                                                                                        cTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac            <  2:802641/150‑1 (MQ=255)
                                                                                                                                        cTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac            >  1:3190966/1‑150 (MQ=255)
                                                                                                                                        cTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac            >  1:4423377/1‑150 (MQ=255)
                                                                                                                                        cTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac            >  1:357933/1‑150 (MQ=255)
                                                                                                                                        cTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac            >  2:252393/1‑150 (MQ=255)
                                                                                                                                          ttGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCg          <  1:1470357/150‑1 (MQ=255)
                                                                                                                                             cAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCa       <  1:4462355/150‑1 (MQ=255)
                                                                                                                                              aaaTAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCaa      <  1:336859/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTAC  >  CP009273/2840898‑2841193

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: