Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,551,361 G→A 8.5% V8V (GTG→GTA)  fdnI → formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,551,3610GA 8.5% 146.5 / 4.0 47V8V (GTG→GTA) fdnIformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (2/2);  total (22/25)
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.96e-01

TGGAAAGGCGCGTTGAAACCGCTGGCAGCGGCTGGCTTTATTGCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGC  >  NC_000913/1551217‑1551506
                                                                                                                                                |                                                                                                                                                 
tgGAAAGGCGCGTTGAAACCGCTGGCAGCGGCTGGCTTTATTGCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTACGTCGAAAATGATTGTgcgc                                                                                                                                                <  6:34610/148‑1 (MQ=255)
               aaaCCGCTGGCAGCGGCTGGCTTTATTGCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGa                                                                                                                                  <  7:263421/147‑1 (MQ=255)
                    gCTGGCAGCGGCTGGCTTTATTGCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCgcgc                                                                                                                            <  7:69315/148‑1 (MQ=255)
                         gaggggctggctTTATTGCCACTTTTGCCGGGTTGATTTTCCACTAAATCGGTATTGGCCAGATTAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTTTc                                                                                                                       <  5:87960/144‑1 (MQ=255)
                                   ctTTATTGCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgt                                                                                                             <  1:57153/148‑1 (MQ=255)
                                          gCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGAtt                                                                                                      >  7:246358/1‑148 (MQ=255)
                                          gCCACTTTTGCCGGGGTGATTTTCCACTACATAGTTATTGGCCAGAATAAGGAAGTGGACGATGACGAGGAGTATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGAt                                                                                                       <  2:71192/147‑1 (MQ=255)
                                           ccACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGGCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGAtt                                                                                                      <  7:322754/147‑1 (MQ=255)
                                            cACTTTTGCCGGGGTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTg                                                                                                    <  8:246358/148‑1 (MQ=255)
                                                tttGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGcttc                                                                                                >  2:187627/1‑148 (MQ=255)
                                                     cGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgtg                                                                                                            <  7:332839/131‑1 (MQ=255)
                                                     cGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgtg                                                                                                            >  8:332839/1‑131 (MQ=255)
                                                        gTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCtggt                                                                                        >  5:136644/1‑148 (MQ=255)
                                                        gTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTACGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgtggtt                                                                                                         <  3:2660/131‑2 (MQ=255)
                                                        gTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTACGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgtggtg                                                                                                         <  3:2657/131‑1 (MQ=255)
                                                        gTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTACGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgtggtg                                                                                                         >  4:2660/1‑131 (MQ=255)
                                                        gTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTACGCACCAAATTTATTGATCGCGCCTGTCACTGGACCgtggtg                                                                                                         >  4:2657/1‑131 (MQ=255)
                                                                tCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGcttctt                                                                                              <  3:35806/134‑1 (MQ=255)
                                                                tCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGcttctt                                                                                              >  4:35806/1‑134 (MQ=255)
                                                                         tCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGAttt                                                                       <  1:187627/148‑1 (MQ=255)
                                                                               ttGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGttct                                                                 >  5:24511/1‑148 (MQ=255)
                                                                                 ggCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGttctt                                                                <  6:136644/147‑1 (MQ=255)
                                                                                 ggCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGttctt                                                                <  6:24511/147‑1 (MQ=255)
                                                                                   cccGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGttcttctt                                                             <  3:209268/148‑1 (MQ=255)
                                                                                             gAAGTGGACGATGACGAGGAGGATCATCATTAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTg                                                   <  4:72533/148‑1 (MQ=255)
                                                                                               aGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAACTGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCa                                                 <  2:164343/148‑1 (MQ=255)
                                                                                                  ggACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATg                                              >  2:22799/1‑148 (MQ=255)
                                                                                                   gACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATgg                                             >  8:280171/1‑148 (MQ=255)
                                                                                                               gaggaTCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAAcc                                 >  1:61040/1‑148 (MQ=255)
                                                                                                                     catcatGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCtt                               <  3:205604/144‑1 (MQ=255)
                                                                                                                     catcatGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCtt                               >  4:205604/1‑144 (MQ=255)
                                                                                                                       tcatGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCa                                                 <  7:50025/124‑1 (MQ=255)
                                                                                                                       tcatGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCa                                                 >  8:50025/1‑124 (MQ=255)
                                                                                                                        catGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTAcg                        >  8:88846/1‑148 (MQ=255)
                                                                                                                          tGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCaaa                                   <  8:329139/135‑1 (MQ=255)
                                                                                                                          tGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCaaa                                   >  7:329139/1‑135 (MQ=255)
                                                                                                                           gAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGCCCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTGCTTCCCGACGCTGCAATGGCTGACGCAAACCTACGGTAcgccg                     >  8:123534/1‑148 (MQ=255)
                                                                                                                           gAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCa                                                 <  3:155109/120‑1 (MQ=255)
                                                                                                                           gAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCa                                                 >  4:155109/1‑120 (MQ=255)
                                                                                                                             gTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGcttctt                                                                                              >  1:20385/1‑73 (MQ=255)
                                                                                                                             gTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGcttctt                                                                                              <  2:20385/73‑1 (MQ=255)
                                                                                                                              tAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTAcgccg                     >  8:176725/1‑145 (MQ=255)
                                                                                                                               aaGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGttctt                                                                <  6:305310/101‑1 (MQ=255)
                                                                                                                               aaGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGttctt                                                                >  5:305310/1‑101 (MQ=255)
                                                                                                                                aGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCTGGGTGATTTGCTTCTTCCTGGTGGCGCGGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGa                 <  8:120330/147‑1 (MQ=255)
                                                                                                                                      aaaTGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACg          >  7:206791/1‑148 (MQ=255)
                                                                                                                                      aaaTGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACg          >  5:33725/1‑148 (MQ=255)
                                                                                                                                         tGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCAt       <  2:61040/148‑1 (MQ=255)
                                                                                                                                             tgtgCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGt                                                                    >  7:246160/1‑83 (MQ=255)
                                                                                                                                             tgtgCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGt                                                                    <  8:246160/83‑1 (MQ=255)
                                                                                                                                               tgCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGc  <  1:22799/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGGAAAGGCGCGTTGAAACCGCTGGCAGCGGCTGGCTTTATTGCCACTTTTGCCGGGTTGATTTTCCACTACATCGGTATTGGCCCGAATAAGGAAGTGGACGATGACGAGGAGGATCATCATGAGTAAGTCGAAAATGATTGTGCGCACCAAATTTATTGATCGCGCCTGTCACTGGACCGTGGTGATTTGCTTCTTCCTGGTGGCGCTGTCCGGGATTTCGTTCTTCTTCCCGACGCTGCAATGGCTGACGCAAACCTTCGGTACGCCGCAGATGGGACGCATTTTGC  >  NC_000913/1551217‑1551506

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