Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,549,859 C→A 13.0% L73M (CTG→ATG)  groL → chaperonin GroEL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,549,8590CA13.0% 580.7 / 63.6 200L73M (CTG→ATG) groLchaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (80/94);  new base A (18/8);  total (98/102)
Fisher's exact test for biased strand distribution p-value = 3.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 201 total aligned reads displayed.

GGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAG  >  NC_002947/1549711‑1550008
                                                                                                                                                    |                                                                                                                                                     
ggCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCt                                                                                                                                                      <  2:2295048/150‑1 (MQ=255)
 gCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTg                                                                                                                                                     >  2:2999969/1‑150 (MQ=255)
   tGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGt                                                                                                                                                   <  2:4739201/150‑1 (MQ=255)
   cgACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTGCGGCGCGCCGACCATCACCAATGACGGCTTTTCCGTCGCCAAAGAAATCGAGCTGATATACGCTTTAGAATACTTGGTAGGCCTGCTGGt                                                                                                                                                   <  2:3883571/149‑1 (MQ=255)
       gcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAg                                                                                                                                               <  1:2184222/150‑1 (MQ=255)
           tAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAg                                                                                                                                           <  2:1877224/150‑1 (MQ=255)
                  gACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTc                                                                                                                                    <  2:5182001/150‑1 (MQ=255)
                     ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGTTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                                                                 >  2:3614934/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)
                           cccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCaa                                                                                                                           <  1:2130657/150‑1 (MQ=255)
                            ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                                                                          <  2:1335135/150‑1 (MQ=255)
                            ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                                                                          <  1:455438/150‑1 (MQ=255)
                            ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                                                                          <  1:272800/150‑1 (MQ=255)
                            ccGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAac                                                                                                                          >  1:4015804/1‑150 (MQ=255)
                             cGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacg                                                                                                                         >  1:2341306/1‑150 (MQ=255)
                                aaGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgacg                                                                                                                      <  1:1066322/150‑1 (MQ=255)
                                     cGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGcc                                                                                                                 <  2:3333790/150‑1 (MQ=255)
                                       tAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCgg                                                                                                               <  2:4469905/150‑1 (MQ=255)
                                       tAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCgg                                                                                                               >  2:3414326/1‑150 (MQ=255)
                                            tGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACg                                                                                                          >  2:1232408/1‑150 (MQ=255)
                                                 cTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcacc                                                                                                     <  2:4987752/150‑1 (MQ=255)
                                                  tGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcacca                                                                                                    >  1:3383552/1‑150 (MQ=255)
                                                          agCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCt                                                                                            <  1:4345416/150‑1 (MQ=255)
                                                           gCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTa                                                                                           <  1:1437285/150‑1 (MQ=255)
                                                           ccTTCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGcga                                                                                           <  1:261480/149‑3 (MQ=255)
                                                           ccTTCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGcga                                                                                           >  1:4034462/2‑148 (MQ=255)
                                                            cTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAc                                                                                          >  1:2369894/1‑150 (MQ=255)
                                                             ttCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGAcc                                                                                         >  2:5341984/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)
                                                                  cgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCGAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                                                                    <  2:4049156/150‑1 (MQ=255)
                                                                  cgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                                                                    <  1:414080/150‑1 (MQ=255)
                                                                  cgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                                                                    <  1:3607028/150‑1 (MQ=255)
                                                                   gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTg                                                                                   >  2:1002034/1‑150 (MQ=255)
                                                                    cgcCGACCATCACCCAGGACGGCGTTTCCGTCGCCAAAGAACTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                                                                  <  2:516871/150‑1 (MQ=255)
                                                                     gcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGc                                                                                 >  1:2908653/1‑150 (MQ=255)
                                                                      ccGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgac                                                                                                                       >  2:4547798/1‑111 (MQ=255)
                                                                        gACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCa                                                                              <  1:1867578/150‑1 (MQ=255)
                                                                        acgcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCa                                                                              >  1:3685903/4‑150 (MQ=255)
                                                                         aCCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAg                                                                             <  2:4189310/150‑1 (MQ=255)
                                                                          gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                                                            >  1:2093184/2‑150 (MQ=255)
                                                                          gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                                                            >  1:4828461/2‑150 (MQ=255)
                                                                          ccATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAgg                                                                            <  1:184726/150‑1 (MQ=255)
                                                                             tCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCa                                                                         >  2:1707257/1‑150 (MQ=255)
                                                                               aCCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATc                                                                       <  1:3305332/150‑1 (MQ=255)
                                                                               aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                       <  1:1967786/150‑1 (MQ=255)
                                                                               aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                       <  2:2783581/150‑1 (MQ=255)
                                                                                  aaGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                                                                    >  2:2556282/1‑134 (MQ=255)
                                                                                  aaGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTc                                                                    >  1:4073065/1‑150 (MQ=255)
                                                                                        ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                                                              >  1:1465819/1‑150 (MQ=255)
                                                                                              tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAACACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg                                                        <  2:842396/150‑1 (MQ=255)
                                                                                              tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg                                                        >  1:4458275/1‑150 (MQ=255)
                                                                                                 gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCt                                                                                    <  1:3210976/119‑1 (MQ=255)
                                                                                                 gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                     <  1:4020978/150‑1 (MQ=255)
                                                                                                    gCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACAACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGcc                                                  >  1:2219153/1‑150 (MQ=255)
                                                                                                      gcaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                <  2:4907492/148‑1 (MQ=255)
                                                                                                      cAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCtc                                                                                   <  2:537774/115‑2 (MQ=255)
                                                                                                      cAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCtc                                                                                   >  1:537774/1‑114 (MQ=255)
                                                                                                      cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                >  1:4284876/1‑150 (MQ=255)
                                                                                                           aaaTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTg                                           >  2:4402455/1‑150 (MQ=255)
                                                                                                                gAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGc                                      >  2:4107260/1‑150 (MQ=255)
                                                                                                                 aGCTGAAAGACGCTTTCGAAAACATGGGCGCCCATCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGTTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCAATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCa                                     <  1:4672526/150‑1 (MQ=255)
                                                                                                                  gCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCAt                                    <  2:2931253/150‑1 (MQ=255)
                                                                                                                   cTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATg                                   >  2:5045400/1‑150 (MQ=255)
                                                                                                                    tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                  >  2:4323395/1‑150 (MQ=255)
                                                                                                                     gAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGaa                                 >  1:3952303/1‑150 (MQ=255)
                                                                                                                     gAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGaa                                 >  1:1805353/1‑150 (MQ=255)
                                                                                                                         gACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCg                             >  2:3768277/1‑150 (MQ=255)
                                                                                                                          acaagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAAcccaa                            >  1:54209/6‑148 (MQ=255)
                                                                                                                          acaagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAAcccaa                            >  1:4172471/6‑148 (MQ=255)
                                                                                                                          aCGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                            >  2:782035/1‑150 (MQ=255)
                                                                                                                           caagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAAt                           <  2:4702801/146‑1 (MQ=255)
                                                                                                                            gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                          <  1:782035/150‑1 (MQ=255)
                                                                                                                            gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                          >  2:3197772/1‑150 (MQ=255)
                                                                                                                            gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                          <  2:2341306/150‑1 (MQ=255)
                                                                                                                            gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                          <  1:2744422/150‑1 (MQ=255)
                                                                                                                            gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                          <  1:55358/150‑1 (MQ=255)
                                                                                                                             agttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATgg                         <  1:4775466/148‑1 (MQ=255)
                                                                                                                              tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGa                        >  2:2863731/1‑150 (MQ=255)
                                                                                                                               ttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAc                       <  2:2131630/150‑1 (MQ=255)
                                                                                                                               ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGaa                                 >  1:1174108/1‑140 (MQ=255)
                                                                                                                               ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGaa                                 <  2:1174108/140‑1 (MQ=255)
                                                                                                                               ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGAc                       >  1:2847167/1‑150 (MQ=255)
                                                                                                                                tCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAcc                      >  2:2087280/1‑150 (MQ=255)
                                                                                                                                  gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                    >  1:1246359/1‑150 (MQ=255)
                                                                                                                                  gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAGGGACCTg                    >  1:932950/1‑150 (MQ=255)
                                                                                                                                  gAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGACCTg                    >  2:1624040/1‑150 (MQ=255)
                                                                                                                                   aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGa                   >  1:2049096/1‑150 (MQ=255)
                                                                                                                                    aaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                  >  2:4595444/1‑150 (MQ=255)
                                                                                                                                    aaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                  <  1:3451200/150‑1 (MQ=255)
                                                                                                                                    aaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                  >  1:3161241/1‑150 (MQ=255)
                                                                                                                                         tGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg             <  1:3414326/150‑1 (MQ=255)
                                                                                                                                         tGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg             <  2:3971086/150‑1 (MQ=255)
                                                                                                                                          gggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCgg            >  2:1793645/1‑150 (MQ=255)
                                                                                                                                              gcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc        >  1:1522263/1‑150 (MQ=255)
                                                                                                                                               cccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCg       >  2:4883305/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAG  >  NC_002947/1549711‑1550008

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