Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,481 G→A 6.6% V200M (GTG→ATG)  ttuD → hydroxypyruvate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,4810GA 6.6% 169.5 / 2.7 61V200M (GTG→ATG) ttuDhydroxypyruvate reductase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

CATCGGCGAGATGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGACCCGGCCCTGGCCCGCAGCCACTT  >  NC_002947/4891333‑4891626
                                                                                                                                                    |                                                                                                                                                 
caTCGGCGAGATGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGt                                                                                                                                                  <  2:135050/150‑1 (MQ=255)
caTCGGCGAGATGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGt                                                                                                                                                  <  2:743399/150‑1 (MQ=255)
 aTCGGCGAGATGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTg                                                                                                                                                 <  2:28100/150‑1 (MQ=255)
      cGAGATGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGa                                                                                                                                            >  1:33121/1‑150 (MQ=255)
          aTGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCg                                                                                                                                        >  2:709539/1‑150 (MQ=255)
            gAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCGACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAg                                                                                                                                      <  2:314172/150‑1 (MQ=255)
            gAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAg                                                                                                                                      >  2:82352/1‑150 (MQ=255)
                    tgcgCAAGCACCTCTCGGCGATCAAGATCGAGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCgg                                                                                                                              >  1:759549/1‑150 (MQ=255)
                         aaGCACCTCTCGGCGATCAAGGGCAGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGTCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGAc                                                                                                                         <  2:319160/150‑1 (MQ=255)
                            cACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGcc                                                                                                                      <  1:787322/150‑1 (MQ=255)
                               ctctCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTg                                                                                                                   >  1:485464/1‑150 (MQ=255)
                                tctcGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTgg                                                                                                                  <  1:82352/150‑1 (MQ=255)
                                 ctcGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCGGGCCGGCGACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCCCCGGCCCCACCGTGGCCGACCCGAGCTCCTCGTCCGAGGCCCTGGc                                                                                                                 <  2:444210/150‑1 (MQ=255)
                                    ggCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCAt                                                                                                              >  1:394108/1‑150 (MQ=255)
                                     gCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCttgatcgc                                                                                                             >  1:232032/1‑143 (MQ=255)
                                     gCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATc                                                                                                             <  1:500638/150‑1 (MQ=255)
                                     gCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATc                                                                                                             <  2:1012877/150‑1 (MQ=255)
                                       gATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCt                                                                                                           >  1:276139/1‑150 (MQ=255)
                                          cAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCGTCCTCaa                                                                                                        >  1:121833/1‑150 (MQ=255)
                                            aGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCGCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAAc                                                                                                      <  1:350501/150‑1 (MQ=255)
                                                cgggcgCCTGGCCAGGCCCGGCTCGCCGGCCCCGGTCTACACCTATGCCCTTTCCGAGGTACCGCGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGCCCGACCCGCGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTa                                                                                                  <  2:234992/150‑1 (MQ=255)
                                                cgggcgCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTa                                                                                                  <  1:389006/150‑1 (MQ=255)
                                                 gggcgCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTac                                                                                                 >  1:221405/1‑150 (MQ=255)
                                                       ctgaacAAGGCCTGCTAGGCGGCCACTGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATc                                                                                           >  2:1123004/6‑150 (MQ=255)
                                                       cTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATc                                                                                           >  1:1100390/1‑150 (MQ=255)
                                                          gCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGaa                                                                                        <  2:485464/150‑1 (MQ=255)
                                                           ccAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAg                                                                                       >  1:480603/1‑150 (MQ=255)
                                                            cAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAgc                                                                                      <  2:132943/150‑1 (MQ=255)
                                                             aaGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAgcg                                                                                     >  2:950840/1‑150 (MQ=255)
                                                             aaGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAgcg                                                                                     >  1:337704/1‑150 (MQ=255)
                                                             aaGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAgcg                                                                                     <  1:267229/150‑1 (MQ=255)
                                                                 ccTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGAAGGACCGGGCGACCTCGCCACGGTAATCCCCTCCGGCACCACCGTGGCTGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACAACGAAGCGCCCa                                                                                 >  1:630099/1‑150 (MQ=255)
                                                                   tgctgGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCaaa                                                                               >  1:638011/1‑150 (MQ=255)
                                                                          cggccACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGCCCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCa                                                                        >  1:730432/1‑150 (MQ=255)
                                                                          cggccACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCa                                                                        <  1:1123004/150‑1 (MQ=255)
                                                                          cggccACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCa                                                                        <  2:221405/150‑1 (MQ=255)
                                                                            gccACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATc                                                                      <  2:638011/150‑1 (MQ=255)
                                                                              cACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGa                                                                    <  2:480603/150‑1 (MQ=255)
                                                                               aCGGTCTACACATATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGAc                                                                   >  2:910625/1‑150 (MQ=255)
                                                                                           tATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGaacaac                                                       >  1:369310/1‑150 (MQ=255)
                                                                                                cATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGc                                                  <  2:1100390/150‑1 (MQ=255)
                                                                                                      cGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGa                                            >  2:349968/1‑150 (MQ=255)
                                                                                                         tgccccGCGCGCCCTCGCCACGCTAACCGCCTCCGTCCCCACCGTGGGCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAACCGCTACAACATCGAAGCGCCCAAAGCGGTCATAGACTGGCTGAACAACCCGGCCTCGGAAAc                                         <  2:121833/146‑1 (MQ=255)
                                                                                                             ccGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAAccc                                                     <  2:839633/134‑1 (MQ=255)
                                                                                                             ccGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAAccc                                                     >  1:839633/1‑134 (MQ=255)
                                                                                                                ggCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAg                                  >  1:386699/1‑150 (MQ=255)
                                                                                                                   gACCTCGCCACGGTAATCGCCTCCGGAGCCGTCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGcc                               >  1:578843/1‑150 (MQ=255)
                                                                                                                     ccTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGCGAACTCATGGCGTCTCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGa                             <  2:578843/150‑1 (MQ=255)
                                                                                                                        cGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGa                          <  2:276139/150‑1 (MQ=255)
                                                                                                                            aCGGTAATCGCCTCCGGCCCCACCATGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAg                                  >  1:476630/1‑138 (MQ=255)
                                                                                                                            aCGGTAATCGCCTCCGGCCCCACCATGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAg                                  >  1:476780/1‑138 (MQ=255)
                                                                                                                            aCGGTAATCGCCTCCGGCCCCACCATGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAg                                  <  2:476630/138‑1 (MQ=255)
                                                                                                                            aCGGTAATCGCCTCCGGCCCCACCATGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAg                                  <  2:476780/138‑1 (MQ=255)
                                                                                                                                 aaTCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAGTCCTTGCTCGGAAACCATCAAGGCCTGGGACGCGGCCCt                 >  1:57626/1‑150 (MQ=255)
                                                                                                                                      ccTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGACCCGGCCCTGGccc            >  1:84465/1‑150 (MQ=255)
                                                                                                                                            gCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGACCCGGCCCTGGCCCGCAGcc      <  1:697423/150‑1 (MQ=255)
                                                                                                                                              cccACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAgg                                 <  2:1085320/121‑1 (MQ=255)
                                                                                                                                              cccACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAgg                                 >  1:1085320/1‑121 (MQ=255)
                                                                                                                                                cACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGACCCGGCCCTGGCCCGCAGCCACtt  <  2:655791/150‑1 (MQ=255)
                                                                                                                                                  ccGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCACCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGa                          <  2:140288/124‑1 (MQ=255)
                                                                                                                                                  ccGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCACCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGa                          >  1:140288/1‑124 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATCGGCGAGATGAACTGCGTGCGCAAGCACCTCTCGGCGATCAAGGGCGGGCGCCTGGCCAAGGCCTGCTGGCCGGCCACGGTCTACACCTATGCCATTTCCGATGTACCGGGCGACCTCGCCACGGTAATCGCCTCCGGCCCCACCGTGGCCGACCCGAGCACCTCGGCCGACGCCCTGGCCATCCTCAAACGCTACAACATCGAAGCGCCCAAAGCGGTCATCGACTGGCTGAACAACCCGGCCTCGGAAACCGTCAAGGCCGATGACCCGGCCCTGGCCCGCAGCCACTT  >  NC_002947/4891333‑4891626

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