Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,191 C→A 21.7% intergenic (+200/+207) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,1910CA21.7% 71.2 / 6.0 49intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (13/24);  new base A (9/3);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 2.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGC  >  NC_000913/1552043‑1552338
                                                                                                                                                    |                                                                                                                                                   
tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                     <  1:1125243/149‑1 (MQ=31)
        aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                             <  1:681427/149‑1 (MQ=34)
                 gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                    <  2:583356/149‑1 (MQ=40)
                  cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTaaaatg                                                                                                                                   >  2:975436/1‑147 (MQ=37)
                      gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTc                                                                                                                               <  1:1109526/149‑1 (MQ=255)
                            gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                         <  1:101404/149‑1 (MQ=255)
                            gCCGGATGCCGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                         <  1:57825/149‑1 (MQ=255)
                                 aTGCGGCGTGAACGCCTAATCCGGCATGAAACCCCTTCCACTCCAAGAGATTGCAGCGAACGTGGAGGCCTGATACGCGTAGCGCAGCAGGCAATGTTGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGGCATTTTCATAtt                                                                                                                    <  2:96842/149‑1 (MQ=14)
                                   gcggctTGAAAGGCTTATCCGGCATGAAAACCCTTCAAATCCCATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGt                                                                                                                  <  1:778908/149‑1 (MQ=17)
                                   gcggcgTGTACGCCTTATCCGGCATGAAAACCCTTCATATCCAATAGATTGCAGTGCACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                  <  1:836273/149‑1 (MQ=255)
                                   gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGGAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                  <  2:553936/149‑1 (MQ=255)
                                              gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                                                                                                       <  2:394932/149‑1 (MQ=255)
                                                cTTATCCGGCATGAAAACCCTTCAAATCCATTATATTTTTGTGAACGTGTAGGTCTGATTAGCGTAGCGAATCAGTAAATGTTTCGTTTTTCTTAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                                                                                                     <  2:835938/149‑1 (MQ=18)
                                                      cGGCATGATAACCCTTCAAATCCAATAGATTGCATTGTCGGTGTAGTCCTGATTAGCGTATCTCTTCAGGGAATGTTGCGTTTGTCATAAGTTTATAATGGAGCTGTAAAAGGCGTCATTTTCATATTGTAGTCAACGTAGGCTTTGtt                                                                                               <  2:322268/149‑1 (MQ=18)
                                                            gAAAACCCTTCAAATCCAATAGATTGCCGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                         <  1:702508/149‑1 (MQ=255)
                                                             aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCg                                                                                        <  2:1165048/149‑1 (MQ=255)
                                                                 cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATg                                                                                    >  1:527198/1‑149 (MQ=255)
                                                                 cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTAATATTGATGACAACGTAGGCTTTGTTCATTCCGGAtc                                                                                    >  2:980097/1‑148 (MQ=18)
                                                                       aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGACATTTTAATATTGTATACAAAGAAGGCTTTTTTCATGCCTGATgatgagt                                                                              >  1:682388/1‑143 (MQ=17)
                                                                         aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                            >  1:144429/1‑149 (MQ=255)
                                                                         aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAACGGGGCATTTTAAGATTGAAGACACCGTATGCGTTTTTCATTCCCGAAgcctcgtta                                                                            >  1:779588/1‑142 (MQ=14)
                                                                             aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTATAATGAGTCATTTTATTATTGTATACAACGTAGGCTTTGTTCATTCCGGATCCGGCGTTAAccc                                                                        >  2:275192/1‑147 (MQ=17)
                                                                              aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTGATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGaaagac                                                                       >  1:13021/1‑145 (MQ=21)
                                                                              aTAGATTGCACTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                       <  2:445967/149‑1 (MQ=255)
                                                                               tAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                      <  2:1029684/149‑1 (MQ=255)
                                                                                       aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCGTCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCTTATTGTAGAAAACGTCGGCTTTGTTCATACAGAATCCGGCTTTAACGCCATATGCGGc                                                              >  1:145667/1‑149 (MQ=18)
                                                                                        gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                             <  1:703019/149‑1 (MQ=255)
                                                                                         tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                            <  2:527198/149‑1 (MQ=255)
                                                                                            aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGCGTCATTTTAATATTGTAGACAACGAAGGCTTTTTTCATGCAGGTAGCGACTTGAACGCATAATCCGGCATGaa                                                         >  2:338905/1‑149 (MQ=21)
                                                                                                   gtcatGATAAGCGTAGCGAATCAGGCAATGTTTCGTTTTGGATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAGGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                  <  2:96151/145‑1 (MQ=17)
                                                                                                    gCCTGATAAGCGTAGCGCTTCAGGCAATGTTGCGATTGTCATCAGTTTCAAAAGGCGCTGTAAAAGCCGTCATTTTCATATTGTAGACAACGTAGGCTTTCCGCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                 <  1:391596/149‑1 (MQ=21)
                                                                                                    gCCGGATAAGCGTAGCTCATCATTCAATGTTGCTTTTGTAATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                 <  2:241457/149‑1 (MQ=17)
                                                                                                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                >  1:956361/1‑149 (MQ=255)
                                                                                                       ttaTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaaa                                              <  2:186552/147‑1 (MQ=255)
                                                                                                            aGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                         <  2:381485/149‑1 (MQ=255)
                                                                                                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                         <  2:144429/149‑1 (MQ=255)
                                                                                                               ggaGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAg                                      <  1:477636/147‑1 (MQ=255)
                                                                                                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATATTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGAATCAAAACCCTTCAAAACCAATAg                                      >  1:651333/1‑149 (MQ=18)
                                                                                                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTAATGCCGGATGCGTCTTGAACGCCTTATACGGCATGAAAAAACTTCAAATCAAATat                                      >  1:442317/1‑148 (MQ=18)
                                                                                                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTACAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGAAGTCTTTGTTCATTCCGAATGCGGCGTGAACGCCATATACGGCATGAAAAACCTTCaaataaaatat                                      >  1:461920/1‑148 (MQ=37)
                                                                                                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGAGTCATTTTAATATTGTAGAAAAAGTAGGCATTTATCATGACGGAGGCGTAGTGAAAGCCTTAACCGGCATGAAAAACCGTCTATTCCAAttt                                      >  2:56797/1‑147 (MQ=12)
                                                                                                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGAAGAAAAAGTAGGCATTGATCATGCCCGATGCGTAGTGAAAGACTTTTCCGGCATGAAAAACCTTCAACTCCAATaa                                      >  2:141174/1‑148 (MQ=21)
                                                                                                                tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTATTATTGTATACAAAGTAGGCTTTTTTCATGCCGGATGCGCCGTGAACGCCTTAAACTGCATTAACAGACTTAAAATCAATTAGa                                     >  1:666904/1‑149 (MQ=21)
                                                                                                                           ttaaTGTTGATTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAc                          <  2:604061/147‑1 (MQ=255)
                                                                                                                             aaTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                        >  1:972109/1‑149 (MQ=255)
                                                                                                                                gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                     <  2:956361/149‑1 (MQ=255)
                                                                                                                                     gTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGACCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGACt                >  2:1182359/1‑149 (MQ=255)
                                                                                                                                            atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc         >  2:300445/1‑149 (MQ=39)
                                                                                                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAg     >  1:1126885/1‑149 (MQ=35)
                                                                                                                                                   tCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAACCGGCATGAAAACCCTTCAAATCCAAAAGAATGCAGTGAACGAGAAGGCCTGATAAGCGTAgcgc  >  2:168101/1‑149 (MQ=11)
                                                                                                                                                    |                                                                                                                                                   
TATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGC  >  NC_000913/1552043‑1552338

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