Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,970,820 G→T 32.9% R193S (CGT→AGT)  motA ← proton conductor component of flagella motor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,970,8200GT32.9% 376.2 / 187.3 231R193S (CGT→AGT) motAproton conductor component of flagella motor
Reads supporting (aligned to +/- strand):  ref base G (94/61);  new base T (50/26);  total (144/87)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 237 total aligned reads displayed.

TTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCG  >  CP009273/1970672‑1970968
                                                                                                                                                    |                                                                                                                                                    
ttGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACg                                                                                                                                                      <  2:2209606/149‑1 (MQ=255)
ttGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACg                                                                                                                                                      <  2:134107/149‑1 (MQ=255)
 tGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGa                                                                                                                                                     <  2:877514/149‑1 (MQ=255)
  gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTAt                                                                                                                                                    <  2:1820283/149‑1 (MQ=255)
  gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGAt                                                                                                                                                    >  1:3347667/1‑149 (MQ=255)
  gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGAt                                                                                                                                                    <  2:3581158/149‑1 (MQ=255)
  gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGAt                                                                                                                                                    <  2:801059/149‑1 (MQ=255)
  gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATACGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGAt                                                                                                                                                    >  1:1149810/1‑149 (MQ=255)
    cGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCg                                                                                                                                                  >  2:3624490/1‑149 (MQ=255)
          gCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGACAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGaa                                                                                                                                            >  2:699066/1‑149 (MQ=255)
          gCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGaa                                                                                                                                            >  1:439460/1‑149 (MQ=255)
           catcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGacc                                                                                                                                           >  2:3458006/1‑147 (MQ=255)
            atcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAAcc                                                                                                                                          <  2:3838168/149‑1 (MQ=255)
            atcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAAcc                                                                                                                                          <  1:3492626/149‑1 (MQ=255)
            atcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAAcc                                                                                                                                          >  2:2957173/1‑149 (MQ=255)
             tcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCt                                                                                                                                         <  2:3347155/149‑1 (MQ=255)
             tcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCt                                                                                                                                         >  2:2733512/1‑149 (MQ=255)
             tcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCt                                                                                                                                         >  2:1556613/1‑149 (MQ=255)
             tcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATACGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCAACGCATGTGCGATAAGCGCACCCAGCACGGCGGCAGGACGATCGGCTGAACCt                                                                                                                                         >  1:1419394/1‑149 (MQ=255)
               atTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTaa                                                                                                                                       <  2:2786253/149‑1 (MQ=255)
                   tGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAAcgcg                                                                                                                                   >  2:1203629/1‑149 (MQ=255)
                   tGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATACGTAAGCCAAAAAAATGCCGAGGAAAGTCCCCAACAAAGCATGTTCGATAATCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAAcccg                                                                                                                                   >  1:84187/1‑149 (MQ=255)
                    gCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGt                                                                                                                                  <  1:2541438/149‑1 (MQ=255)
                     cTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTg                                                                                                                                 >  2:199563/1‑149 (MQ=255)
                     cTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTg                                                                                                                                 >  2:198131/1‑149 (MQ=255)
                        gTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAAc                                                                                                                              >  2:686843/1‑149 (MQ=255)
                        gTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAAc                                                                                                                              <  2:3485002/149‑1 (MQ=255)
                        gTTTCGGCGCTTTTCTGACGTAAAACAGTAGCTCATGTGGAAATAAATCCCTAAGCCAATAAAATGCCGACGAAAGACCCCCCAATCGCATTTGAGATAAGCGCACCCACCTCCGCGGCATGAAGATCGGCTGAACCTAACGCGTTaaa                                                                                                                              >  2:1670719/1‑148 (MQ=255)
                          ttCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGa                                                                                                                            <  1:3530452/149‑1 (MQ=255)
                              gcgcTTTTCTGACGTAAAACAGTCGCTAATGGGTAAATAAATCCGTAAGCAAATAAAATGCCGAGTAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGAcccc                                                                                                                        <  1:1154381/149‑1 (MQ=255)
                               cgcTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGAccccc                                                                                                                       <  1:2678165/149‑1 (MQ=255)
                                gcTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCa                                                                                                                      <  2:1625369/149‑1 (MQ=255)
                                   tttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTa                                                                                                                   <  2:567682/149‑1 (MQ=255)
                                   tttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTa                                                                                                                   >  1:2381858/1‑149 (MQ=255)
                                     tCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACa                                                                                                                 >  1:215618/1‑149 (MQ=255)
                                          cGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGc                                                                                                            <  1:2965585/149‑1 (MQ=255)
                                                aCAGTCGCTAATGGGGCAATAAATACGTAAGCCAAAAAAATGCCGAGGAAAGTCCCCACCATCGCATTTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTAAAGCCGCAACAAt                                                                                                      >  1:1823908/1‑149 (MQ=255)
                                                aCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAAt                                                                                                      >  2:4059413/1‑149 (MQ=255)
                                                 cAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATa                                                                                                     <  1:2458099/149‑1 (MQ=255)
                                                  aGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATAc                                                                                                    <  1:3594532/149‑1 (MQ=255)
                                                   gTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATAcc                                                                                                   >  1:2740616/1‑149 (MQ=255)
                                                     cGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaa                                                                                                 >  1:2640220/1‑149 (MQ=255)
                                                     cGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaa                                                                                                 >  1:573128/1‑149 (MQ=255)
                                                      gCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaaa                                                                                                <  2:4089497/149‑1 (MQ=255)
                                                       cTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAAc                                                                                               >  2:3832639/1‑149 (MQ=255)
                                                           tGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCg                                                                                           >  1:3285699/1‑149 (MQ=255)
                                                              ttaaataaatCCTTAAGCAAATAAAATGCCGAGTAAAGTCCCAACTATGGCTTGTGCGATAAGGGCACCCTGCTGGGCGGCAGGACTATCGGCTTAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGaa                                                                                        <  1:75470/147‑1 (MQ=255)
                                                                    aaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGt                                                                                  <  2:2040049/149‑1 (MQ=255)
                                                                     aatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTc                                                                                 >  2:3067148/1‑149 (MQ=255)
                                                                       tCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGATCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTccc                                                                               >  1:3753452/1‑149 (MQ=255)
                                                                        ccGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATGGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTcccc                                                                              <  1:348390/149‑1 (MQ=255)
                                                                        ccGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTcccc                                                                              <  2:1609888/149‑1 (MQ=255)
                                                                          gTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGa                                                                            >  2:362469/1‑149 (MQ=255)
                                                                             aGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCa                                                                         >  1:1012110/1‑149 (MQ=255)
                                                                             aGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCa                                                                         <  2:227140/149‑1 (MQ=255)
                                                                                 aaTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcgc                                                                     >  2:2219649/1‑149 (MQ=255)
                                                                                 aaTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcgc                                                                     <  2:3478917/149‑1 (MQ=255)
                                                                                  aTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCTCACCCAGCTAGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGcc                                                                    >  1:4071358/1‑149 (MQ=255)
                                                                                  aTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGcc                                                                    >  1:1809768/1‑149 (MQ=255)
                                                                                      aaTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGAc                                                                <  1:3655054/149‑1 (MQ=255)
                                                                                      aaTGCCGAGGAAAGTCCCCACCATCGCATGTGCGAGAAGCGCACCCAGCTCGGCTGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGAc                                                                >  1:3832606/1‑149 (MQ=255)
                                                                                       aggCCGAGTAACGGCCCCACCATCGCTTGTGCTATACGCGCACCCAGCTCGGCGGCAGCACTATCGGCTGACCCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATCCCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                               <  1:1068674/147‑1 (MQ=255)
                                                                                       aTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                               >  1:1536510/1‑149 (MQ=255)
                                                                                            gAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCg                                                          >  1:4088625/1‑149 (MQ=255)
                                                                                             aGGAAAGTCCCCACCATCGCATGTGCGATACGCGCACCCCGCTGGGCGGCAGGACGATCGGCTGACCCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGc                                                         <  1:2760436/149‑1 (MQ=255)
                                                                                              ggAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGcc                                                        >  1:3881904/1‑149 (MQ=255)
                                                                                              ggAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGcc                                                        >  2:1901850/1‑149 (MQ=255)
                                                                                              ggAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGcc                                                        >  2:329167/1‑149 (MQ=255)
                                                                                                aaaGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCgg                                                      <  1:1203674/149‑1 (MQ=255)
                                                                                                   gTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGAc                                                   >  1:3250162/1‑149 (MQ=255)
                                                                                                    tCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACt                                                  >  1:1304151/1‑149 (MQ=255)
                                                                                                    tCCCCACAATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCTTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACt                                                  >  1:3132020/1‑149 (MQ=255)
                                                                                                      cccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTc                                                <  1:3476101/149‑1 (MQ=255)
                                                                                                         accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGc                                             >  1:3436735/1‑149 (MQ=255)
                                                                                                         accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGc                                             >  2:1215063/1‑149 (MQ=255)
                                                                                                         accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGc                                             >  2:2631639/1‑149 (MQ=255)
                                                                                                          ccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGcc                                            >  2:678802/1‑149 (MQ=255)
                                                                                                           caTCGCATGTGCGATAAGCGCACCCCGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCt                                           <  1:874798/149‑1 (MQ=255)
                                                                                                           caTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCt                                           >  1:1223771/1‑149 (MQ=255)
                                                                                                            aTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTc                                          >  2:2459735/1‑149 (MQ=255)
                                                                                                            aTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTc                                          >  1:3416453/1‑149 (MQ=255)
                                                                                                             tCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCg                                         >  1:931201/1‑149 (MQ=255)
                                                                                                                cATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTCACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCtt                                      <  2:487319/149‑1 (MQ=255)
                                                                                                                  tgtgCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGACCGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                    >  2:3503627/1‑149 (MQ=255)
                                                                                                                  tgtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                    <  1:3328220/149‑1 (MQ=255)
                                                                                                                  tgtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                    <  1:809816/149‑1 (MQ=255)
                                                                                                                  tgtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                    >  1:1793821/1‑149 (MQ=255)
                                                                                                                        aTAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGt                              <  2:3840306/149‑1 (MQ=255)
                                                                                                                         tAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtc                             >  1:2074074/1‑149 (MQ=255)
                                                                                                                                  ccAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATctct                    >  2:108164/1‑149 (MQ=255)
                                                                                                                                  ccAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATctct                    <  1:1318071/149‑1 (MQ=255)
                                                                                                                                       tCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTcatc               >  1:1911043/1‑149 (MQ=255)
                                                                                                                                       tCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTcatc               <  2:2416149/149‑1 (MQ=255)
                                                                                                                                         ggcggcAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCACACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTcatcca             <  2:599994/149‑1 (MQ=255)
                                                                                                                                            ggcAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATccacca          <  1:751769/149‑1 (MQ=255)
                                                                                                                                             gcAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCag         >  2:1427921/1‑149 (MQ=255)
                                                                                                                                                  aCGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCtt    <  2:1242572/149‑1 (MQ=255)
                                                                                                                                                    gATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCg  >  1:4008438/1‑149 (MQ=255)
                                                                                                                                                    gATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCg  >  2:1562681/1‑149 (MQ=255)
                                                                                                                                                    gATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCg  >  1:1191330/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCG  >  CP009273/1970672‑1970968

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: