Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,159,200 A→G 8.9% L357L (CTT→CTC)  U712_16220 ← Purine catabolism protein PucG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,159,2000AG 8.9% 121.1 / 4.4 45L357L (CTT→CTC) U712_16220Purine catabolism protein PucG
Reads supporting (aligned to +/- strand):  ref base A (24/17);  new base G (2/2);  total (26/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 = 4.27e-01

AGCGCTTTTCCTGCTTCGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCAAACAATCT  >  CP006881/3159072‑3159337
                                                                                                                                |                                                                                                                                         
atcgcTTTTCCTGCTTCGATTCCTGAGGTATGTCTGAGAAGCACAGCCTCTAACCCCGCAAGCACGAATAGCAAGTTTTCTTTTCTGCAGCTGTACCCCATCGTTCCGATCGTCCAGATCTTCCCTGCAAGCGGGCCGa                                                                                                                                 <  7:294219/137‑1 (MQ=255)
aGCGCTTTTCCTGCTTCGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCg                                                                                                                                  >  4:181252/1‑138 (MQ=255)
     ttttCCTGCTTCGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAg                                                                                                                            >  8:156183/1‑139 (MQ=255)
     ttttCCTGCTTCGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAg                                                                                                                            >  3:260446/1‑139 (MQ=255)
       ttCCTGCTTCGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAg                                                                                                                            <  4:82806/137‑1 (MQ=255)
                cGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGAttt                                                                                                                  <  8:162246/138‑1 (MQ=255)
                 gATTCCTGCGGTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTGTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCa                                                                                                                <  2:154147/139‑1 (MQ=255)
                                      aaGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCagcagc                                                                                           <  7:156183/139‑1 (MQ=255)
                                        gCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCAt                                                                                         >  4:282306/1‑139 (MQ=255)
                                         cacaGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATg                                                                                        >  5:8043/1‑139 (MQ=255)
                                              ccTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACg                                                                                   >  6:81120/1‑139 (MQ=255)
                                                  tAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGCTCCTCCAGATCTTCCCGGCAAGCGGGCCGAATGAGCTGGCGATTTCCATGCCGAATTCTGCCAGCAGACTGTCACGGa                                                                                 >  7:301820/1‑137 (MQ=255)
                                                   aaCCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGa                                                                                 >  5:133569/1‑136 (MQ=255)
                                                        ggCAAGCACGAATAGCACGTTTTCGTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGCTCTTCCCTGCACGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCg                                                                         >  7:71052/1‑139 (MQ=255)
                                                        ggCAAGCACGAATAGCACGTTTTATTTTCTGAAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCg                                                                         >  7:81144/1‑139 (MQ=255)
                                                          cAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGCCAGACTCGcc                                                                       >  7:260692/1‑139 (MQ=255)
                                                              cACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGAGTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTc                                                                    >  5:76454/1‑138 (MQ=255)
                                                                         cGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCggg                                                        <  8:61860/139‑1 (MQ=255)
                                                                             ttcttttctGCAGCTGTACCCCATCGTGCCGAACCTCCAGATCTTCACTCCAAGCGGGCCCAATGAGCTGGCTAATTCAATTCCGCATTGATCCACCAGCATGTCACGGCCAGACTCGCCGTCAATGACGCCCGGAAtt                                                    >  3:69798/1‑139 (MQ=255)
                                                                                  ttctGCCGCTGTACCCCATCGTGCCGATACTACAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGTATTGTGCGAGCAGCATGTAACGGAAAGACTCGCGGTCAATGCCGCCGGGAATTTcta                                                <  4:69798/138‑3 (MQ=255)
                                                                                   agacaggCTGTACCCCATCGTGACGATCCTCCAGATCTTCCCTGCGAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGaca                                              <  5:175889/133‑1 (MQ=255)
                                                                                   agacaggCTGTACCCCATCCTGCCGATCCTCCAGATCTTCCCTGCGAGCGGGCCGAATGAGCTGGCGAATTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGaca                                              <  5:171795/133‑1 (MQ=255)
                                                                                    ctGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGaca                                              <  5:151128/138‑1 (MQ=255)
                                                                                    ctGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGaca                                              <  1:76002/138‑1 (MQ=255)
                                                                                        aGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCATGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGt                                         <  6:109395/139‑1 (MQ=255)
                                                                                        aGCTGTACCCCATCGTCCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGt                                         <  6:133569/139‑1 (MQ=255)
                                                                                         gcggTACCCCATCGTGCCGATCCTGCAGATCTTCCCTGCATGCGGGCAGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCTGCTGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGt                                         <  5:30720/135‑1 (MQ=255)
                                                                                                    tCGTGCCGATCCTCCAGATCTTCCCTGCGAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATTTCACGGACAGACTCGCCGTCAATTACGCCGGGAATTTCGCCACATGTGacaacctgtc                              >  6:175889/1‑133 (MQ=255)
                                                                                                    tCGTGCCGATCCTCCAGATCTTCCCTGCGAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCGCGGACAGACTCGCCGTCAATGCCGCCGGGTATTTAGACACAGGTGCCAACCggtct                             >  6:171795/1‑136 (MQ=255)
                                                                                                    tCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGATTTTCGACACAGGTGCCCACAGGCAt                             >  3:101891/1‑139 (MQ=255)
                                                                                                    tCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCAt                             >  8:88746/1‑139 (MQ=255)
                                                                                                    tCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCAt                             >  4:255494/1‑139 (MQ=255)
                                                                                                    tCGTGCCGATCCTCCAGATCTACCCTGCAAGCGGGCCGAATGAGCTGGCGATTTAAATGCCGAATTGAGCCAGCAGCATGTCACGGACAGCCTCGCCGTCAAAGCCGCCGTGAATTTAGACGACCGTGacaaccagca                              >  6:30720/1‑133 (MQ=255)
                                                                                                       tGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCAtttt                          >  5:275878/1‑139 (MQ=255)
                                                                                                         ccGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTg                         >  1:20284/1‑138 (MQ=255)
                                                                                                            agcctccAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCACTGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGc                     >  8:76579/3‑139 (MQ=255)
                                                                                                                  cAGCTCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCAt               <  8:81144/139‑1 (MQ=255)
                                                                                                                   aGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCAtc              <  5:81120/139‑1 (MQ=255)
                                                                                                                      tCTTCCCTGCAAGCGGGCCGAATGCGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCAGTTTGCAGCTGTCAtctc            >  8:241369/1‑138 (MQ=255)
                                                                                                                          cccTGCAAGCGGGCCGAATGCGCTGGCGATTTCAATGCCGACTTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCaaa        <  7:241369/138‑1 (MQ=255)
                                                                                                                          cccTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCaaa        <  3:184724/138‑1 (MQ=255)
                                                                                                                          cccTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCaaa        <  3:181252/138‑1 (MQ=255)
                                                                                                                           ccTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCAAAc       >  1:106796/1‑138 (MQ=255)
                                                                                                                           ccTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGAAAACAGGCATTTTGCAGCTGTCATCTCCAAAc       >  1:165841/1‑138 (MQ=255)
                                                                                                                               cAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGCCAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCAAACAATCt  >  5:134659/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                         
AGCGCTTTTCCTGCTTCGATTCCTGCGTTATGTCTGAGAAGCACAGCCTCTAACCCGGCAAGCACGAATAGCACGTTTTCTTTTCTGCAGCTGTACCCCATCGTGCCGATCCTCCAGATCTTCCCTGCAAGCGGGCCGAATGAGCTGGCGATTTCAATGCCGAATTGTGCCAGCAGCATGTCACGGACAGACTCGCCGTCAATGCCGCCGGGAATTTCGACACAGGTGACAACAGGCATTTTGCAGCTGTCATCTCCAAACAATCT  >  CP006881/3159072‑3159337

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