Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,041,643 C→T P111S (CCA→TCA)  ygfZ → folate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,6430CT95.7% 54.4 / ‑5.7 23P111S (CCA→TCA) ygfZfolate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (10/12);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGCTGACCGATAAACTGCGCGG  >  NC_000913/3041400‑3041848
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                             
gCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCCGATGTCAGCAAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCCTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAgacg                                                                                                                                                                                                         >  1:269291/1‑251 (MQ=255)
 ccaccaTTACTGACGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGCTGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAgacgac                                                                                                                                                                                                       >  1:391508/1‑251 (MQ=255)
                                          cAGGTGACAGCAGAGGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgcgcgn                                                                                                                                                              <  1:146628/251‑2 (MQ=255)
                                                   gcagATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCACTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGc                                                                                                                                                     >  1:11122/1‑250 (MQ=255)
                                                    cagATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGttn                                                                                                                                                                                 <  1:101868/222‑2 (MQ=255)
                                                           cAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGGAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTgg                                                                                                                                                                                      <  1:26129/210‑1 (MQ=255)
                                                              ccAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGgtgt                                                                                                                                                                                   >  1:216225/1‑210 (MQ=255)
                                                                                                        cGCCAAAGGTAAAATGTGGAGCTATTGACGTCTGTTCCGCGACGGCGATGGGTTTGCATGGATTGAACGGTGCAGCGGGGGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATGCTCGAAAGTGCTCATCGCGTCAGTCGTCGAGCGGGTGCTGCTTGGGGTTGCCGGTTTGCGGGCGCGCTCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAGGAAGgcg                                                                                                <  1:329646/251‑1 (MQ=255)
                                                                                                        cGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                                                                                                                                                       >  1:500088/1‑180 (MQ=255)
                                                                                                         gCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAgac                                                                                                                                                                                                          <  1:328218/144‑1 (MQ=255)
                                                                                                                                      tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGGGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTgg                                                                                                                                                                                      <  1:428638/135‑1 (MQ=255)
                                                                                                                                             cgcgACTGCGATGGCTTTGCATGGATGGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAGCAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACg                                                           <  1:465431/251‑1 (MQ=255)
                                                                                                                                                      gATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                                                                                                                                                       <  1:6466/134‑1 (MQ=255)
                                                                                                                                                         ggCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCTCCGCGCTGGCAAATCTCTTTa                                                                                                                                          >  1:410531/1‑160 (MQ=255)
                                                                                                                                                                   ggATTGAACGGCGCAGCGTGCGGGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCNCTGNNNNNtc                                                                                                                                                >  2:319068/1‑144 (MQ=255)
                                                                                                                                                                        gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt                                                                                                                                                                           >  1:418843/1‑112 (MQ=255)
                                                                                                                                                                                                 aGNNNNNTGANCTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAAc                                                                                                                  <  2:311907/144‑1 (MQ=255)
                                                                                                                                                                                                     gactgaACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGCTGACCGATAAACTGcgcg   <  1:42381/251‑1 (MQ=255)
                                                                                                                                                                                                      actgaactgaAAAAATATGCGGTATTCTCTAAAGTGACCAGCGGGTCAGACGACGTGCGTGTGCTGCGTGGTGTTGCCGGTTTGCAGGCGCGCGCCGCGCTGGCAAAACTCTTTATCGAACTGCCTTCGAAAGAAAAACAGGGATTCAAAGATGGCGCGACCGCTTTGCTATGGGTGGAACACCCGGCAGAATGTTGCCTGAGCGTAACCGATGGATCGTCTTCTAATATGCTGACCGATAAACTGCGCgg  <  1:423365/251‑1 (MQ=255)
                                                                                                                                                                                                       ctgaactgaaAAAATATGCGGTATTCTCTAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgcgc                                                                                                                                                                >  1:198754/1‑92 (MQ=255)
                                                                                                                                                                                                             ngncanaaTATNCGGTNTTCTCTAAAGTGACCNTCGCGTCAGACGACGAGCGNNTNCTGCTNGgtgt                                                                                                                                                                                   <  2:216225/63‑1 (MQ=255)
                                                                                                                                                                                                                                 ctctAAAGTGACCATCGCGTCAGACGACGAGCGTGTGCTGCTCGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGc                                                                 >  1:79006/1‑162 (MQ=255)
                                                                                                                                                                                                                                              agcgcgTCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCCGTTTTCAGGCGCGTTCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaag                                                                                                                        <  1:76039/92‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                             
GCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGTTTGAACACCCGGCAGAACGTTTCCTGATCGTAACCGATGAAGCTACTGCTAATATGCTGACCGATAAACTGCGCGG  >  NC_000913/3041400‑3041848

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