Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 1,256,820 A→G I171T (ATC→ACC)  prs ← phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,256,8200AG92.5% 121.9 / 0.2 40I171T (ATC→ACC) prsphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (20/17);  total (22/18)
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

AGTCACGACCTGCAACGTCACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCC  >  CP009273/1256686‑1256948
                                                                                                                                      |                                                                                                                                
aGTCACGACCTGCAACGTCACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTc                                                                                                                              >  1:2689599/1‑139 (MQ=255)
 gTCACGACCTGCAACGTCACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTcc                                                                                                                             >  1:690728/1‑139 (MQ=255)
    aCGACCTGCAACGTCACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGga                                                                                                                          <  1:530431/139‑1 (MQ=255)
            cAACGTCACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACcac                                                                                                                  >  1:55421/1‑139 (MQ=255)
                  cACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTgg                                                                                                            >  1:1062127/1‑139 (MQ=255)
                    ccGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAACTGGGt                                                                                                          >  1:2297851/1‑139 (MQ=255)
                            atGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGa                                                                                                  <  1:1985604/139‑1 (MQ=255)
                              gCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTACCATATCGGTATAGTTCAGAAGCTTAGCGATAGCGCGGGCACGCACAATGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGAt                                                                                                 >  1:1194538/1‑138 (MQ=255)
                                 tCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAg                                                                                             <  1:1883052/139‑1 (MQ=255)
                                  cACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCa                                                                                              >  1:801173/1‑137 (MQ=255)
                                  cACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCa                                                                                              >  1:561449/1‑137 (MQ=255)
                                    ccTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCt                                                                                           >  1:1430038/1‑138 (MQ=255)
                                       gtgAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATc                                                                                                      <  1:283631/124‑1 (MQ=255)
                                        tgAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAAc                                                                                                                     <  1:1419677/108‑1 (MQ=255)
                                         gAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTgcagca                                                                                     <  1:1266536/139‑1 (MQ=255)
                                            aCGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGt                                                                                  <  1:5316/139‑1 (MQ=255)
                                                 cGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGc                                                                                            >  1:1708005/1‑124 (MQ=255)
                                                  gCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTggg                                                                                                           >  1:2842533/1‑108 (MQ=255)
                                                              cgTTTGTCGATGATTGCTATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCAGTCACGCACATCACCGCCGATGTCCGGAGAATCCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCggc                                                                 <  1:2617716/138‑2 (MQ=255)
                                                                              ccATATCGGTATCGTTCATCAGCTGAGCGATAGCGCGGGCACGCACAACGCCGCCGATGGCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATc                                                >  1:2626771/1‑139 (MQ=255)
                                                                                atatCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCa                                               >  1:1744278/1‑138 (MQ=255)
                                                                                  atCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCggg                                                                 <  1:2089721/118‑1 (MQ=255)
                                                                                         tCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACg                                      >  1:1226448/1‑138 (MQ=255)
                                                                                            ttCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAAc                                                                                                                     <  1:368206/56‑1 (MQ=255)
                                                                                            ttCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATa                                                       >  1:665692/1‑118 (MQ=255)
                                                                                                    cTTAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCCGGATCGGGCTACCAAATACGTTATCAACCGGACCGTCGCAGAAAc                            >  1:2353392/1‑137 (MQ=255)
                                                                                                     ttAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCTGGc                                                                <  1:1395882/100‑1 (MQ=255)
                                                                                                      tAGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTggg                                                                                                           >  1:1792476/1‑56 (MQ=255)
                                                                                                       aGCGATAGCGCGGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTggg                                                                                                           <  1:413807/55‑1 (MQ=255)
                                                                                                               cgcgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATAc                                                      >  1:2818768/1‑100 (MQ=255)
                                                                                                                gcgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGtt                <  1:7000/137‑1 (MQ=255)
                                                                                                                gcgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCa              <  1:319939/139‑1 (MQ=255)
                                                                                                                gcgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCa              <  1:2649577/139‑1 (MQ=255)
                                                                                                                gcgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCa              <  1:2364764/139‑1 (MQ=255)
                                                                                                                 cgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGAt                                                                                                 >  1:1838087/1‑55 (MQ=255)
                                                                                                                 cgGGCACGCACAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCa              >  1:672336/1‑138 (MQ=255)
                                                                                                                    gcacgcacAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCg           <  1:929774/138‑1 (MQ=255)
                                                                                                                        gcacAACGCCGCCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCACATACGTTATCAACCgg                                          <  1:2633594/103‑1 (MQ=255)
                                                                                                                              cgccgcCGGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATcc  >  1:1864355/1‑137 (MQ=255)
                                                                                                                                  gccgGTGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCgaag                                >  1:1210441/1‑103 (MQ=255)
                                                                                                                                      |                                                                                                                                
AGTCACGACCTGCAACGTCACCGATGATATGCATCACCTGTGAAACGTTCGCACGCGGACGACGTTTGTCGATGATTGCCATATCGGTATCGTTCAGCAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCC  >  CP009273/1256686‑1256948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: