Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,377 T→C R161R (CGT→CGC)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,3770TC92.9% 68.0 / ‑2.8 29R161R (CGT→CGC) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (18/8);  total (20/9)
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.78e-01

CGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACAT  >  NC_000913/4625150‑4625622
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
cGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTg                                                                                                                                                                                                                                 >  1:108742/1‑250 (MQ=255)
            cTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGAAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGGGGCGACGCTGAAAGGCGCTGACGCCACTATTCt                                                                                                                                                                                                                    >  1:276405/1‑251 (MQ=255)
                   agaaggagcactgggATGCCGCCCGGTTGGGG‑AAACCCCCCCCCCTGCGCACGC‑TGGTTTGTCTGGTGATGGACATGAATTCCCCCGCCATCCAGATTGAATGTATGGATAAAATTGCCAAACTGGC‑CGGAACGGGCGAGTGG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAg                                                                                                                                                                                                             <  2:69716/238‑1 (MQ=255)
                             acTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGaa                                                                                                                                                                                                    >  1:82004/2‑250 (MQ=255)
                                gtaTGTCGCCCCG‑TGGGGTAAAATCCCGCTCCGGCGCCCGCCGGGTTTG‑CTGGTGATGGGCATGGATTCCCCCGCCATCCAGAGTGATTGTATTGATGAAATTGCCAAACGGGC‑CGGAACGGGCGAGAGG‑GTGGCGGAAGTAACCG‑AACGGCCGATGCGCGGTGACCTCGATTTTCCCTCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGGGAAAATc                                                                                                                                                                                                <  1:242192/249‑1 (MQ=255)
                                  aTGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTg                                                                                                                                                                                              >  1:281338/1‑251 (MQ=255)
                                               ggaggatttcccctactcTGCGCACGTCTGGTTTG‑CTGGTGAGTGACATGATATCCGTCGCCATCCAGATTGAATCTTTGGATAAAATTGCCAAATGGGC‑CGGACCGGGCGTGTTG‑GTGGCGGAAGTACCCG‑AACTGGCGTTGCGCTGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCa                                                                                                                                                                                  <  2:10222/234‑1 (MQ=255)
                                                ggg‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCA‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGGGGCGACGCTGAAATGCGCTGACGCCAATATGCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGGCAgg                                                                                                                                                                                >  2:287385/1‑251 (MQ=255)
                                                           cGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCa                                                                                                                                                                       >  1:110475/1‑250 (MQ=255)
                                                               ccTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTgg                                                                                                                                                                  <  1:331062/251‑1 (MQ=255)
                                                                                      gTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctggg                                                                                                                                            >  1:263125/1‑251 (MQ=255)
                                                                                       tGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAACATCTGCCGCTGCTGCCCGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggc                                                                                                                                           >  2:198825/1‑251 (MQ=255)
                                                                                             acatGGATTCCACGGCCAGCCAGATTGAATGTATGGATGAAATTGCCAAACTGGC‑CGGAACGGGAGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTGACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGaaa                                                                                                                                     <  1:281777/249‑1 (MQ=255)
                                                                                               atGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGCAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCCGCCTGCGCAGCCGCGTGGCGTCGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGATAGt                                                                                                                                   >  1:18625/1‑251 (MQ=255)
                                                                                                           gCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCAACGCTGCAAGGCGCTGACGACAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTATAGCAACTGGTGCTCAAACTGGCAACGCTGGGATGGAAAgggggggttgactc                                                                                                                       >  2:12867/1‑238 (MQ=255)
                                                                                                            ccctcaaGATTGAATGTATTGATGAAATTTCCAA‑CTCGCACGGAACGGGCGAG‑TGTGCTGCGTACGTA‑CCGTATCGGGCGATGCGCGGCGAACCTTATTTTCCCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTTGTGCTCAAGCTGGAAACGCTTGGCTGGAAAGTTGCGCTTGCCt                                                                                                                        <  2:177047/243‑1 (MQ=255)
                                                                                                                ccAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcggcg                                                                                                                  >  1:146462/1‑251 (MQ=255)
                                                                                                                                          tcAAACTGGC‑CGGAACGGGCGAAATG‑GAGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCTATTTTACCGGCAGCCTGCGCAGCCGCGTGTGGACGCTGCAAGGCGATGACGCCAATAGTCTGAAACAGGTGCGTGAAAAGCTGCAGCGTCTTCCAGGCTTAACTCAAGTGGAGCTCAAGCTGGTACCGCGGAGATGGAAGGGGGCTGTTGCCGCCGGCGGCTTTATTTTCTTTGCttactatctg                                                                                        >  2:17171/2‑242 (MQ=255)
                                                                                                                                          ccAAACTGGC‑CGGACCGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTg                                                                                        >  1:299465/1‑251 (MQ=255)
                                                                                                                                                   c‑cGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACa                                                                                 >  1:161059/1‑249 (MQ=255)
                                                                                                                                                          cGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGAGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGAGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCACATGGTGCTCAAGCTGGAAACGCTGGGCTGGAGCGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGAGCGACAAGTTgcgc                                                                         >  1:362267/1‑251 (MQ=255)
                                                                                                                                                               gggaTG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCATGCGCAGCCGCGTGGCGACGCTGAAAGGAGCAGACGCCAATATTCTGCAACAAGTGCGTGAAAATATGACGCTGATACCAGGCGTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAACGCGGAGATTGACTCCGGGGGGCTTACTTTCTTTGCTGTATATCTGCGCGACAAGTTGCGCATGAc                                                                    >  2:107381/3‑251 (MQ=255)
                                                                                                                                                                agaTG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGAcc                                                                   >  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)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
CGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGC‑CGGAACGGGCGAGATG‑GTGGCGGAAGTAACCG‑AACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACAT  >  NC_000913/4625150‑4625622

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