Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,966,807 T→G 16.1% L319R (CTG→CGG) ‡ ydcP → peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,966,8070TG16.1% 82.4 / 9.1 39L319R (CTG→CGG) ‡ydcPpeptidase
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base G (2/4);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCG  >  NC_002947/3966664‑3966950
                                                                                                                                                |                                                                                                                                               
gTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGt                                                                                                                                           >  2:218562/1‑150 (MQ=255)
    gCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTgaag                                                                                                                                       <  1:510635/150‑4 (MQ=255)
    gCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTgaag                                                                                                                                       <  1:345803/150‑4 (MQ=255)
      ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTgca                                                                                                                                        <  2:429181/147‑3 (MQ=255)
      ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTgaa                                                                                                                                        >  1:429181/1‑145 (MQ=255)
             tCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaa                                                                                                                              >  2:925/1‑150 (MQ=255)
                    gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg                                                                                                                       <  1:925/150‑1 (MQ=255)
                    gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg                                                                                                                       <  1:333447/150‑1 (MQ=255)
                             aaCACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGCCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTg                                                                                                              >  1:128731/1‑150 (MQ=255)
                              acacTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGa                                                                                                             >  1:126774/1‑150 (MQ=255)
                                  ttcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGt                                                                                                               <  2:56212/144‑1 (MQ=255)
                                  ttcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGt                                                                                                               >  1:56212/1‑144 (MQ=255)
                                      tcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAg                                                                                                     <  1:176632/150‑1 (MQ=255)
                                       ccTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGt                                                                                                    >  1:269114/1‑150 (MQ=255)
                                           cccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTc                                                                                                <  1:486613/150‑1 (MQ=255)
                                            ccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAg                                                                                                                           >  1:95523/1‑122 (MQ=255)
                                            ccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAg                                                                                                                           <  2:95523/122‑1 (MQ=255)
                                                   ggATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAg                                                                                                     >  2:94987/1‑137 (MQ=255)
                                                   ggATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAg                                                                                                     <  1:94987/137‑1 (MQ=255)
                                                         aaCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCt                                                                                  >  1:78533/1‑150 (MQ=255)
                                                         aaCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCt                                                                                  >  2:527282/1‑150 (MQ=255)
                                                         aaCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACGCGCCAACCCACGCCACTCCAGCAAGGTTACGTTGGTAAAAGGTCGGCAAGCGTGACATGCAAATGGACACCGACGTACCGCt                                                                                  >  1:258407/1‑150 (MQ=255)
                                                           ccTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGa                                                                                <  1:218562/150‑1 (MQ=255)
                                                                cttctcgGCAGCACCGTCTACTTTGTCACCGCACGCAAGTTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGTCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAAc                                                                           <  2:60520/145‑1 (MQ=255)
                                                                               gACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAAt                                                            <  2:269114/150‑1 (MQ=255)
                                                                                 ctactTCGTGACCGACCGCAAGGCCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAAtgt                                                          <  2:128731/150‑1 (MQ=255)
                                                                                     ttCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCAGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTg                                                      >  1:440622/1‑150 (MQ=255)
                                                                                        gTGACCGACCGTAAGGTCGATATCGTCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTc                                                   <  1:561801/150‑1 (MQ=255)
                                                                                          gaccgaccGCAAGGTCGATATCGGCGCCTTCGGCTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCaa                                                 <  1:527282/150‑1 (MQ=255)
                                                                                          gaccgaccGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCaa                                                 >  2:121953/1‑150 (MQ=255)
                                                                                                       gTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGc                                                                                                                 >  1:632692/1‑73 (MQ=255)
                                                                                                       gTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGc                                                                                                                 <  2:632692/73‑1 (MQ=255)
                                                                                                                                   ccTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt                                                       <  2:25356/103‑1 (MQ=255)
                                                                                                                                   ccTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt                                                       >  1:25356/1‑103 (MQ=255)
                                                                                                                                      ccACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCACAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACcgcgcg     >  1:24700/2‑148 (MQ=255)
                                                                                                                                        aCCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt                                                       >  2:643751/1‑98 (MQ=255)
                                                                                                                                        aCCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt                                                       <  1:643751/98‑1 (MQ=255)
                                                                                                                                        aCCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGc   <  2:78533/150‑1 (MQ=255)
                                                                                                                                         ccGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCACAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGCGGca  <  2:24700/150‑2 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCG  >  NC_002947/3966664‑3966950

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: