Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,909,227 A→C 33.0% 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,2270AC33.0% 32.7 / 4.0 34G146G (GGT→GGGpyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/12);  new base C (5/7);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGATTGCGTTAGTGATGTGCGGAATAACCTGCACGGTTGCGCCGAGGTAGTCACCGCGGCGTTCTTTACGCAGAACGTCAGAGTAGATACGACCCGTGGTGAAGTTGTTGCGGCGGCTCATTTTGGTACGAATGAAACGCTCGTAGTGCCCCAGGTCCAGG  >  NC_000913/2909000‑2909455
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
gCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCa                                                                                                                                                                                                                                      >  1:548225/1‑228 (MQ=255)
 caggtgcgttcgCCGGAACAGCGGGATGTGAACGAAAAAGCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGGGTGGTGAGTCGGTTTGGTTTTGACTGACCCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCAcc                                                                                                                                                                                                                                    <  2:273681/223‑1 (MQ=255)
  cggttggtcgggcGGAAAAGCGCGTCCTGAACTATAATGCAGGATATAAGGCTGGATACCGAAGGAGAGCAGCTCTTTTACAGAGTGGTAACTCGGTTTGGTTTTCACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGTATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCAcc                                                                                                                                                                                                                                    <  2:431964/216‑1 (MQ=255)
           cggcggaatgggccgATGTGAACGAAAAGTCAGGATGGGCGGCTGGATGCCGATGGGGAGAAGCTCTTGTACCGAGTGCTGAGTCGGTTTGGTTTTGACCTCACCAGACGCTGCCATGTACGGCCCCGGCGTCAGGTGCATAAACAGTGTGGTCTCAGGCCCAATTTCAAAAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTCccg                                                                                                                                                                                                                          <  2:750512/216‑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)
                                                             tttgtaaaggggacagctgtgaCAGAGTGCTGGATCGGTTTGGTTTTGACTTCCCCAGACGCTGCCATGTACGGCACCAGCGGAAGGTTCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                          <  1:314386/207‑1 (MQ=255)
                                                                      agcagcTCTTTTACAGAGTTCTGGTTCGGTTTGTTTTGGACTTCACCAGACGCTGCCATTTACGCAACCAGCGTCAGTTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAAGCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCACGCGCTCTTTGa                                                                                                                                                                <  1:581675/229‑1 (MQ=255)
                                                                                      agTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGGGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCACCACGCGCTCTTTGATTGCGCTAGTGCtgt                                                                                                                                                 >  1:640213/1‑228 (MQ=255)
                                                                                        tGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTAACGCCGATTTCTACCCGTACTACGGCATGACCTTCGCCCCCTTCCCGCACCCGCTCTTTGtttcgactcttataaaca                                                                                                                                               >  2:314386/1‑213 (MQ=255)
                                                                                        tGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTCAGTCATGACCTTCCCCCCCTTCCAGCCCCGCCTCTTTCATTGCGTTAGGAAAGTGCg                                                                                                                                              >  2:425672/1‑229 (MQ=255)
                                                                                                  gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCCCCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                          <  1:364850/191‑1 (MQ=255)
                                                                                                  gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGATATCACCTACTGTACCGCCGATTTCTACCAGTACTACGTCATGACCTTCGCCACCTTCCAGCAcg                                                                                                                                                                          >  2:364850/1‑191 (MQ=255)
                                                                                                  gtttggtttTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTTCGAGGAACGGCAAGGATTCGAAATCCCCTACTGGACCGCCGCTTTCTACCAGTACTAAGTCATGACCTTCCCCACCTTCCAGCCGGCGCTCTTTGATTGCGTTAgagtgtttcgcaaaaacctc                                                                                                                                    >  2:198966/1‑210 (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/2909000‑2909455

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