Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,132,627 (T)6→7 coding (12/891 nt) metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,132,6211.T97.1% 116.0 / NA 35S2S (AGC→AGT) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (23/11);  minor base A (0/1);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGC‑TTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCA  >  NC_000913/4132493‑4132748
                                                                                                                                 |                                                                                                                               
aCCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGCCGTCTAAACGGATAGATATGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGGCTGAGC‑tttt                                                                                                                             >  2:255770/1‑133 (MQ=255)
                          cGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCt                                                                                                    <  2:404231/133‑1 (MQ=255)
                               ttCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATc                                                                                               >  1:312763/1‑133 (MQ=255)
                                              aTCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAg                                                                                >  1:429868/1‑133 (MQ=255)
                                                   tACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAg                                                                           >  2:259675/1‑133 (MQ=255)
                                                      aTCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAgggg                                                                        <  1:259675/133‑1 (MQ=255)
                                                                       gATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTttcgtt                                                       <  1:112545/133‑1 (MQ=255)
                                                                          aGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGa                                                    >  2:477698/1‑133 (MQ=255)
                                                                                 cacaacacaaCATATAACTACAAGCGATTGATGAGGAAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTc                                             >  2:425516/1‑133 (MQ=255)
                                                                                      cacaacaTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGACGTCCAGGGGCAGGTTACCGTTTCGTACGAGATTTTCccgcc                                        >  2:134509/1‑133 (MQ=255)
                                                                                         aacaTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCcgcg                                     >  2:422187/1‑133 (MQ=255)
                                                                                           caTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCcgcg                                     <  2:57194/131‑1 (MQ=255)
                                                                                           caTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCcgcg                                     >  1:57194/1‑131 (MQ=255)
                                                                                           caTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTa                                   >  2:60625/1‑133 (MQ=255)
                                                                                            atatAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCGCGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGGGTTTTTCCCGCCGCGTAc                                  <  2:113033/133‑1 (MQ=255)
                                                                                              ataACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCa                                >  1:117648/1‑133 (MQ=255)
                                                                                                aaCTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTGTCGTTCGAGTTATTCCCGCCGCGTACCAGt                              >  1:172802/1‑133 (MQ=255)
                                                                                                     cAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGaaa                          <  1:45693/132‑1 (MQ=255)
                                                                                                     cAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGaaa                          >  2:45693/1‑132 (MQ=255)
                                                                                                     cAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAAt                         >  1:283537/1‑133 (MQ=255)
                                                                                                     cAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAAt                         >  1:56431/1‑133 (MQ=255)
                                                                                                      aaGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAAt                         >  2:314095/1‑132 (MQ=255)
                                                                                                      aaGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAAt                         <  1:314095/132‑1 (MQ=255)
                                                                                                      aaGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATg                        >  2:452236/1‑133 (MQ=255)
                                                                                                      aaGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATg                        >  1:73842/1‑133 (MQ=255)
                                                                                                      aaGCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGATCGGGTTATTCCCGCCCCGTAACAGTGAAATg                        >  2:159186/1‑133 (MQ=255)
                                                                                                        gCGATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGa                      >  1:355645/1‑133 (MQ=255)
                                                                                                          gATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGTAGc                    <  2:56431/133‑1 (MQ=255)
                                                                                                          gATTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGc                    <  1:422187/133‑1 (MQ=255)
                                                                                                           aTTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGTCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCa                   >  1:383125/1‑133 (MQ=255)
                                                                                                           aTTGATGAGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCa                   >  2:493835/1‑133 (MQ=255)
                                                                                                             tgatgaGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGa                 <  1:402158/133‑1 (MQ=255)
                                                                                                               atgaGGTAAGGTATGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGAcc               <  1:477698/133‑1 (MQ=255)
                                                                                                                  aggtaaggcaTGAGCATTCTTTCACGCCAGTCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCtg            <  1:134509/133‑1 (MQ=255)
                                                                                                                            tGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCa  >  2:265260/1‑133 (MQ=255)
                                                                                                                            tGAGCTTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCa  >  2:286374/1‑133 (MQ=255)
                                                                                                                                 |                                                                                                                               
ACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGC‑TTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCA  >  NC_000913/4132493‑4132748

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