Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,549,855 C→A 9.9% A71A (GCC→GCA groL → chaperonin GroEL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,549,8550CA 9.9% 591.6 / 42.1 191A71A (GCC→GCAgroLchaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (83/89);  new base A (8/11);  total (91/100)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 193 total aligned reads displayed.

GTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGA  >  NC_002947/1549706‑1550004
                                                                                                                                                     |                                                                                                                                                     
gTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGcc                                                                                                                                                       <  2:4099859/150‑1 (MQ=255)
        tGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGt                                                                                                                                               >  1:478822/1‑150 (MQ=255)
           cgcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCaa                                                                                                                                            <  2:4540669/150‑1 (MQ=255)
            gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                                                           <  2:986892/150‑1 (MQ=255)
            gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                                                           >  1:2265336/1‑150 (MQ=255)
             cgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAgg                                                                                                                                          >  2:2084577/1‑150 (MQ=255)
               gTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGaa                                                                                                                                        <  2:157984/150‑1 (MQ=255)
                tAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAg                                                                                                                                       >  2:2328306/1‑150 (MQ=255)
                tAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAg                                                                                                                                       >  1:826195/1‑150 (MQ=255)
                tAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAg                                                                                                                                       <  1:2676993/150‑1 (MQ=255)
                  aaaGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGtt                                                                                                                                     <  2:1128911/150‑1 (MQ=255)
                   aaGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCCAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTg                                                                                                                                    >  2:2694434/1‑150 (MQ=255)
                    aGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCCTCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGc                                                                                                                                   >  1:3239919/1‑150 (MQ=255)
                     gCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCt                                                                                                                                  >  1:3780929/1‑150 (MQ=255)
                        aCCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTcc                                                                                                                               <  2:3447224/150‑1 (MQ=255)
                          ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                                                             >  2:1841530/1‑150 (MQ=255)
                          ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                                                             <  1:5116843/150‑1 (MQ=255)
                          ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                                                             <  1:102829/150‑1 (MQ=255)
                           cTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAg                                                                                                                            >  1:3004371/1‑150 (MQ=255)
                            tGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAgg                                                                                                                           >  1:4043571/1‑150 (MQ=255)
                             gggCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGc                                                                                                                          <  2:1666617/150‑1 (MQ=255)
                               gCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCa                                                                                                                        >  2:825511/1‑150 (MQ=255)
                               gCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCa                                                                                                                        <  2:3296486/150‑1 (MQ=255)
                                 ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                                                                      <  2:3131783/150‑1 (MQ=255)
                                    aaaGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgac                                                                                                                   <  2:4265654/150‑1 (MQ=255)
                                           gTAACGTGGTACTGGACAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCg                                                                                                            <  2:652320/150‑1 (MQ=255)
                                             aaCGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGt                                                                                                          <  1:4671214/150‑1 (MQ=255)
                                               cGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGa                                                                                                        <  1:3081179/150‑1 (MQ=255)
                                                gTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGAc                                                                                                       >  2:3553132/1‑150 (MQ=255)
                                                gTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGAc                                                                                                       <  1:2373756/150‑1 (MQ=255)
                                                     aCTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcac                                                                                                  <  2:2783539/150‑1 (MQ=255)
                                                     aCTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcac                                                                                                  <  1:2007647/150‑1 (MQ=255)
                                                      cTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcacc                                                                                                 <  2:2931514/150‑1 (MQ=255)
                                                       tGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcacca                                                                                                >  2:571734/1‑150 (MQ=255)
                                                       tGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCAAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTACCAAGGCCAACGACGCTGCCGGTGACGGcacca                                                                                                >  1:3354161/1‑150 (MQ=255)
                                                           caagtccTTCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACTGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGcaccaccac                                                                                            >  2:464935/7‑150 (MQ=255)
                                                           cAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGTCGCTGCCGGTGACGGcaccaccac                                                                                            >  2:3169502/1‑150 (MQ=255)
                                                           cAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcacaaccac                                                                                            >  1:4223407/1‑150 (MQ=255)
                                                            aagtccTTCGGCGCTCCACGCATCACCCAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACACGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGcaccaccact                                                                                           <  2:2655268/145‑2 (MQ=255)
                                                            aaGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCaccaccacc                                                                                           >  2:4874990/1‑150 (MQ=255)
                                                              gagCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGc                                                                                         <  2:2439735/150‑1 (MQ=255)
                                                                gCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTa                                                                                       <  1:73660/150‑1 (MQ=255)
                                                                 cTTCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGAc                                                                                      >  2:3330156/1‑150 (MQ=255)
                                                                 cTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAc                                                                                      >  1:4455866/1‑150 (MQ=255)
                                                                     ggCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGtt                                                                                  <  1:5017356/150‑1 (MQ=255)
                                                                       cgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                                                                <  2:5063796/150‑1 (MQ=255)
                                                                         cgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                                                              <  1:946084/150‑1 (MQ=255)
                                                                             gACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCa                                                                          <  2:181912/150‑1 (MQ=255)
                                                                             acgcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCa                                                                          <  1:1565016/147‑1 (MQ=255)
                                                                                  tCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGcta                                                                     <  2:4427236/150‑3 (MQ=255)
                                                                                   cACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTAt                                                                    <  2:4975839/150‑1 (MQ=255)
                                                                                   cACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCAt                                                                    >  2:3666998/1‑150 (MQ=255)
                                                                                     ccAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTCCCGTTCTGGCTCAGGCCATCg                                                                  >  2:4501491/1‑150 (MQ=255)
                                                                                     ccAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCg                                                                  <  1:432185/150‑1 (MQ=255)
                                                                                      cAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                                                                 <  1:444202/150‑1 (MQ=255)
                                                                                      cAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                                                                 >  1:2617048/1‑150 (MQ=255)
                                                                                         ggACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttca                                                              <  2:701142/150‑4 (MQ=255)
                                                                                           aCGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttcagg                                                            <  1:1412907/150‑6 (MQ=255)
                                                                                             ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                                                          <  1:3036372/150‑1 (MQ=255)
                                                                                             ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                                                          <  2:2430891/150‑1 (MQ=255)
                                                                                              gCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                         <  2:2636610/150‑1 (MQ=255)
                                                                                              gCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                         <  2:2920658/150‑1 (MQ=255)
                                                                                                gTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                       <  1:2503456/150‑1 (MQ=255)
                                                                                                gTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                       <  1:4143943/150‑1 (MQ=255)
                                                                                                 tttCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGcc                                                      <  2:4059209/150‑1 (MQ=255)
                                                                                                  ttCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt                                                     <  2:3587681/150‑1 (MQ=255)
                                                                                                    ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                   <  2:3366777/150‑1 (MQ=255)
                                                                                                    ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                   <  1:768683/150‑1 (MQ=255)
                                                                                                         gCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGcc                                              >  2:271570/1‑150 (MQ=255)
                                                                                                          ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCg                                             >  2:1142793/1‑150 (MQ=255)
                                                                                                          ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTAGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCg                                             >  1:1478632/1‑150 (MQ=255)
                                                                                                           cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                            >  2:1456333/1‑150 (MQ=255)
                                                                                                           cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                            >  2:209639/1‑150 (MQ=255)
                                                                                                           cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                            >  2:1748621/1‑150 (MQ=255)
                                                                                                           cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                            >  1:2304575/1‑150 (MQ=255)
                                                                                                            aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                           >  1:4873434/1‑150 (MQ=255)
                                                                                                            aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                           <  1:11089/150‑1 (MQ=255)
                                                                                                            aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                           >  2:2236052/1‑150 (MQ=255)
                                                                                                                         tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                              >  2:509266/1‑150 (MQ=255)
                                                                                                                         tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                              >  2:1698551/1‑150 (MQ=255)
                                                                                                                         tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGATCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                              >  2:2357838/1‑150 (MQ=255)
                                                                                                                             aGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAAccc                          >  1:1012024/1‑150 (MQ=255)
                                                                                                                                 gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                      <  2:2304508/150‑1 (MQ=255)
                                                                                                                                   tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGa                    >  2:730900/1‑150 (MQ=255)
                                                                                                                                   taacGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGACCGTGACGGAACAACCACCGCTACCGTTCTGGCTCAGGACATCGTCAACGAAGTCCTGAAAGCCGACGCTGCCGGCATGAACCCGAAGGa                    >  2:1219206/4‑150 (MQ=255)
                                                                                                                                    ttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAc                   >  1:3464578/1‑150 (MQ=255)
                                                                                                                                    ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGAc                   >  1:4834328/1‑150 (MQ=255)
                                                                                                                                     tCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAcc                  >  2:3803403/1‑150 (MQ=255)
                                                                                                                                      cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCt                 <  2:2498129/150‑1 (MQ=255)
                                                                                                                                       gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                >  2:470752/1‑150 (MQ=255)
                                                                                                                                       gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                >  2:1010671/1‑150 (MQ=255)
                                                                                                                                       gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                >  2:358350/1‑150 (MQ=255)
                                                                                                                                         aaaCATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGACCTGaa              >  1:2115202/1‑150 (MQ=255)
                                                                                                                                             aTGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAgcgc          <  1:1708743/150‑1 (MQ=255)
                                                                                                                                              tGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg         <  2:4633368/150‑1 (MQ=255)
                                                                                                                                               gggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCgg        <  2:3017406/150‑1 (MQ=255)
                                                                                                                                                ggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGt       <  2:478822/150‑1 (MQ=255)
                                                                                                                                                ggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGt       <  2:977552/150‑1 (MQ=255)
                                                                                                                                                   gcCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTTGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc    <  1:3876519/150‑1 (MQ=255)
                                                                                                                                                     ccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGa  <  1:2328306/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGA  >  NC_002947/1549706‑1550004

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: