Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,133,252 C→T P213S (CCG→TCG)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,133,2520CT96.4% 70.3 / ‑4.3 28P213S (CCG→TCG) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (5/22);  total (5/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

GTATGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGAC  >  NC_000913/4133005‑4133495
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gatttgttCTGACCTGGTGACGCTGGTAAAAGAAGTGGCAGATTTCGATATCGCAGTGGCGGCGTATCCTGAAGTTCACACGGAAGCAAAAAGCGCTCAGGCGGATGTGCTTAATCTGAAACGCAAAGTTGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCgg                                                                                                                                                                                                                                                   <  2:229632/245‑1 (MQ=255)
          gACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTaa                                                                                                                                                                                                                                            <  1:198951/248‑1 (MQ=255)
                  gACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAgg                                                                                                                                                                                                                                 <  1:493220/251‑1 (MQ=255)
                                 agTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCg                                                                                                                                                                                                                  <  2:257668/251‑1 (MQ=255)
                                  gTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGGTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGa                                                                                                                                                                                                                 <  2:528330/251‑1 (MQ=255)
                                                             gcggagcctgAAGTTCACCCGGAAGCAAAAAGCGCTAAGGCGGAGTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACGCAGTTCTTCTTCGATGGCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTg                                                                                                                                                                                     <  2:216939/242‑1 (MQ=255)
                                                             gcgTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGGAATTTTTCCTGGAAATTTGTCGGGGTCTAACTTTTAAAAGGCGAAGAAATTTTTCGGTATGACCCAAC‑TGCGTATTCCGGCGTg                                                                                                                                                                                     >  1:390889/1‑251 (MQ=255)
                                                                        aGTTCCCCCGTAAGAAAAAAGCGTACAGGCGGGTTTGCTGAATCTGAAGCGTAAAGGGGATGCCGGAGCCTACCGCGCGGTTACTCAGGTCTTCTGCAATGTCGAAAGCTACCTGCGTTGTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCTATATGACC‑AACGTGCGTAATCCGGCGTGGATGGCGCaaa                                                                                                                                                                          <  2:431775/251‑1 (MQ=255)
                                                                                cGGAAGCAAAAAGAGCTAAGGCGGATTTGCTTAATCTTAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTTCTCAGTTCTTCTTCGATTTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAAGTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGTCGTGGATGGCGCAAATGTTCGAc                                                                                                                                                                  <  2:366234/251‑1 (MQ=255)
                                                                                      cAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCGTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTg                                                                                                                                                            <  1:435938/251‑1 (MQ=255)
                                                                                              gcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGgatgatga                                                                                                                                                    <  1:461822/251‑1 (MQ=255)
                                                                                                 cAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCGGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGc                                                                                                                                                 <  1:322535/251‑1 (MQ=255)
                                                                                                            gCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCa                                                                                                                                      <  2:390889/251‑1 (MQ=255)
                                                                                                                     gcaactcAAAGTTGATGCCGGAGCCAACCGCGCGCGTACGCAGGTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGTGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCTATATGCCC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGTAAACta                                                                                                                                 <  2:72751/241‑2 (MQ=255)
                                                                                                                           cAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCaa                                                                                                                                     <  1:138136/237‑1 (MQ=255)
                                                                                                                                                                                      aCCTGCGTTTTCGTGACCGCGGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACGTTAAACAGGCGAAGAAATTTGCTGATATGACC‑AACGTGCGGATTCCGGCGTGGATGGCGCAAATTTTCGACGGTCTGGAGGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACt                                                              <  2:411113/249‑1 (MQ=255)
                                                                                                                                                                                          gCGTTTTCGTGACCGCTGTGTATCTGCGGGCATTGAGGTGGAAATTATTCCGGGAAGTTTGTCTGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGTCGTTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTTCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGGGAAAGAGTTCCACTTCtata                                                        <  2:403293/251‑1 (MQ=255)
                                                                                                                                                                                          gCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTTGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCtat                                                         <  2:72210/250‑1 (MQ=255)
                                                                                                                                                                                              tttCGTGACCGCTGTGTATCGGCGGGCATTTGTGTGGGAATTATTTCGGGGATTTTGTCGGGATCTAACTTTAAACAGGCGGAGAAAATTTCCGGTATGACC‑AACGTGCGGTTTCCGGGGGGGGTGGCGCAAATGTTTGACGGTCTGGGTGGTGATGCCGGAACCCCCAAACTGGGTTGCGGGGATATTTCCATGGGTATGGTGAAGATTTTTAGCCGTGGAGGGGTGGAAGATTTCCACTTCtataacct                                                    >  1:44127/1‑247 (MQ=255)
                                                                                                                                                                                                                   gcggGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTa                               <  1:581565/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggtgGAAATTATTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGcc                    <  1:273368/250‑1 (MQ=255)
                                                                                                                                                                                                                                tgGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCAt                  >  1:495228/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaaTTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACg               <  1:93939/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaTTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGCATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAAGCGATTTGCAATACGc              >  2:659931/1‑251 (MQ=255)
                                                                                                                                                                                                                                     attattCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGata                                                                                                        >  2:380707/1‑160 (MQ=255)
                                                                                                                                                                                                                                     attattCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGata                                                                                                        <  1:380707/160‑1 (MQ=255)
                                                                                                                                                                                                                                              ggAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCg    <  1:470871/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aaTTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGAc  <  1:325267/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GTATGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACC‑AACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGAC  >  NC_000913/4133005‑4133495

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