Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 664,812 C→T 18.1% G127R (GGG→AGG)  rlpA ← septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913664,8120CT18.1% 50.4 / 4.1 22G127R (GGG→AGG) rlpAseptal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (2/2);  total (12/10)
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

CAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCC  >  NC_000913/664680‑664951
                                                                                                                                    |                                                                                                                                           
caaTAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGa                                                                                                                            <  2:58204/150‑1 (MQ=255)
                aCGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTggg                                                                                                            >  2:107922/1‑150 (MQ=255)
                       ttGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGa                                                                                                     >  2:115428/1‑150 (MQ=255)
                             ttgttTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGtt                                                                                               <  2:194126/150‑1 (MQ=255)
                                tttGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGa                                                                                            >  1:193751/1‑150 (MQ=255)
                                tttGACGTGTTAAGACGGTCAGCTGCCGAGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATACTGGCGTAGCTGGGGATCGGACGCGTTGGa                                                                                            >  2:188256/1‑150 (MQ=255)
                                      gtgtTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCt                                                                                      >  1:190588/1‑150 (MQ=255)
                                      gtgtTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCt                                                                                      >  1:80437/1‑150 (MQ=255)
                                      gtgtTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTAGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGACGCGTTGGATGGGCt                                                                                      >  1:190579/1‑150 (MQ=255)
                                       tgtTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTg                                                                                     <  2:248338/150‑1 (MQ=255)
                                                  tCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCTGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCtgt                                                                          >  2:51472/1‑150 (MQ=255)
                                                  tCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCTGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCtgt                                                                          >  2:147704/1‑150 (MQ=255)
                                                                         aaTAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCaga                                                   <  1:146837/150‑1 (MQ=255)
                                                                         aaTAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCaga                                                   <  1:146840/150‑1 (MQ=255)
                                                                           tAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAgg                                                 <  1:135229/150‑1 (MQ=255)
                                                                           tAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAgg                                                 <  1:135221/150‑1 (MQ=255)
                                                                                                  ccGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCg                          >  1:222049/1‑150 (MQ=255)
                                                                                                     cgATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGca                       <  2:222049/150‑1 (MQ=255)
                                                                                                          aTTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATa                  >  1:48576/1‑150 (MQ=255)
                                                                                                                      cacGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAg      >  2:67128/1‑150 (MQ=255)
                                                                                                                          aTCATTCGCCTGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTccc  <  1:51472/150‑1 (MQ=255)
                                                                                                                          aTCATTCGCCTGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTccc  <  1:147704/150‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                           
CAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTGATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCC  >  NC_000913/664680‑664951

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