Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 709,068 G→A G245G (GGG→GGA chiP → chitobiose outer membrane channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913709,0680GA93.3% 41.5 / ‑5.9 15G245G (GGG→GGAchiPchitobiose outer membrane channel
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (6/8);  minor base T (1/0);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGGTATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATG  >  NC_000913/708940‑709195
                                                                                                                                |                                                                                                                               
cATCTGGAGATGGATAAGTTTTCTCAGAACCATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGAtata                                                                                                                             <  1:196642/133‑1 (MQ=255)
     ggAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGAt                                                                                                                        <  1:404744/133‑1 (MQ=255)
                    ttATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGCCGGAAGGTTATATCGAACAATGTTTTGCCaaa                                                                                                         >  2:194601/1‑133 (MQ=255)
                           aaCGATAAAACCACCAAAGTTGATTATCTGCACCCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCt                                                                                                  >  2:8819/1‑133 (MQ=255)
                                        ccaAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGatat                                                                                     >  2:390377/1‑133 (MQ=255)
                                                          aCTCCTTTGGGGCGAAATACGCCTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAac                                                                   <  2:398039/133‑1 (MQ=255)
                                                          aCTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAac                                                                   <  1:390377/133‑1 (MQ=255)
                                                          aCTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAac                                                                   <  1:8819/133‑1 (MQ=255)
                                                               tttGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAAccacca                                                              <  1:492209/133‑1 (MQ=255)
                                                                  gaggCGAAATCCGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCACTATTCTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGca                                                           <  2:432142/131‑2 (MQ=255)
                                                                    ggCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTAc                                                         <  1:81899/133‑1 (MQ=255)
                                                                        aaaTACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGt                                                     >  1:96341/1‑133 (MQ=255)
                                                                                   cAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTAcc                                          >  2:393541/1‑133 (MQ=255)
                                                                                                            gCATTTGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACCGTACCCGAGATAAAGTTGACGTTCGCAGCg                 >  1:197544/1‑133 (MQ=255)
                                                                                                                 tGGTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAAc            >  2:5928/1‑133 (MQ=255)
                                                                                                                  ggTCAGGCGGAAGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACg           <  2:96341/133‑1 (MQ=255)
                                                                                                                           taaGGATATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCGCCAGCTACCAGTTCTACGGTACCCGAGATAAAGGTGACGATCGCAGCGTCAACGACCTTTATg  <  2:453263/132‑1 (MQ=255)
                                                                                                                                |                                                                                                                               
CATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGACTTCAAAAATAACTTCGTACTGGAAGCGGCATTTGGTCAGGCGGAAGGGTATATCGATCAATATTTTGCCAAAGCCAGCTACAAATTTGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATG  >  NC_000913/708940‑709195

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