Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 65,812 T→C 9.3% intergenic (‑32/+43) polB ← / ← araD DNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091365,8120TC 9.3% 130.2 / 4.8 43intergenic (‑32/+43)polB/araDDNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATAT  >  NC_000913/65667‑65959
                                                                                                                                                 |                                                                                                                                                   
cGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtat                                                                                                                                                   <  3:937462/148‑1 (MQ=255)
                tGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                   >  1:67941/1‑148 (MQ=255)
                         ccGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                         >  4:203106/1‑133 (MQ=255)
                         ccGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                         <  3:203106/133‑1 (MQ=255)
                                  gTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGCATAGCCTgg                                                                                                                                            >  1:1085888/1‑121 (MQ=255)
                                  gTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGCATAGCCTgg                                                                                                                                            <  2:995730/121‑1 (MQ=255)
                                  gTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGCATAGCCTgg                                                                                                                                            <  2:1085888/121‑1 (MQ=255)
                                  gTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGCATAGCCTgg                                                                                                                                            >  1:995730/1‑121 (MQ=255)
                                              aaGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                         >  2:532111/1‑112 (MQ=255)
                                              aaGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                         <  1:532111/112‑1 (MQ=255)
                                              aaGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTg                                                                                                     <  2:67941/148‑1 (MQ=255)
                                                 gagaCTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                     >  4:594309/1‑113 (MQ=255)
                                                 gagaCTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                     <  3:594309/113‑1 (MQ=255)
                                                  agaCTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtat                                                                                                                                                   >  1:988816/1‑98 (MQ=255)
                                                  agaCTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtat                                                                                                                                                   <  2:988816/98‑1 (MQ=255)
                                                                tGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGcc                                                                                                                                               <  1:464578/88‑1 (MQ=255)
                                                                tGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGcc                                                                                                                                               >  2:464578/1‑88 (MQ=255)
                                                                       tgtCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGt                                                                                                                       >  4:1053985/1‑105 (MQ=255)
                                                                       tgtCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGt                                                                                                                       <  3:1053985/105‑1 (MQ=255)
                                                                            cGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGa                                                                                                                                  >  4:953936/1‑89 (MQ=255)
                                                                            cGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGa                                                                                                                                  <  3:953936/89‑1 (MQ=255)
                                                                              ccAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGATAGTCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGa                                                                                                                                  <  4:701947/87‑1 (MQ=255)
                                                                              ccAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTCGa                                                                                                                                  >  3:701947/1‑87 (MQ=255)
                                                                                          tAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTaa                                                                                              <  3:968482/111‑1 (MQ=255)
                                                                                          tAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTaa                                                                                              >  4:968482/1‑111 (MQ=255)
                                                                                             gATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtata                                                                                                                                                  <  3:667670/56‑1 (MQ=255)
                                                                                             gATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtata                                                                                                                                                  >  4:667670/1‑56 (MQ=255)
                                                                                               tAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAg                                                                                                                                                 <  4:810433/55‑1 (MQ=255)
                                                                                               tAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAg                                                                                                                                                 >  3:810433/1‑55 (MQ=255)
                                                                                               tAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                         >  1:954734/1‑63 (MQ=255)
                                                                                               tAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                         <  2:954734/63‑1 (MQ=255)
                                                                                                               cACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGccc                                                                                                  >  3:889164/1‑86 (MQ=255)
                                                                                                               cACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGccc                                                                                                  <  4:889164/86‑1 (MQ=255)
                                                                                                                           ccATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATAc                                                                                                                >  3:747149/1‑60 (MQ=255)
                                                                                                                           ccATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATAc                                                                                                                <  4:747149/60‑1 (MQ=255)
                                                                                                                              tCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCatat                                    <  1:609560/133‑1 (MQ=255)
                                                                                                                              tCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCatat                                    >  2:609560/1‑133 (MQ=255)
                                                                                                                               cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                         >  1:1114823/1‑63 (MQ=255)
                                                                                                                               cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                         <  2:1114823/63‑1 (MQ=255)
                                                                                                                               cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagc                                                    <  4:50272/116‑1 (MQ=255)
                                                                                                                               cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagc                                                    >  3:50272/1‑116 (MQ=255)
                                                                                                                                 aaaaaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtg                                            <  4:1121441/122‑1 (MQ=255)
                                                                                                                                 aaaaaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtg                                            >  3:1121441/1‑122 (MQ=255)
                                                                                                                                     aCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTa                                                                         <  4:687636/89‑1 (MQ=255)
                                                                                                                                     aCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTa                                                                         >  3:687636/1‑89 (MQ=255)
                                                                                                                                        aGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACg           >  4:510604/1‑148 (MQ=255)
                                                                                                                                             tGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTa                                                                                                        <  4:674412/50‑1 (MQ=255)
                                                                                                                                             tGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTa                                                                                                        >  3:674412/1‑50 (MQ=255)
                                                                                                                                                 tataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGcc                                                                                                   <  3:841335/51‑1 (MQ=255)
                                                                                                                                                 tataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGcc                                                                                                   >  4:841335/1‑51 (MQ=255)
                                                                                                                                                 tataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAatat  <  3:510604/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATAT  >  NC_000913/65667‑65959

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