Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,769,067 A→G 16.7% intergenic (‑382/‑7) ydiJ ← / → ydiK putative FAD‑linked oxidoreductase/putative transporter YdiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,0670AG16.7% 58.5 / 5.0 24intergenic (‑382/‑7)ydiJ/ydiKputative FAD‑linked oxidoreductase/putative transporter YdiK
Reads supporting (aligned to +/- strand):  ref base A (8/12);  new base G (2/2);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GTGAACTCCAGGAGCTTACAATAGTAAGTGACTGGGGTGAACGAACGTAGCCGCAGCACATGCAACTTGAAATACGACGAGTAAATCGTTTGCGTGTTGCCTGAGTTGTTGTACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATGGGCCGGTACGGTGGTTATCGCCACCTGGCCGGTATTGTTACGT  >  NC_000913/1768815‑1769235
                                                                                                                                                                                                                                                            |                                                                                                                                                                        
ttGAACTCCAGGGTCTTAAAATATTAAGGGACTGGGGGGAACGAACGTAGCCGCAGAACATGCAACTTGAAATACGACGAGTAAATCGTTTGCTTGTTGCCTGAGTTGTTGTACCACATTTTTTTCTACCACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGATTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGAACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCg                                                                                                                                                      <  1:159446/272‑1 (MQ=255)
               ttACAATAGTAAGTGACTGGGGTGAACGAACGTAGCCGCAGCACATGCAACTTGAAATACGACGAGTAAATCGTTTGCGTGTTGCCTGAGTTGTTGTACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCCCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGCTTTGACGCACAAAAAAGGGGAAAGGGAGGTAAGGTAAATGTTCGTTAGCACAgaggggg                                                                                                                                       >  2:128665/1‑267 (MQ=255)
                  aaaTAGTAAGTGACTGGGGTGAACGAACGTAGGCGCAGCACATGCAACTGGAAATACGACGAGTAAATCGTTTGCGTGTTGCCTGAGTTTTTGTACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTAATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAGGTAGTTATGGTAAATGTTCg                                                                                                                                                      <  1:128665/254‑1 (MQ=255)
                                                       gcaCATGCAACTTGAAATACGACGAGTAAATCGTTTGCGTGTTGCCTGAGTTGTTGTACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGa                                                                                              <  1:153457/274‑1 (MQ=255)
                                                                                                        gttgttGTACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTTACGCACAAAAGAGGTGAAAcgta                                                                                                                                                                       >  2:258170/1‑152 (MQ=255)
                                                                                                           gttgtACCACATTTTTTTCTACCACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGtt                                                                                                                                                                    <  1:258170/152‑1 (MQ=255)
                                                                                                                    caTTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGt                                                                                                                                                               >  2:108663/1‑148 (MQ=255)
                                                                                                                     aTTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATggg                                           <  1:51312/263‑1 (MQ=255)
                                                                                                                       ttttttCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTaaa                                                                                                                                                            <  1:108663/148‑1 (MQ=255)
                                                                                                                                                   gggCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATTGGCCGGGACGGTGGGTATCGCACCCCGGTCGGTTTTGtttatg  >  2:101685/1‑270 (MQ=255)
                                                                                                                                                        gAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCcatca                                                                                                 >  1:64131/1‑174 (MQ=255)
                                                                                                                                                           tCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGa                                                                                              <  2:64131/174‑1 (MQ=255)
                                                                                                                                                                  gTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTGCGCGTGTTCACGTCGCGTCGACGATGTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTAGTCTCGGCTTTGc                                                <  1:187459/213‑1 (MQ=255)
                                                                                                                                                                              tGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTc                                                                                                                                                    >  2:204580/1‑101 (MQ=255)
                                                                                                                                                                              tGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATGGGCCGGTACggtggt                              >  1:37623/1‑219 (MQ=255)
                                                                                                                                                                                 aaTTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGc                                                                                                                                                 <  1:204580/101‑1 (MQ=255)
                                                                                                                                                                                 aaTTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATGGGCCGGTACGGTGGTTAt                           <  2:37623/219‑1 (MQ=255)
                                                                                                                                                                                                 gCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAGGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGAtt                                                                                            >  2:235513/1‑138 (MQ=255)
                                                                                                                                                                                                    ggTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAGGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATtgtg                                                                                         <  1:235513/138‑1 (MQ=255)
                                                                                                                                                                                                    ggTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCaa                                                                    >  2:220702/1‑159 (MQ=255)
                                                                                                                                                                                                      tGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGt                                                                                                                                                         >  1:212287/1‑72 (MQ=255)
                                                                                                                                                                                                      tGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATGGGCCGGTACgg                                  >  2:39361/1‑191 (MQ=255)
                                                                                                                                                                                                       gAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAAccc                                                                 <  1:220702/159‑1 (MQ=255)
                                                                                                                                                                                                         gTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCg                                                                                                                                                      <  2:212287/72‑1 (MQ=255)
                                                                                                                                                                                                         gTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATGGGCCGGTACggtgg                               <  1:39361/191‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                        
GTGAACTCCAGGAGCTTACAATAGTAAGTGACTGGGGTGAACGAACGTAGCCGCAGCACATGCAACTTGAAATACGACGAGTAAATCGTTTGCGTGTTGCCTGAGTTGTTGTACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGAAAAGTAGTTATGGTAAATGTTCGTCAGCCCAGGGATGTCGCACAAATTCTGCTTTCGGTGCTGTTTTTAGCCATCATGATTGTGGCATGTCTGTGGATTGTTCAACCCTTTATTCTCGGCTTTGCATGGGCCGGTACGGTGGTTATCGCCACCTGGCCGGTATTGTTACGT  >  NC_000913/1768815‑1769235

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑