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→AAT) frdBfumarate 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
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
acaGTTGCTGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCTTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTACTCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACcat                                                                                                                                                                                                                              >  2:196242/1‑251 (MQ=255)
      gATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGc                                                                                                                                                                                                                         <  2:124375/250‑1 (MQ=255)
                     cTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAAt                                                                                                                                                                                                         >  1:44069/1‑251 (MQ=255)
                           aTACTTCGCCATCTGCGCCGGGGCCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCCCGTTGTTAACCATCATGCCGCAGGAACCCCAAACCGcccc                                                                                                                                                                                                   >  2:207293/1‑249 (MQ=255)
                                         gcgcCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGTGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGTTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGAcc                                                                                                                                                                                        <  1:216886/248‑1 (MQ=255)
                                                  tCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGCTGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTAATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTAATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTc                                                                                                                                                                             >  1:227400/1‑250 (MQ=255)
                                                          ttAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTACGGt                                                                                                                                                                    >  1:97712/1‑251 (MQ=255)
                                                                        cgcgGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGCCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCttt                                                                                                                                                      >  2:161351/1‑251 (MQ=255)
                                                                          cgGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTg                                                                                                                                                     >  2:107045/1‑250 (MQ=255)
                                                                           ggTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGTCTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATTCCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCTCGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGAt                                                                                                                                                   <  1:42365/251‑1 (MQ=255)
                                                                                                       ttGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAACCCAGCGCATCCAGTAATGAGGt                                                                                                                        >  1:77990/1‑250 (MQ=255)
                                                                                                            cGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCTGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGGCAGGTTGTCTTTGATGTAGCCCAGCGCTTCCAGTAATGTGGTAGTTGc                                                                                                                  <  2:100899/251‑1 (MQ=255)
                                                                                                                tCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGt                                                                                                                <  1:102565/249‑1 (MQ=255)
                                                                                                                 ccAGTCTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTTTCTTTTATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCAt                                                                                                             <  1:213632/251‑1 (MQ=255)
                                                                                                                        ttCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGc                                                                                                        <  2:16848/249‑1 (MQ=255)
                                                                                                                          cGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTATCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGGAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACtt                                                                                                    >  2:156190/1‑251 (MQ=255)
                                                                                                                                aGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGa                                                                                              <  1:107045/251‑1 (MQ=255)
                                                                                                                                                          cgTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTTTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACt                                                                    <  1:99809/251‑1 (MQ=255)
                                                                                                                                                              cAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTcc                                                                 <  1:123046/250‑1 (MQ=255)
                                                                                                                                                                                   cTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCa                                            >  1:90485/1‑250 (MQ=255)
                                                                                                                                                                                       aaCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAAtttt                                       >  2:126062/1‑251 (MQ=255)
                                                                                                                                                                                        aCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTc                                      <  1:97403/251‑1 (MQ=255)
                                                                                                                                                                                            tCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAgg                                   <  1:30442/250‑1 (MQ=255)
                                                                                                                                                                                               tACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     <  1:87733/117‑1 (MQ=255)
                                                                                                                                                                                               tACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     >  2:87733/1‑117 (MQ=255)
                                                                                                                                                                                                           aatcaaaaacaaaGGATTTACATGCCCATTTAG‑CACGTTGTTAACCATCAACCCGCAGGCACCACAAATCGCCATACGGCAGGCCCAGCGGTACCTCACGTCCGGAGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGCCATAAGGCACTTCCTAGCATGCGCTATGCGGTGCGGCATCGCCTTCCGGGTTATAGCCCACCACCTCAATTTTCAGGTTTTTCATCTCAGcc                    >  2:205528/12‑250 (MQ=255)
                                                                                                                                                                                                                             tACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCcttctt  <  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: ‑