Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,156,379 C→T 6.1% R264C (CGC→TGC)  PP_4539 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,156,3790CT 6.1% 227.4 / 2.9 78R264C (CGC→TGC) PP_4539LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base T (2/3);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GAGCAGCGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGACC  >  NC_002947/5156232‑5156520
                                                                                                                                                   |                                                                                                                                              
gAGCAGCGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAgccgc                                                                                                                                              <  1:253105/150‑1 (MQ=255)
     gcgcTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGccaccc                                                                                                                                         <  1:1115369/150‑1 (MQ=255)
           gACCGTCGGCAGCCTGGGCACCGCCGCCCCCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGtt                                                                                                                                   <  2:233758/150‑1 (MQ=255)
               gTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTAc                                                                                                                               <  1:621641/150‑1 (MQ=255)
                tCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTAcc                                                                                                                              >  2:398389/1‑150 (MQ=255)
                  ggCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTGCTACCtt                                                                                                                            >  2:212027/1‑150 (MQ=255)
                     aGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGACAGCCGCCACCCACTGTTATACCTTTAt                                                                                                                         <  1:931869/150‑1 (MQ=255)
                       ccTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTc                                                                                                                       <  1:1434744/150‑1 (MQ=255)
                              aCCGCCGCCAGCTTCGTCCGCAGCGGCCTGGGCTTTGCCTGGCAGCCCAGTAACATGATCGAGCGCGAGCTGCGCGACGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCTGCCACCCACTGTTCTACCTTTATTCGAGCaaa                                                                                                                <  1:1489189/150‑1 (MQ=255)
                              aCCGCCGCCACCTTCGTCACCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCaaa                                                                                                                <  2:1588771/150‑1 (MQ=255)
                                cgccgcCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAaga                                                                                                              <  1:1387570/150‑1 (MQ=255)
                                cgccgcCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAaga                                                                                                              >  2:496688/1‑150 (MQ=255)
                                 gccgccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAagag                                                                                                             >  1:756621/1‑150 (MQ=255)
                                 gccgccACCTTCGTCAGCAGCGGCCAGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAagag                                                                                                             <  2:2275648/150‑1 (MQ=255)
                                  ccgccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCGGCTGccac                                                                                                                                           <  2:447896/119‑1 (MQ=255)
                                  ccgccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCGGCTGccac                                                                                                                                           >  1:447896/1‑119 (MQ=255)
                                  ccgccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAgaga                                                                                                            >  1:316951/1‑150 (MQ=255)
                                  ccgccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAgaga                                                                                                            >  2:2226187/1‑150 (MQ=255)
                                    gccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCt                                                                                                                             <  2:37231/131‑1 (MQ=255)
                                    gccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCt                                                                                                                             >  1:37231/1‑131 (MQ=255)
                                    gccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCAAGCGATTGGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGagaag                                                                                                          <  1:2349820/150‑1 (MQ=255)
                                     ccACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGagaaga                                                                                                         <  1:1947482/150‑1 (MQ=255)
                                       aCCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGAcc                                                                                                       <  1:765880/150‑1 (MQ=255)
                                        ccTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGAccc                                                                                                      >  2:1154283/1‑150 (MQ=255)
                                          ttCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCGGCTGccaccc                                                                                                                                         >  2:1367176/1‑113 (MQ=255)
                                          ttCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCGGCTGccaccc                                                                                                                                         <  1:1367176/113‑1 (MQ=255)
                                                agcagcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCg                                                                                             <  1:1781997/150‑1 (MQ=255)
                                                 gcagcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCgg                                                                                            >  1:2075536/1‑150 (MQ=255)
                                                 gcagcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCgg                                                                                            <  1:1489537/150‑1 (MQ=255)
                                                  cagcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGc                                                                                           >  1:2059925/1‑150 (MQ=255)
                                                   agcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGcc                                                                                          >  2:1668517/1‑150 (MQ=255)
                                                    gcGGCTTGGGCTTTGCCTGGCTGCCCCGTCAGATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTCCCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCa                                                                                         <  2:1135033/150‑1 (MQ=255)
                                                         cTGGGCTTTGCCTGGCTGCCCCGTCACATGATCAAGCGCGAGCTGCGCGAAGGCGTGCTCAAACCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTATACCTTTATTCGAGCCAAGAGAAGACCCTGGGCCCCGCCACCCAg                                                                                    >  2:154038/1‑150 (MQ=255)
                                                          tGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGa                                                                                   <  2:315638/150‑1 (MQ=255)
                                                               tttGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATtctc                                                                              >  1:368253/1‑150 (MQ=255)
                                                               tttGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATtctc                                                                              >  1:301209/1‑150 (MQ=255)
                                                                 ggCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCAt                                                                            <  2:2255952/149‑1 (MQ=255)
                                                                  gCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATc                                                                           <  1:188068/150‑1 (MQ=255)
                                                                  gCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCACGGCGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATc                                                                           <  1:294905/150‑1 (MQ=255)
                                                                    cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGa                                                                         >  1:1784611/1‑150 (MQ=255)
                                                                    cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGa                                                                         <  2:2057418/150‑1 (MQ=255)
                                                                    cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGa                                                                         >  2:2339423/1‑150 (MQ=255)
                                                                         tGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGa                                                                                                                     >  2:1705458/1‑102 (MQ=255)
                                                                         tGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGa                                                                                                                     <  1:1705458/102‑1 (MQ=255)
                                                                             ccGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCCTTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTgcg                                                                >  1:1025191/1‑150 (MQ=255)
                                                                             ccGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTgcg                                                                >  1:840705/1‑150 (MQ=255)
                                                                             ccGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTgcg                                                                >  1:444990/1‑150 (MQ=255)
                                                                             ccGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGATgcg                                                                >  1:1011834/1‑150 (MQ=255)
                                                                                 cacaTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGc                                                                                      >  2:979007/1‑124 (MQ=255)
                                                                                 cacaTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGc                                                                                      <  1:979007/124‑1 (MQ=255)
                                                                                  acaTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGACGACACTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAAtt                                                           >  1:2338955/1‑150 (MQ=255)
                                                                                  acaTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAAtt                                                           >  2:1746539/1‑150 (MQ=255)
                                                                                   caTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAAttt                                                          <  2:368253/150‑1 (MQ=255)
                                                                                    aTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTc                                                         <  2:165198/150‑1 (MQ=255)
                                                                                    aTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTc                                                         <  1:789694/150‑1 (MQ=255)
                                                                                    aTGATCGAGCGCGAGCTGCGCCAAATCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGTCAGCCGCCACCCACTGTCACACCTTTATTCGAACCCAGAGACATCGTGGTGCTGGGCCACGTAGATTCTCATCGACCTGCTGCGCAATTTc                                                         <  1:421103/150‑1 (MQ=255)
                                                                                      gATCGAGCGCGGGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGa                                                       >  2:2023744/1‑150 (MQ=255)
                                                                                      gATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGa                                                       <  2:2059925/150‑1 (MQ=255)
                                                                                         cGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGacac                                                    >  1:430825/1‑150 (MQ=255)
                                                                                              gcgAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCgcgc                                               <  2:430825/150‑1 (MQ=255)
                                                                                                 aGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCAc                                            <  2:486361/150‑1 (MQ=255)
                                                                                                 aGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCAc                                            <  1:1154283/150‑1 (MQ=255)
                                                                                                  gCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACt                                           <  1:2397925/150‑1 (MQ=255)
                                                                                                         gAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGa                                                       <  2:1450901/131‑1 (MQ=255)
                                                                                                         gAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGa                                                       >  1:1450901/1‑131 (MQ=255)
                                                                                                         gAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTg                                    >  2:212926/1‑150 (MQ=255)
                                                                                                            ggCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGccc                                 <  1:2047686/150‑1 (MQ=255)
                                                                                                            ggCGTGCTCAAGACATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGccc                                 <  1:496688/150‑1 (MQ=255)
                                                                                                               gTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGAAACCGCGCCACTGGACGTGCCCTTc                              >  2:797795/1‑150 (MQ=255)
                                                                                                                  cTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCa                           <  2:2075536/150‑1 (MQ=255)
                                                                                                                     aaGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCtgctgc                                                                 <  1:591377/109‑1 (MQ=255)
                                                                                                                     aaGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCtgctgc                                                                 >  2:591377/1‑109 (MQ=255)
                                                                                                                        ccATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGccccc                     <  2:444990/150‑1 (MQ=255)
                                                                                                                              ccGCTGGATGAGGGTGGCAGCGGCTAGCCACTGTTCTACCTTTATATGAGTAAAGAGAAGAGCCTGGGTCCGGCCATGGAGGTGGTCATCGGCCTGCTGCGCAATTTCGACACGGCGCCACCGGACGTGCCCTTCGCGGCCCCCCcacga               <  2:1357734/150‑3 (MQ=255)
                                                                                                                                      tCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGAc                                                                        >  1:1633071/1‑85 (MQ=255)
                                                                                                                                      tCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGAc                                                                        <  2:1633071/85‑1 (MQ=255)
                                                                                                                                        aGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAgg     >  2:2032753/1‑150 (MQ=255)
                                                                                                                                         gggTGGTAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAACGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGa    >  2:242759/1‑150 (MQ=255)
                                                                                                                                           gTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGAcc  >  1:1051113/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GAGCAGCGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGACC  >  NC_002947/5156232‑5156520

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: