Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,854 C→T intergenic (‑64/‑270) sthA ← / → fabR pyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,8540CT100.0% 159.5 / NA 50intergenic (‑64/‑270)sthA/fabRpyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/27);  total (23/27)

CCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAA  >  NC_000913/4160715‑4160992
                                                                                                                                           |                                                                                                                                          
ccgggCCCTTTGCAGCGCCTTCGCGGCCGGGGCCGGAACCTATTACGATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCaa                                                                                                                                         <  2:355456/140‑1 (MQ=255)
ccAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAAGCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCaa                                                                                                                                         <  1:919871/143‑1 (MQ=255)
 cAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAt                                                                                                                                        <  2:466909/143‑1 (MQ=255)
 cAGGCCCATTGCAGCCCCTTCGCCGCCGGGGCCGCAACCTATTACTATGGCATCGTAATCGTAGGAATGTGTCATGCTAAGTCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAt                                                                                                                                        <  1:385254/143‑1 (MQ=255)
  aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGCTAAGTAAAtt                                                                                                                                       >  2:685032/1‑143 (MQ=255)
  aGGCCCATTGCAGCGCCCTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAtt                                                                                                                                       <  1:161986/143‑1 (MQ=255)
      ccATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCCGGTAAGTCAATTGATc                                                                                                                                   >  2:166240/1‑143 (MQ=255)
                      gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTgg                                                                                                                   <  1:453391/143‑1 (MQ=255)
                      gccgccGGCGCGGCAACCTTTTACTAACGCATCGTGATCGGAGGTATTTGGCATGTTAGGGCTTACATTCTCTTATTTATATTAGCAACAGAATAGTAACATTTTATCGCGGGTAAGTCAATTGTTCCCCGTCATTTATCTgg                                                                                                                   <  2:1032176/143‑1 (MQ=255)
                        cgccgGGGCCGGAACCTATTACTATTGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCt                                                                                                                 <  1:64852/143‑1 (MQ=255)
                         gccgGGGCCGGAACCTAGTACTATGGCATTGTAATCGTAGGAAGGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCta                                                                                                                <  2:722807/143‑1 (MQ=255)
                          ccgGGGCCGGAATCTAGTACTGTGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCtat                                                                                                               <  2:243222/143‑1 (MQ=255)
                             gggCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATc                                                                                                            <  2:431300/143‑1 (MQ=255)
                               gCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAGGTGGCATGGTAGGGCTTACCTGTTCTTATAGATAAAAGCAACAGAATGGTAACAGTTTATCGCGGGTAAGTACATTGATCCCCGTCATTGATCTGGCTATATCCt                                                                                                          >  2:722299/1‑143 (MQ=255)
                                 cGGAACCCATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGa                                                                                                        <  1:96737/143‑1 (MQ=255)
                                  ggcaCCTATTACTATGGCATCGTAATCGTAGGAATATGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAg                                                                                                       <  1:1032712/140‑1 (MQ=255)
                                   gAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGCAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTCTCGCGGGTACGTCAATTGATCCCCGTCATTTATCCGGCTACATCCTGacc                                                                                                      >  2:104386/1‑141 (MQ=255)
                                                        taatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCCGTTTc                                                                                 >  1:788990/1‑143 (MQ=255)
                                                           tcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTa                                                                              <  2:624618/143‑1 (MQ=255)
                                                            cgtaGGAATGTGTCATGGTATGGCTTACCTTTTCTTAGACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTATGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTAc                                                                             <  2:901519/143‑1 (MQ=255)
                                                              taGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACtt                                                                           >  2:1069071/1‑143 (MQ=255)
                                                              taGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTCGTCTGTTTCTTACtt                                                                           >  2:198181/1‑143 (MQ=255)
                                                                       agCATGATACGGCTTACCTGTGCTTATACATAAAAGCAACAGAATGGTAACATGTTATCGCGGGTAAGTCAATTGAGCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGa                                                                  <  1:730804/142‑1 (MQ=255)
                                                                           tGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGtt                                                              >  1:775176/1‑143 (MQ=255)
                                                                                    ttACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTTGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAt                                                     <  1:520446/143‑1 (MQ=255)
                                                                                    ttACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGAGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAt                                                     >  2:568085/1‑143 (MQ=255)
                                                                                        cTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTGTTGGAttttt                                                 >  1:933197/1‑143 (MQ=255)
                                                                                        cTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTCTTGGCgtttt                                                 >  1:11997/1‑143 (MQ=255)
                                                                                        cTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGGttttt                                                 >  1:683145/1‑143 (MQ=255)
                                                                                        cTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAttttt                                                 >  1:742161/1‑143 (MQ=255)
                                                                                          gTTCTTCTACATAAACGCAACTGAATTGTAACATTTTATCGCGGGTAAGTCAATTGATCCACGTCATTTGTGTTGCCATGTCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTAt                                               <  2:72406/143‑1 (MQ=255)
                                                                                                ataCATAAAAGCAACAGAATGGTAACATGTTATCGCGTGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATcgacga                                         <  2:717356/143‑1 (MQ=255)
                                                                                                    aTAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATAct                                     <  1:630402/143‑1 (MQ=255)
                                                                                                       aaaGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATc                                              <  1:734757/131‑1 (MQ=255)
                                                                                                       aaaGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATc                                              >  2:734757/1‑131 (MQ=255)
                                                                                                            aaCAGAGTGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGCCCGGCCTTTGCTTTGTCTGTTTTTTACTTTTGCCCTGACGTTTTATTGGCTTTTTATCGCCGATACTCTCCGttt                             >  1:421984/1‑143 (MQ=255)
                                                                                                              tctaaTTGTCACATTTTATCGCGCGTAACTCAATTGATCCCCGTCATTTAACTTGCTATATCCTGAGCGGCTTTTTCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTaa                           <  2:943486/140‑1 (MQ=255)
                                                                                                                    ggcaaCCTTCTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCa                     <  1:830815/140‑1 (MQ=255)
                                                                                                                     gTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAg                    <  1:986119/143‑1 (MQ=255)
                                                                                                                      tAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAgg                   >  1:838700/1‑143 (MQ=255)
                                                                                                                       aaCATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGCCGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGt                  >  2:786813/1‑143 (MQ=255)
                                                                                                                        aCATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATAGACGATACTCTCCGTTTAAGCGGCAGGtt                 >  1:691686/1‑143 (MQ=255)
                                                                                                                          aTTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTc               >  1:1042689/1‑143 (MQ=255)
                                                                                                                                tCGTGGGTAAGTCAATTGACCCCCGTCATTAATCTGGCTATATCCTGAGCGGCCTTTGCTTTTTCTGTTTCTTACTTTTGCCCTGACGTGTTATTGGATTTTTATCGACGATATTCTCCGTTTAAGCGGCAGGTTTCCGCTGt         <  1:403566/143‑1 (MQ=255)
                                                                                                                                  gcgTGTAAGTCAATTGATCCCCGTTATTTATCTGGCTATATCCTGAGCGGCCTGTGCGCTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTAc       <  1:43208/143‑1 (MQ=255)
                                                                                                                                  gcgGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCGCTGGCGTTTTTTTGGGTTTTTATCGTCGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGCAc       >  1:846052/1‑143 (MQ=255)
                                                                                                                                  gcgGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTCCCCTGTCGTTTTGTTGGCTTTTTATCGACGTTCCTCTCCGTTTAAGCGACAGGTCTCAGCTGTAc       >  1:288088/1‑143 (MQ=255)
                                                                                                                                   cgGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGCTTCTTATCGCCGATTCTCTCCGTTTCCGCGGCAGTTTTTCGCTGTACg      >  1:525484/1‑143 (MQ=255)
                                                                                                                                     ggTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTa    >  1:1082031/1‑143 (MQ=255)
                                                                                                                                       tAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaaa  <  1:49052/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAA  >  NC_000913/4160715‑4160992

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: ‑