Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,509,497 C→G intergenic (‑71/‑133) glk ← / → yfeO glucokinase/putative ion channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,4970CG97.9% 142.9 / NA 48intergenic (‑71/‑133)glk/yfeOglucokinase/putative ion channel protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (19/28);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAT  >  NC_000913/2509358‑2509631
                                                                                                                                            |                                                                                                                                      
ggCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGcc                                                                                                                                     <  2:125230/143‑1 (MQ=255)
   aaTATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGa                                                                                                                                  >  1:875299/1‑143 (MQ=255)
        cacacaGAGCAAGACGTGCGTTGGTGCCGCCCACATCCCCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAAGAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAaggc                                                                                                                             <  2:337070/143‑1 (MQ=255)
         acacaGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAaggca                                                                                                                            >  2:159801/1‑143 (MQ=255)
           acaGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAaggcagg                                                                                                                          >  1:849669/1‑143 (MQ=255)
               agCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATCACACCAAGGCCGGAACGCCGGGCTcc                                                                                                                      >  1:48617/1‑143 (MQ=255)
                   aGACGTGCGTTGGTGCCGCCCACATTACCGACTAATGCATACTCTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAACTACCTGGGGACATAACAACCAGGCCGGACAGGCAGGCTCCCTGt                                                                                                                  >  2:687132/1‑143 (MQ=255)
                       gTGCGTTGGTGCCGCCCCCACCACCGACTAATGCACACTCTGTCATTCCTCAACTGCTCCGCCAAAG‑CCAAAATAACGCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAat                                                                                                              <  1:307154/143‑1 (MQ=255)
                        tGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAata                                                                                                             <  1:18562/143‑1 (MQ=255)
                         gCGTTGGTGCCGCCCACACCACCGACTAATGCATACTTTCTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAatat                                                                                                            <  1:650224/143‑1 (MQ=255)
                           gTTGGTGCCTCCCACAGCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAAGTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCg                                                                                                          <  2:697401/143‑1 (MQ=255)
                             tGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGAt                                                                                                        <  1:808883/143‑1 (MQ=255)
                                 gcggccCAGATCACCGACTGAGGCATACGTGGTCAGTGTTCAACGGCTCCGCTAAAG‑TCAAGATAATTCTGGCTCAGACTGTATATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTgg                                                                                                    <  1:684017/140‑1 (MQ=255)
                                     cccACATCACCTACTAATGCATACTTTGTCATTCTTCAACTGCTCTG‑TATTGATCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTca                                                                                                <  1:216862/143‑1 (MQ=255)
                                           tccccGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAAt                                                                                          <  1:805550/140‑1 (MQ=255)
                                           tacccGACTAGTGCATACTTTGTCATTCTTCAACTGCGCCGCTAAAG‑TCAAAATAATTCTTTCTCACACGGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAAt                                                                                          <  2:436224/140‑1 (MQ=255)
                                               cGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACt                                                                                      <  1:467377/143‑1 (MQ=255)
                                                 aCTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACttt                                                                                    >  1:116473/1‑143 (MQ=255)
                                                  cTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAGCAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTa                                                                                   >  1:442111/1‑143 (MQ=255)
                                                       gCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGtt                                                                              >  1:474246/1‑143 (MQ=255)
                                                           aCTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAg                                                                          <  1:765493/143‑1 (MQ=255)
                                                                 tCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCaa                                                                    >  2:297685/1‑143 (MQ=255)
                                                                 tCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCaa                                                                    >  1:164331/1‑143 (MQ=255)
                                                                   aTTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAAtt                                                                  >  1:103068/1‑143 (MQ=255)
                                                                   aTTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGACACACGATTACTTTATGGTCTCCGCagctcta                                                                  >  1:943206/1‑137 (MQ=255)
                                                                      cttcAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCa                                                               >  2:373523/1‑143 (MQ=255)
                                                                        tcAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGTCCGGAAAGGCAGGCTCCCTGTAAATATAGATCTGGGTCACACAATTACCTTATCGTTTCAGCACCAATTGCAGc                                                             >  2:352161/1‑143 (MQ=255)
                                                                          caCTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGa                                                           <  1:477334/142‑1 (MQ=255)
                                                                             tGCCCCGCTAAAG‑CCACAACAAGTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGTCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGc                                                        <  1:15385/143‑1 (MQ=255)
                                                                                tCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCtt                                                     >  1:840889/1‑143 (MQ=255)
                                                                                   gCTAATG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCttttt                                                  <  2:253964/143‑1 (MQ=255)
                                                                                    cTAAAG‑TCAAAATAATTCTTTCTCACACAGCCCGGAGCTGCGCGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTg                                                 <  2:925091/143‑1 (MQ=255)
                                                                                           tCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCt                                           <  1:159801/143‑1 (MQ=255)
                                                                                           tCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGACGCCTCTTCGCAAGCt                                           >  2:7837/1‑143 (MQ=255)
                                                                                              acataatTGTGTCTCACCCGTTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTggg                                        <  2:813257/141‑1 (MQ=255)
                                                                                              aaaTAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGGTggg                                        >  1:736974/1‑143 (MQ=255)
                                                                                                      tattctCACATGGTAAATTCCTGGGGGCATAACGACACGGCCGGAATGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTaa                                <  2:624132/141‑1 (MQ=255)
                                                                                                       tttctCACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAg                               <  2:474246/143‑1 (MQ=255)
                                                                                                          ctcACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTAt                            <  1:373523/143‑1 (MQ=255)
                                                                                                           tcACACTGTAAATACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCCCCCTGTAAATATCGATCTGGGTCACACAATTACTTTACCGTTTCAGCACAAATCGCAGCGATGCCTCTTTGCAAGCTGGGCACACTAAGTagc                           >  2:248804/1‑141 (MQ=255)
                                                                                                              cattGTAAATTCCTTGGGGCATAACAACAAGTCTGTAAAGGTAGGCTCCCTGTAAATATCGATCTGGGTCACATAATTACTTTATCTTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGa                        <  2:663947/140‑1 (MQ=255)
                                                                                                                cTGGAAAGACCCGGGTGCATAACCACTGGGCCGGAAATGTGGTCGTCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGAcc                      <  2:925701/143‑1 (MQ=255)
                                                                                                                       tACCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATa               <  2:116473/143‑1 (MQ=255)
                                                                                                                        aCCTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATAGATATGGGGCACAAACTTACTTTACTGATTCAGCACCAATTGCAGCGATGCCATTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATac              >  2:833012/1‑142 (MQ=255)
                                                                                                                         tcTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAg             <  2:103068/142‑1 (MQ=255)
                                                                                                                          cTGGGGGCATAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATaatg            >  2:760521/1‑141 (MQ=255)
                                                                                                                               ggCATAACACCAGGGCCGGAAAGGAAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAg       <  1:157955/143‑1 (MQ=255)
                                                                                                                                   tAACAACAAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa   <  1:7837/143‑1 (MQ=255)
                                                                                                                                    aacaacaagGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAAAAt  >  1:5866/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                      
GGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAG‑TCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAT  >  NC_000913/2509358‑2509631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑