Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,951,015 T→G 6.6% E666D (GAA→GAC ppk ← polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,951,0150TG 6.6% 163.0 / 2.7 60E666D (GAA→GACppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base G (2/2);  total (28/32)
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.18e-01

AGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCA  >  NC_002947/5950866‑5951157
                                                                                                                                                     |                                                                                                                                              
agcagGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCt                                                                                                                                                <  2:634120/150‑1 (MQ=255)
agcagGGTCGCCTGAACGCTGCGCGGGTTCTCGTTGCCGCTCGGCGTGCTGCGCACGTAGCGCCCGTCTGACTTCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCt                                                                                                                                                <  1:1048334/150‑1 (MQ=255)
      gttgCCTGAACGCTGCGCGGGTTTTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCAcc                                                                                                                                          <  1:273254/147‑1 (MQ=255)
        cGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCgg                                                                                                                                        <  1:490846/150‑1 (MQ=255)
        cGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCgg                                                                                                                                        <  1:490747/150‑1 (MQ=255)
           cTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGaa                                                                                                                                     >  2:1270774/1‑150 (MQ=255)
            tGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAg                                                                                                                                    <  2:335640/150‑1 (MQ=255)
             gAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGc                                                                                                                                   >  1:255358/1‑150 (MQ=255)
              aaCGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCa                                                                                                                                  <  2:1250030/150‑1 (MQ=255)
              aaCGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCa                                                                                                                                  >  1:1078671/1‑150 (MQ=255)
                         ggTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGacac                                                                                                                       >  2:770863/1‑150 (MQ=255)
                         ggTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGacac                                                                                                                       <  2:936399/150‑1 (MQ=255)
                          gTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACg                                                                                                                      <  1:152906/150‑1 (MQ=255)
                               ggTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGt                                                                                                                 <  1:770863/150‑1 (MQ=255)
                                     cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCCGAGCTACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGt                                                                                                           <  2:816220/150‑1 (MQ=255)
                                      ggTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          <  1:613458/150‑1 (MQ=255)
                                      ggTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          >  1:564925/1‑150 (MQ=255)
                                          ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCTTTACCACGGTAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                      <  1:1210471/150‑1 (MQ=255)
                                                ctgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatc                                                                                                <  1:791142/150‑1 (MQ=255)
                                                        gTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgc                                                                                        >  2:1168785/1‑150 (MQ=255)
                                                        gTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgc                                                                                        >  2:1168732/1‑150 (MQ=255)
                                                           gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCt                                                                                     <  1:1168785/150‑1 (MQ=255)
                                                           gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCt                                                                                     <  1:1168732/150‑1 (MQ=255)
                                                              cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          >  2:355921/1‑126 (MQ=255)
                                                              cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          <  1:355921/126‑1 (MQ=255)
                                                                  tCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGCTCCACAGATCGGAAGGGGAAGCATGTCTCGACACGCATGTCGCCGGACGGCTCTTTCCAGTCGGCGCTGGACAgg                                                                              >  1:665385/1‑150 (MQ=255)
                                                                   cTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGt                                                                             <  1:1177429/150‑1 (MQ=255)
                                                                     ggCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAg                                                                           <  1:1270774/150‑1 (MQ=255)
                                                                            aGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGACGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGct                                                                    <  2:922599/150‑1 (MQ=255)
                                                                            aGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGct                                                                    <  2:624407/150‑1 (MQ=255)
                                                                              gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                     <  1:303552/99‑1 (MQ=255)
                                                                              gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                     >  2:303552/1‑99 (MQ=255)
                                                                                           tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                     <  2:564925/150‑1 (MQ=255)
                                                                                           tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATATGCTCTTCGCCGCCATTGa                                                     >  2:733923/1‑150 (MQ=255)
                                                                                            aTTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAg                                                    <  2:1078671/150‑1 (MQ=255)
                                                                                                   gTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAg                                             >  1:659655/1‑150 (MQ=255)
                                                                                                    tCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTACTTGCATTCCACCAGGAAGCAGGTCTGGGCAGGCTTATCGACAGAAGTCGCCAACCGGTCGCCGCGGTACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGa                                            >  2:786121/1‑150 (MQ=255)
                                                                                                       ggTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGt                                               <  2:1180625/144‑1 (MQ=255)
                                                                                                       ggTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGt                                               >  1:1180625/1‑144 (MQ=255)
                                                                                                                 cAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCCCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGtctg                               >  1:67412/1‑150 (MQ=255)
                                                                                                                  aGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   >  2:1152669/1‑81 (MQ=255)
                                                                                                                  aGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   <  1:1152669/81‑1 (MQ=255)
                                                                                                                   gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   <  2:1052902/80‑1 (MQ=255)
                                                                                                                   gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   >  1:1052902/1‑80 (MQ=255)
                                                                                                                   gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCt                             >  1:1239969/1‑150 (MQ=255)
                                                                                                                      ccTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                          <  1:1189357/150‑1 (MQ=255)
                                                                                                                           ttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGATTGCGCTCGATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGTCGCCATTGAGGACGTAGAACACCCGCGTGTGCTCGAGGAAg                     >  1:998929/1‑150 (MQ=255)
                                                                                                                             cACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCg                   >  2:606225/1‑150 (MQ=255)
                                                                                                                                ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                 >  2:995122/1‑85 (MQ=255)
                                                                                                                                ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                 <  1:995122/85‑1 (MQ=255)
                                                                                                                                  gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAaca                                         >  2:72513/1‑123 (MQ=255)
                                                                                                                                  gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAaca                                         <  1:72513/123‑1 (MQ=255)
                                                                                                                                  gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCaa              <  2:659655/150‑1 (MQ=255)
                                                                                                                                   cagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                <  2:277298/147‑1 (MQ=255)
                                                                                                                                   cagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                >  1:277298/1‑147 (MQ=255)
                                                                                                                                       aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTTGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc         >  1:380272/1‑150 (MQ=255)
                                                                                                                                       aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc         >  2:1000241/1‑150 (MQ=255)
                                                                                                                                       aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGATGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc         <  2:255358/150‑1 (MQ=255)
                                                                                                                                        aCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCg        >  1:810533/1‑150 (MQ=255)
                                                                                                                                              tCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAACTGCCTTCTACGCCGCCATTGAGGAAGTAGAACACCAGCGTGTGCTCGATGAAGCGGCCAATGATCGAGCGCa  >  1:988189/1‑150 (MQ=255)
                                                                                                                                                  gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACAcc                                       >  2:1270486/1‑109 (MQ=255)
                                                                                                                                                  gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACAcc                                       <  1:1270486/109‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCA  >  NC_002947/5950866‑5951157

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