Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,546 C→T L218L (CTG→TTG)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,5460CT86.4% 46.0 / ‑1.7 22L218L (CTG→TTG) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (13/6);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAAAGGCGGAAGT  >  NC_000913/4625301‑4625796
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
cGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGAGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGAGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCACATGGTGCTCAAGCTGGAAACGCTGGGCTGGAGCGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGAGCGACAAGTTgcgc                                                                                                                                                                                                                                                       >  1:362267/1‑251 (MQ=255)
     gggaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCATGCGCAGCCGCGTGGCGACGCTGAAAGGAGCAGACGCCAATATTCTGCAACAAGTGCGTGAAAATATGACGCTGATACCAGGCGTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAACGCGGAGATTGACTCCGGGGGGCTTACTTTCTTTGCTGTATATCTGCGCGACAAGTTGCGCATGAc                                                                                                                                                                                                                                                  >  2:107381/3‑251 (MQ=255)
      agaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGAcc                                                                                                                                                                                                                                                 >  1:225988/1‑251 (MQ=255)
                                             aacggtttttcctccccAGCCTGCGCTCCCGGGTGGCAACGGTGAAAGGCGCTGGCGCTAATATTCTGTACCAGTTTCGGGAAAACGTGCGGCGGATCCCAGGCAAACCGCAACCGGGGCTCAAGCTGAGAACGCTGGGCGGGAAGTTGGGGATGGCTCCCGGCGTCTGTACTTTTTTTGCTAAATATCTTGGCGTCAAGTGGCGCCTGACCGCCGTGGTAGCCAATGAACTGTAGAGCATGGACGGTaaa                                                                                                                                                                                                          <  2:299465/236‑1 (MQ=255)
                                               cTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAtt                                                                                                                                                                                                        >  1:383967/1‑251 (MQ=255)
                                                       tACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                                                 >  1:268364/1‑250 (MQ=255)
                                                          cGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCAAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAtg                                                                                                                                                                                             <  1:12867/251‑1 (MQ=255)
                                                                   gcgcAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCgg                                                                                                                                                                                      <  1:64393/249‑1 (MQ=255)
                                                                       aGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGAGACAt                                                                                                                                                                                >  1:369522/1‑251 (MQ=255)
                                                                                       tGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGAAATGTGATCGGCGACATCGTAGACGCGCAGTAc                                                                                                                                                                >  1:355988/1‑251 (MQ=255)
                                                                                       tGAAAGGCGCTGACGCCAATATTCTGCAACAGGAGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAACCGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAAATCATGGACGGTAAATTTACCTGCAATGTGAACGGCGGCATCGTCGACGCGCAgaac                                                                                                                                                                >  1:355991/1‑248 (MQ=255)
                                                                                                  gACGCCAATATACTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAAc                                                                                                                                                     <  1:244177/251‑1 (MQ=255)
                                                                                                                  aaCAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGCAAACTCTGACTCGCCTcgcg                                                                                                                                     >  1:129821/1‑251 (MQ=255)
                                                                                                                     aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCACATGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATAGCCTCAGGCGGCTTTACTTTCTTTGCTGAATCCCTGCGCGACAAGCTGCGCCTGACCGCAGTGGTAGCCAATGAACTGTAGATCATGGACCGTAAATTTACAGGAAATGTGATCGGCGTCATCGTAGACGAGCCGTACAAAGCGACAACTCTGACTCGCCTcgcggtg                                                                                                                                  >  2:364833/1‑248 (MQ=255)
                                                                                                                                        tcTGATACCAGACTTAACACAACTGTTGCTCAAGCTGGAAACGCTGGTCTGAAAAGTGGAGATTGCCTCCGTCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACACGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGAAAATTTACCGCCAATGAGCACGGCCTCACCGTCGACGCGCAGTACAACGCGAATACCCTGACTCGCCTCACGCAGGAGTCTGAAATACCGCttg                                                                                                               >  2:223090/2‑249 (MQ=255)
                                                                                                                                                          gCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattg                                                                                             <  1:198825/251‑1 (MQ=255)
                                                                                                                                                                                                     aTTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGcg                                                  >  1:110396/1‑250 (MQ=255)
                                                                                                                                                                                                                   cTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTAcc                                    >  1:8670/1‑251 (MQ=255)
                                                                                                                                                                                                                      tcctttctttGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCAt                                  <  1:161433/248‑1 (MQ=255)
                                                                                                                                                                                                                          ttctttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGTGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgcca                              <  1:287385/250‑1 (MQ=255)
                                                                                                                                                                                                                              ttGCTTAATCCCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagcc                          <  1:364833/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGtgaatg               >  1:167047/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     cTGCGCCTGACCGCCGTGGTAGCCATTGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAAAGGCGGAAGt  >  1:125708/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
CGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAAAGGCGGAAGT  >  NC_000913/4625301‑4625796

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: