Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,549,816 A→G 8.9% K58K (AAA→AAG)  groL → chaperonin GroEL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,549,8160AG 8.9% 529.6 / 28.8 179K58K (AAA→AAG) groLchaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (70/93);  new base G (11/5);  total (81/98)
Fisher's exact test for biased strand distribution p-value = 6.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 186 total aligned reads displayed.

GGTGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTC  >  NC_002947/1549667‑1549963
                                                                                                                                                     |                                                                                                                                                   
ggTGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCaaa                                                                                                                                                     >  2:455438/1‑150 (MQ=255)
ggTGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTAAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCaaa                                                                                                                                                     <  1:2597224/150‑1 (MQ=255)
   gATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGaa                                                                                                                                                  >  2:2855697/1‑150 (MQ=255)
    aTTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGaaa                                                                                                                                                 <  2:5158207/150‑1 (MQ=255)
      tCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATc                                                                                                                                               >  1:2234334/1‑150 (MQ=255)
        cGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGa                                                                                                                                             <  1:2189931/150‑1 (MQ=255)
         gCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAg                                                                                                                                            <  1:5381635/150‑1 (MQ=255)
           tCGTAAGAAAATGCTGGTTGGTGACAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCt                                                                                                                                          >  2:4672526/1‑150 (MQ=255)
              tAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGaa                                                                                                                                       <  1:4775279/150‑1 (MQ=255)
                            ttggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGaaa                                                                                                                         <  2:412519/150‑1 (MQ=255)
                            ttggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGaaa                                                                                                                         <  2:1914139/150‑1 (MQ=255)
                              ggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAAc                                                                                                                       <  2:3809530/150‑1 (MQ=255)
                                 gtCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATg                                                                                                                    >  2:2199491/1‑150 (MQ=255)
                                 gtCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATg                                                                                                                    >  2:3621565/1‑150 (MQ=255)
                                  tCAAGGTTGGGGCTGAGGCGGTAAATGCGACCCTGGGCCCGAACGGCCGTAACGTGGTACTGGCCAAGAGCTTGGGCGCGCGGACCATCACCACGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                   <  2:4984089/150‑1 (MQ=255)
                                  tCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                   <  2:3087618/150‑1 (MQ=255)
                                  tCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCGTTCGAAAACATgg                                                                                                                   >  1:216537/1‑150 (MQ=255)
                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGc                                                                                                                                 >  1:3730827/1‑134 (MQ=255)
                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTc                                                                                                                             >  2:5335130/1‑138 (MQ=255)
                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGgc                                                                                                                 >  1:3339704/1‑150 (MQ=255)
                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGgc                                                                                                                 >  2:2958613/1‑150 (MQ=255)
                                            ggCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCt                                                                                                         <  2:2295048/150‑1 (MQ=255)
                                             gCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTg                                                                                                        >  2:2999969/1‑150 (MQ=255)
                                              cTGACGCGGTAACAGCGACCCGGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTgg                                                                                                       <  2:3912215/150‑1 (MQ=255)
                                               tGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGt                                                                                                      >  1:4685471/1‑150 (MQ=255)
                                               tGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGt                                                                                                      <  2:4739201/150‑1 (MQ=255)
                                               tGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGt                                                                                                      <  1:2999969/150‑1 (MQ=255)
                                               cgACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTGCGGCGCGCCGACCATCACCAATGACGGCTTTTCCGTCGCCAAAGAAATCGAGCTGATATACGCTTTAGAATACTTGGTAGGCCTGCTGGt                                                                                                      <  2:3883571/149‑1 (MQ=255)
                                                   gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                  <  1:2115988/150‑1 (MQ=255)
                                                   gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                  >  1:1730879/1‑150 (MQ=255)
                                                   gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                  <  1:2184222/150‑1 (MQ=255)
                                                   gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                  >  1:4078260/1‑150 (MQ=255)
                                                      gTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGaa                                                                                               >  2:5373765/1‑150 (MQ=255)
                                                         aaaGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGtt                                                                                            <  1:3450256/150‑1 (MQ=255)
                                                           aGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGc                                                                                          >  2:2087441/1‑150 (MQ=255)
                                                              gACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTc                                                                                       <  2:5182001/150‑1 (MQ=255)
                                                                 ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGTTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                    >  2:3614990/1‑150 (MQ=255)
                                                                 ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                    >  2:1492818/1‑150 (MQ=255)
                                                                   tGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAgg                                                                                  >  2:3451200/1‑150 (MQ=255)
                                                                    gggCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAACGACGCTGTCGAAATCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGc                                                                                 >  1:5289559/1‑150 (MQ=255)
                                                                       cccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGa                                                                                                                                             >  1:1249326/1‑87 (MQ=255)
                                                                       cccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGa                                                                                                                                             <  2:1249326/87‑1 (MQ=255)
                                                                        ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                             >  1:4015804/1‑150 (MQ=255)
                                                                        ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                             <  1:455438/150‑1 (MQ=255)
                                                                         cGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacg                                                                            >  1:2341306/1‑150 (MQ=255)
                                                                            aaGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgacg                                                                         <  2:2318470/150‑1 (MQ=255)
                                                                            aaGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgacg                                                                         <  1:1066322/150‑1 (MQ=255)
                                                                                ccGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGc                                                                     <  1:960279/150‑1 (MQ=255)
                                                                                ccGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGc                                                                     >  2:782928/1‑150 (MQ=255)
                                                                                 cGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGcc                                                                    <  2:3333790/150‑1 (MQ=255)
                                                                                  gTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCg                                                                   <  1:1824792/150‑1 (MQ=255)
                                                                                  gTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCg                                                                   <  2:4389113/150‑1 (MQ=255)
                                                                                   tAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCgg                                                                  >  2:3414326/1‑150 (MQ=255)
                                                                                        tGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACg                                                             >  2:1232257/1‑150 (MQ=255)
                                                                                          gTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGc                                                           <  1:3403398/150‑1 (MQ=255)
                                                                                              tGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcacca                                                       >  1:2720040/1‑150 (MQ=255)
                                                                                                  cAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcaccaccac                                                   <  1:4171823/150‑1 (MQ=255)
                                                                                                       gCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTa                                              <  1:1437285/150‑1 (MQ=255)
                                                                                                       ccTTCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGcga                                              >  1:4034462/2‑148 (MQ=255)
                                                                                                        cTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAc                                             >  1:2369894/1‑150 (MQ=255)
                                                                                                        cTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAc                                             >  1:569454/1‑150 (MQ=255)
                                                                                                         ttCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGAcc                                            >  2:5278307/1‑150 (MQ=255)
                                                                                                         ttCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAcc                                            >  2:1416183/1‑150 (MQ=255)
                                                                                                         ttCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAcc                                            >  1:1352239/1‑150 (MQ=255)
                                                                                                            ggCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGtt                                         >  1:3990380/1‑150 (MQ=255)
                                                                                                             gcgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTc                                        <  2:2231213/150‑1 (MQ=255)
                                                                                                              cgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                       <  1:3607078/150‑1 (MQ=255)
                                                                                                               gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCg                                                                                                                            >  2:4962963/1‑64 (MQ=255)
                                                                                                               gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCg                                                                                                                            <  1:4962963/64‑1 (MQ=255)
                                                                                                               gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTg                                      >  2:1002034/1‑150 (MQ=255)
                                                                                                                cgcCGACCATCACCCAGGACGGCGTTTCCGTCGCCAAAGAACTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                     <  2:516871/150‑1 (MQ=255)
                                                                                                                    gACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCa                                 <  1:1867578/150‑1 (MQ=255)
                                                                                                                     aCCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAg                                <  2:4189310/150‑1 (MQ=255)
                                                                                                                     aCCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAg                                <  1:5339013/150‑1 (MQ=255)
                                                                                                                      gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                               >  1:4828461/2‑150 (MQ=255)
                                                                                                                      ccATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAgg                               <  2:4027879/150‑1 (MQ=255)
                                                                                                                      ccATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAgg                               <  1:184726/150‑1 (MQ=255)
                                                                                                                           aCCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATc                          <  1:3305332/150‑1 (MQ=255)
                                                                                                                           aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                          >  2:1390641/1‑150 (MQ=255)
                                                                                                                           aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                          >  1:274733/1‑150 (MQ=255)
                                                                                                                           aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                          <  1:1967786/150‑1 (MQ=255)
                                                                                                                           aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                          <  2:2783581/150‑1 (MQ=255)
                                                                                                                               aGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCa                      <  2:1730879/150‑1 (MQ=255)
                                                                                                                                ggACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttca                     <  2:1493699/150‑4 (MQ=255)
                                                                                                                                 gACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttcag                    >  2:4901322/1‑146 (MQ=255)
                                                                                                                                 gACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttcag                    >  1:1087694/1‑146 (MQ=255)
                                                                                                                                    ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                 >  1:1465819/1‑150 (MQ=255)
                                                                                                                                     gCGTTTCCGTCGCCAAAGAAATCGAGCTGAACGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                <  2:1310292/150‑1 (MQ=255)
                                                                                                                                         ttCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt            <  2:5018223/150‑1 (MQ=255)
                                                                                                                                          tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAACACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg           <  2:842396/150‑1 (MQ=255)
                                                                                                                                          tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg           >  1:4458275/1‑150 (MQ=255)
                                                                                                                                           ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa          <  1:2738230/150‑1 (MQ=255)
                                                                                                                                            cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGag         <  1:2855697/150‑2 (MQ=255)
                                                                                                                                             gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                       >  2:3210976/1‑119 (MQ=255)
                                                                                                                                             gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                       <  1:3210976/119‑1 (MQ=255)
                                                                                                                                             gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa        <  1:4020978/150‑1 (MQ=255)
                                                                                                                                             gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa        <  1:2074508/150‑1 (MQ=255)
                                                                                                                                                  cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCGCCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt   >  2:4522324/1‑150 (MQ=255)
                                                                                                                                                  cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt   >  1:4284876/1‑150 (MQ=255)
                                                                                                                                                   aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc  <  1:4070381/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTC  >  NC_002947/1549667‑1549963

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: ‑