Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,757 G→T 16.0% intergenic (+350/‑149) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,7570GT16.0% 72.8 / 6.5 25intergenic (+350/‑149)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (11/10);  new base T (2/2);  total (13/12)
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.72e-01

TAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAAT  >  NC_000913/4235654‑4235901
                                                                                                       |                                                                                                                                                
tAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTtat                                                                                                                                          >  1:913331/1‑112 (MQ=255)
tAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTtat                                                                                                                                          <  2:913331/112‑1 (MQ=255)
                                       gCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGa                                                                                                                                    <  4:442236/79‑1 (MQ=255)
                                       gCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGa                                                                                                                                    >  3:442236/1‑79 (MQ=255)
                                       gCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTgtgt                                                              >  3:1393030/1‑149 (MQ=255)
                                         aGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGta                                                             <  4:826437/148‑1 (MQ=255)
                                         aGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGta                                                             >  3:826437/1‑148 (MQ=255)
                                                   cTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAAATcccc                                                                                           >  1:1212884/1‑108 (MQ=255)
                                                   cTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAAATcccc                                                                                           <  2:1212884/108‑1 (MQ=255)
                                                    ttttGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAg                                                                    <  4:532738/130‑1 (MQ=255)
                                                    ttttGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAg                                                                    >  3:532738/1‑130 (MQ=255)
                                                                       ttATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTtgtg                                                               >  3:1402955/1‑116 (MQ=255)
                                                                       ttATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTtgtg                                                               <  4:1402955/116‑1 (MQ=255)
                                                                                        tatcaaaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt                                   <  3:153706/123‑1 (MQ=255)
                                                                                        tatcaaaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt                                   >  4:153706/5‑127 (MQ=255)
                                                                                         tttAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTaa                                                                                                 >  4:1137269/1‑64 (MQ=255)
                                                                                         tttAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTaa                                                                                                 <  3:1137269/64‑1 (MQ=255)
                                                                                         tttAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAAt            >  3:425409/1‑149 (MQ=255)
                                                                                            aaaTATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt                                   <  3:597063/123‑1 (MQ=255)
                                                                                            aaaTATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt                                   >  4:597063/1‑123 (MQ=255)
                                                                                             aaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAAtctc                        >  3:456339/1‑133 (MQ=255)
                                                                                             aaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAAtctc                        <  4:456339/133‑1 (MQ=255)
                                                                                                   aGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAAt  <  4:1393030/149‑1 (MQ=255)
                                                                                                    gTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCa                                                                            <  4:1134603/74‑1 (MQ=255)
                                                                                                    gTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCa                                                                            >  3:1134603/1‑74 (MQ=255)
                                                                                                       |                                                                                                                                                
TAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAAT  >  NC_000913/4235654‑4235901

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