Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,380,292 2 bp→A coding (225‑226/735 nt) frdB ← fumarate reductase iron‑sulfur protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,380,2920G.96.3% 117.4 / ‑3.2 27coding (226/735 nt)frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (12/14);  minor base A (1/0);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_0009134,380,2930TA96.3% 88.9 / ‑3.2 27K75N (AAA→AATfrdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (12/14);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTT  >  NC_000913/4380062‑4380533
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
acaGTTGCTGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCTTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTACTCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACcat                                                                                                                                                                                                                              >  2:196242/1‑251 (MQ=255)
      gATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGc                                                                                                                                                                                                                         <  2:124375/250‑1 (MQ=255)
                     cTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAAt                                                                                                                                                                                                         >  1:44069/1‑251 (MQ=255)
                           aTACTTCGCCATCTGCGCCGGGGCCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCCCGTTGTTAACCATCATGCCGCAGGAACCCCAAACCGcccc                                                                                                                                                                                                   >  2:207293/1‑249 (MQ=255)
                                         gcgcCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGTGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGTTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGAcc                                                                                                                                                                                        <  1:216886/248‑1 (MQ=255)
                                                  tCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGCTGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTAATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTAATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTc                                                                                                                                                                             >  1:227400/1‑250 (MQ=255)
                                                          ttAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTACGGt                                                                                                                                                                    >  1:97712/1‑251 (MQ=255)
                                                                        cgcgGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGCCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCttt                                                                                                                                                      >  2:161351/1‑251 (MQ=255)
                                                                          cgGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTg                                                                                                                                                     >  2:107045/1‑250 (MQ=255)
                                                                           ggTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGTCTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATTCCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCTCGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGAt                                                                                                                                                   <  1:42365/251‑1 (MQ=255)
                                                                                                       ttGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAACCCAGCGCATCCAGTAATGAGGt                                                                                                                        >  1:77990/1‑250 (MQ=255)
                                                                                                            cGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCTGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGGCAGGTTGTCTTTGATGTAGCCCAGCGCTTCCAGTAATGTGGTAGTTGc                                                                                                                  <  2:100899/251‑1 (MQ=255)
                                                                                                                tCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGt                                                                                                                <  1:102565/249‑1 (MQ=255)
                                                                                                                 ccAGTCTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTTTCTTTTATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCAt                                                                                                             <  1:213632/251‑1 (MQ=255)
                                                                                                                        ttCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGc                                                                                                        <  2:16848/249‑1 (MQ=255)
                                                                                                                          cGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTATCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGGAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACtt                                                                                                    >  2:156190/1‑251 (MQ=255)
                                                                                                                                aGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGa                                                                                              <  1:107045/251‑1 (MQ=255)
                                                                                                                                                          cgTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTTTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACt                                                                    <  1:99809/251‑1 (MQ=255)
                                                                                                                                                              cAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTcc                                                                 <  1:123046/250‑1 (MQ=255)
                                                                                                                                                                                   cTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCa                                            >  1:90485/1‑250 (MQ=255)
                                                                                                                                                                                       aaCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAAtttt                                       >  2:126062/1‑251 (MQ=255)
                                                                                                                                                                                        aCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTc                                      <  1:97403/251‑1 (MQ=255)
                                                                                                                                                                                            tCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAgg                                   <  1:30442/250‑1 (MQ=255)
                                                                                                                                                                                               tACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     <  1:87733/117‑1 (MQ=255)
                                                                                                                                                                                               tACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     >  2:87733/1‑117 (MQ=255)
                                                                                                                                                                                                           aatcaaaaacaaaGGATTTACATGCCCATTTAGCACGTTGTTAACCATCAACCCGCAGGCACCACAAATCGCCATACGGCAGGCCCAGCGGTACCTCACGTCCGGAGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGCCATAAGGCACTTCCTAGCATGCGCTATGCGGTGCGGCATCGCCTTCCGGGTTATAGCCCACCACCTCAATTTTCAGGTTTTTCATCTCAGcc                    >  2:205528/12‑250 (MQ=255)
                                                                                                                                                                                                                             tACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCcttctt  <  1:136326/250‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTT  >  NC_000913/4380062‑4380533

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