Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,909,229 A→C 32.3% G146G (GGT→GGG pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,909,2290AC32.3% 28.7 / 2.9 29G146G (GGT→GGGpyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base C (3/7);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGG  >  NC_000913/2909002‑2909457
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
gCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCa                                                                                                                                                                                                                                      >  1:548225/1‑228 (MQ=255)
 caggtgcgttcgCCGGAACAGCGGGATGTGAACGAAAAAGCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGGGTGGTGAGTCGGTTTGGTTTTGACTGACCCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCAcc                                                                                                                                                                                                                                    <  2:273681/223‑1 (MQ=255)
                gAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGATTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAACCAGAGTGGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCCCTTCGAGGCAAGGCAAGGATTCGAAATCCCCTACTGTActcgcatt                                                                                                                                                                                                                     >  2:249159/1‑222 (MQ=255)
                 accAGCGCGATCTGAACGACAAACCAAGATGTCAGGCTTGATCACTATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTGAACTTCCCCAGCCGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGGGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGAAAGGATTCGATATCCCCTACTGTACCGCCGAtt                                                                                                                                                                                                                    <  1:683323/227‑1 (MQ=255)
                      tgctgtatGAACGATAAATCGGGATGTGAGGCTGGAGACCGATGGAGACCAGCTCTTTTACAGAGTGCTGAGTCGGTGTGGTTTTGACCTCCCCAGACGCTGCCATGGACGGCCCTCGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTCGCAAATCTCTTCAAGGAACGGCAAGGATTCTATATCCCCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                               <  2:250850/222‑1 (MQ=255)
                      gtggatcgaaaCGACAAAGCAGGGTGTCAGGCTGGAACCCGATGGAGGGCAGCCTTTTTACAGAGTTCTGAGTGGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGTGGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTCCTGTACCGCCGATTTCTAc                                                                                                                                                                                                               <  1:308604/220‑1 (MQ=255)
                      ggttgtCTTAACGACAAATCAGAGGTACGGGTTGGATACCGTTGGGGGGCAGCCTTTTTACAGAGTGCTGAGGCGGTTTGGTTTTGCTTACACCAGCCGTGGCCATTGACGGCACCAGCGTCAGGTGCTAAAACAGGGTGTGCTCACGGCCAATTTCAACAGCCATCGGGCGAATCGCTTCGAGGAACGGCAAGAATTCGATATCACCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                               <  2:638557/224‑1 (MQ=255)
                      ggcgATCTGAACGACAAATCCGGATGGCAGGCTGGATTCCGATGGAGAGCAGCTCTTTTACAGAGTGCTGGGGCGGTTTGGTTTGGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                               <  1:733796/228‑1 (MQ=255)
                      cgcgATCGGAACGAAAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTAc                                                                                                                                                                                                               <  1:564212/229‑1 (MQ=255)
                              gAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCCGGTGCATAAACAGAGTGGGCTCACGGCCAATTTCAAAAGCCATCTGGCGAATCGCCTCGGGGAACGGCCAGGATTCGAATTACCCCACTTTACCGCCGCtttttccccatac                                                                                                                                                                                                         >  2:184990/1‑217 (MQ=255)
                                    cAAATCAGGATTGCAGGCTGGTACCCGATGGAGAGCCGCTTTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCCCCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCa                                                                                                                                                                                                  <  1:229881/228‑1 (MQ=255)
                                                            aCCGATGGAGAGCAGTCTTTTTCCAGATTGCTTAGTCGGTTTGGTTTTGCCTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACGGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCTCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                          <  1:2303/229‑1 (MQ=255)
                                                                      agcagcTCTTTTACAGAGTTCTGGTTCGGTTTGTTTTGGACTTCACCAGACGCTGCCATTTACGCAACCAGCGTCAGTTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAAGCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGa                                                                                                                                                                <  1:581675/229‑1 (MQ=255)
                                                                                      agTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCACCACGCGCTCTTTGATTGCGCTAGTGCtgt                                                                                                                                                 >  1:640213/1‑228 (MQ=255)
                                                                                        tGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTCAGTCATGACCTTCCCCCCCTTCCAGCCCCGCCTCTTTCATTGCGTTAGGAAAGTGCg                                                                                                                                              >  2:425672/1‑229 (MQ=255)
                                                                                                  gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                          <  1:364850/191‑1 (MQ=255)
                                                                                                  gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                          >  2:364850/1‑191 (MQ=255)
                                                                                                       gtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCAATCTGGCGAATCGCTTCGAGGAACGGCAAGGCTTCGATATCACCTACTGTACCGCCGATTTCTACCAGTCCTACGTCATGAACTTCCCCCCCTTCCCGCACGCGCTCTTTTATTTCGTTTGTCATGTTCGGAAACAACTGCaacg                                                                                                                               >  2:487911/1‑226 (MQ=255)
                                                                                                                                       cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTCCTACGCCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTCCGTTATTGATGTGCGGAATAACCTGCACGGTTCGCCCCAGGCAGTCACCGCGGCGTTCttt                                                                                                >  2:231018/1‑228 (MQ=255)
                                                                                                                                       cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTCGGGATGTGCGGAATAACCTGCACGGTTGCGCCGTGGTAGTCACCGCGGCGGTCgtt                                                                                                >  1:588823/1‑228 (MQ=255)
                                                                                                                                       cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCttt                                                                                                >  1:161828/1‑228 (MQ=255)
                                                                                                                                       cACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACAACGTCATGACCTTCGCCACCTTCCACCCCGCGCTCTTTGATTCCGTTAGGGACTGGCGGACTACCCTCCACGCTGCCGCCCTGGTAGTTCACGCGGCGGTAttt                                                                                                >  1:354108/1‑228 (MQ=255)
                                                                                                                                            gcgcCAGGGCCATAAACAGAGTGGGCTCCCGGCCAATTTCAACAGCCATCGGGCGAATCGCTTCGAGAAACGGCAAGGATTCGATATCACCTACTGTCCCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAg                                                                                          <  1:700517/229‑1 (MQ=255)
                                                                                                                                            gCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGGGGCGTTCTTTACGCAg                                                                                          >  1:577277/1‑229 (MQ=255)
                                                                                                                                                 gggTGGAAAAAAAGAGTGTGCCACCCGCCCATTTCACCCTCCTTGTGGTGAATGCTTTAGGTGAACGGCCTGGGTTCGTTCTCACCTACGTTACCGCCGATTTCTACCAGTTCTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTATTTACGCAGAACGt                                                                                     <  2:640212/228‑1 (MQ=255)
                                                                                                                                                  ggTGCAAAAAAAGAGTGTGCCCTCGGCCAATTCCAACAGCCATCTGGCGAACCGCTTCGAGGAACGGCAAGGATTCGATATCACCTCCTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTc                                                                                    <  1:249159/229‑1 (MQ=255)
                                                                                                                                                   gTGCATAAACAGATTGTGCCCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTCCTGTTCCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTc                                                                                    <  1:585380/228‑1 (MQ=255)
                                                                                                                                                              gagTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGCATTCGATATCACCGACTGTAACGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGCGATGTGCGGACTAACCTGCACGGTTGGGCCGCGGTAGACACCGCCGCGTTCGTTACGCAAACAGTCAGGGTAGATAca                                                                        >  2:636658/1‑228 (MQ=255)
                                                                                                                                                               agTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCCCGCGCTCTTTGATTGCGTTAGTGATGTGGGGAATAACCTGCACGGTTGCGCCGGGGTAGTCACCGCGGCGTTCTTTACGCAAAACGCCAGGGTAGATAcc                                                                        >  2:585379/1‑227 (MQ=255)
                                                                                                                                                                                                           cGAGGAACGGCAAGGATTCGATATCCCCTACTGTTCCGCCGATTTCTGACAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGa                            <  1:363388/228‑1 (MQ=255)
                                                                                                                                                                                                                                    aCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCTGGGTGAAGTCGTTGCGGCGGCGCCATTTGGTCAGAAGGGAACGCTCGTGGTGCCcccgggaccgc  >  1:689371/1‑223 (MQ=255)
                                                                                                                                                                                                                                    aCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCAATTTGGGACCAAAGAAACGCGCGAAGTGGCccatgtcccg   >  1:444381/1‑226 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
GCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGG  >  NC_000913/2909002‑2909457

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