Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,027,999 G→C 16.6% A248A (GCC→GCG dppA‑III ← dipeptide ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,027,9990GC16.6% 51.8 / 2.4 23A248A (GCC→GCGdppA‑IIIdipeptide ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base C (2/2);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACGTTGTAGGCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGATCGGCTTCTGGT  >  NC_002947/1027850‑1028141
                                                                                                                                                     |                                                                                                                                              
aTCACGTTGTAGGCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAg                                                                                                                                                >  2:427264/1‑150 (MQ=255)
  cACGTCGTAGGCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCCCGATGTCGGCCGGGCGTGGTAACAGGGTGACCTGCCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGc                                                                                                                                              <  2:189158/150‑1 (MQ=255)
       tGTAGGCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCgg                                                                                                                                         <  2:577815/150‑1 (MQ=255)
           ggCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGACGCCTCggtgc                                                                                                                                      <  1:76961/149‑2 (MQ=255)
           ggCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGACGCCTCggtgc                                                                                                                                      >  2:76961/1‑148 (MQ=255)
                      cGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACGGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGa                                                                                                                          <  1:577201/150‑1 (MQ=255)
                          gTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAg                                                                                                                      <  1:597246/150‑1 (MQ=255)
                                         ctgcGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCAc                                                                                                       <  2:417341/150‑1 (MQ=255)
                                              gCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCt                                                                                                  >  2:33015/1‑150 (MQ=255)
                                                aTTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCg                                                                                                <  1:18985/150‑1 (MQ=255)
                                                 tttGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCgg                                                                                               >  1:177151/1‑150 (MQ=255)
                                                        gTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCa                                                                                        <  1:308321/150‑1 (MQ=255)
                                                          tCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCagt                                                                                      >  2:265512/1‑150 (MQ=255)
                                                                    ttttAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCttgtt                                                                            <  1:225972/150‑1 (MQ=255)
                                                                      ttAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCGTGATGTGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGc                                                                          <  1:319747/150‑1 (MQ=255)
                                                                            gCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCATGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGa                                                                    >  2:574176/1‑150 (MQ=255)
                                                                                   gTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGACGCGTCGGTGGTGATGGCAAAGATCAGGGTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGAt                                                             <  2:296895/150‑1 (MQ=255)
                                                                                                                     tCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAAc                           >  2:45476/1‑150 (MQ=255)
                                                                                                                                 aGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGcgg               >  1:160538/1‑148 (MQ=255)
                                                                                                                                  gCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGa              <  1:265512/150‑1 (MQ=255)
                                                                                                                                       tGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGATCGGc         <  2:177151/150‑1 (MQ=255)
                                                                                                                                       tGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGATCGGc         <  1:427264/150‑1 (MQ=255)
                                                                                                                                             accACCGACGCCTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGATCGGCTTCTgg   >  1:477496/3‑150 (MQ=255)
                                                                                                                                              gcACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTAACATACTGTGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGATCGGCTTCTGGt  >  1:601699/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
ATCACGTTGTAGGCGATATAGCCGAGGTTGAAACCGGCCTGCTGCGGCATTTGCAGGTTCTTGTCGGCTTTTAGCGGCTCGATGTCGGCCGGGCGTGGGAACAGGGTGACCTGGCACTCGTTCTTCTTGAGCTTCTGCATGCGCACCGAGGCGTCGGTGGTGATGGCAAAGATCAGGTTGTCGATCTTCACGTCCTCGGGCTGCCAGTAGTCCTTGTTGCCCTTGAAGCGGATCTGCGCGTCCTTCTGGTATTTGCTGAACACGAACGGGCCGGTGCCGATCGGCTTCTGGT  >  NC_002947/1027850‑1028141

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: