Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,532,632 G→A 9.5% intergenic (‑193/‑186) yhgE ← / → pck DUF4153 family putative inner membrane protein/phosphoenolpyruvate carboxykinase [ATP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,532,6320GA 9.5% 120.3 / 4.7 42intergenic (‑193/‑186)yhgE/pckDUF4153 family putative inner membrane protein/phosphoenolpyruvate carboxykinase [ATP]
Reads supporting (aligned to +/- strand):  ref base G (23/15);  new base A (2/2);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-02

TTGCTTTTTTCAGAAAATAATCAAAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAACACCCCCAAAAAGACTTTA  >  NC_000913/3532495‑3532780
                                                                                                                                         |                                                                                                                                                    
ttGCTTTTTTCAGAAAATAATCAAAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTc                                                                                                                                           >  3:122641/1‑149 (MQ=255)
ttGCTTTTTTCAGAAAATAATCAAAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTc                                                                                                                                           <  4:287690/149‑1 (MQ=255)
    ttttttCAGAAAATAATCATAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATAAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCcaca                                                                                                                                       >  4:300071/1‑149 (MQ=255)
    ttttttCAGAAAATAATCAAAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCcaca                                                                                                                                       >  4:463488/1‑149 (MQ=255)
                    tCAAAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTc                                                                                                                       >  1:44556/1‑149 (MQ=255)
                       aaaaaaGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGTGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGc                                                                                                                    <  1:85633/149‑1 (MQ=255)
                              ttAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGataata                                                                                                             <  2:44556/149‑1 (MQ=255)
                              ttAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGataata                                                                                                             >  2:235307/1‑149 (MQ=255)
                                      gtgAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGt                                                                                                     >  2:36288/1‑149 (MQ=255)
                                        gAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGttt                                                                                                   >  2:389833/1‑149 (MQ=255)
                                           tCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGACAACATTTTCCACACCTAAAATACTATTTCTGCGATAATAGCAACCGTTTCGt                                                                                                >  3:363303/1‑149 (MQ=255)
                                                  tttACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGa                                                                                         >  3:152673/1‑149 (MQ=255)
                                                  tttACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCAGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTTAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGa                                                                                         >  1:388400/1‑149 (MQ=255)
                                                       cGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCAc                                                                                    >  3:202568/1‑149 (MQ=255)
                                                        ggTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACg                                                                                   <  3:300071/149‑1 (MQ=255)
                                                         gTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCAc                                                                                    <  2:185023/147‑1 (MQ=255)
                                                         gTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCAc                                                                                    >  1:185023/1‑147 (MQ=255)
                                                         gTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACgg                                                                                  >  4:1236/1‑149 (MQ=255)
                                                          ttGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGa                                                                                 <  1:235307/149‑1 (MQ=255)
                                                            gAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGc                                                                                                                    >  1:391740/1‑112 (MQ=255)
                                                            gAATTTGCATCAAGTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGc                                                                                                                    <  2:391740/112‑1 (MQ=255)
                                                                         tttCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAata                                                                  >  3:172645/1‑149 (MQ=255)
                                                                              ttcaGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAAt                                                             <  4:363303/149‑1 (MQ=255)
                                                                                 aGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTtc                                                          >  2:415504/1‑149 (MQ=255)
                                                                                 aGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTtc                                                          >  4:256422/1‑149 (MQ=255)
                                                                                 aGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTtc                                                          >  2:99752/1‑149 (MQ=255)
                                                                                   gAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTtctc                                                        <  2:344051/149‑1 (MQ=255)
                                                                                      tGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAg                                                     <  3:1236/149‑1 (MQ=255)
                                                                                          aTTCCACTCACAATATTCCCGCCATATAAACCAACATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATGACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATAc                                                 <  1:389833/149‑1 (MQ=255)
                                                                                               aCTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTaaa                                            >  2:418135/1‑149 (MQ=255)
                                                                                                  cacaATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCt                                         <  1:415504/149‑1 (MQ=255)
                                                                                                            ccGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTg                               <  4:152673/149‑1 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAacac                  >  3:446417/1‑149 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAacac                  >  1:289593/1‑149 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTGTGTCAAATATGAATTTCTTCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAacac                  <  2:388400/149‑1 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTAAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAaca                   <  3:130373/148‑1 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTAAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAaca                   >  2:467549/1‑148 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTAAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAaca                   >  4:130373/1‑148 (MQ=255)
                                                                                                                         cAAGATTTAACCTTTTAAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAaca                   <  1:467549/148‑1 (MQ=255)
                                                                                                                                 aaCCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTAt                                       <  2:21988/120‑1 (MQ=255)
                                                                                                                                 aaCCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTAt                                       >  1:21988/1‑120 (MQ=255)
                                                                                                                                     ttttGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAACACCCCCAAAAAGAc      >  2:77126/1‑149 (MQ=255)
                                                                                                                                         gagaACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAACACCCCCAAAAAGACTTTa  <  1:36288/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TTGCTTTTTTCAGAAAATAATCAAAAAAAGTTAGCGTGGTGAATCGATACTTTACCGGTTGAATTTGCATCAATTTCATTCAGGAATGCGATTCCACTCACAATATTCCCGCCATATAAACCAAGATTTAACCTTTTGAGAACATTTTCCACACCTAAAATGCTATTTCTGCGATAATAGCAACCGTTTCGTGACAGGAATCACGGAGTTTTTTGTCAAATATGAATTTCTCCAGATACGTAAATCTATGAGCCTTGTCGCGGTTAACACCCCCAAAAAGACTTTA  >  NC_000913/3532495‑3532780

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