Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,381,782 G→A G179G (GGC→GGT frdA ← fumarate reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,7820GA97.3% 116.5 / NA 37G179G (GGC→GGTfrdAfumarate reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/22);  minor base C (1/0);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGTGTGCAGCATATGG  >  NC_000913/4381653‑4381917
                                                                                                                                  |                                                                                                                                       
aaCGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCcatca                                                                                                                               >  1:451092/1‑140 (MQ=255)
       tGGTTTAGCGCCATCCACATCCCGTCACCGGCAGCCATGCCGCCGTTGTTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCAt                                                                                                                        <  2:352082/140‑1 (MQ=255)
             aGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTAc                                                                                                                  >  2:131320/1‑140 (MQ=255)
                   atacccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTCCCTACCTTCATGTCCATTGCTACCAGGcc                                                                                                            >  1:264339/1‑140 (MQ=255)
                    tacccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTGATTGCTACCAGGCcg                                                                                                           >  1:323309/1‑140 (MQ=255)
                         accccGTCCCCGGAAACGAGGCCGCCTTGGGTGTTGAAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAAc                                                                                                      <  2:455406/137‑1 (MQ=255)
                                 aCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcat                                                                                              >  2:825920/1‑140 (MQ=255)
                                 aCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcat                                                                                              >  1:767912/1‑140 (MQ=255)
                                     gTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCAAACATGACCatcatc                                                                                          >  1:277080/1‑140 (MQ=255)
                                      tAACGATGCCGCCGTTGGTGTTGTAACGATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCAtcatca                                                                                         <  1:166229/140‑1 (MQ=255)
                                       aaCGATGCGGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCaa                                                                                        <  1:446770/140‑1 (MQ=255)
                                           aggccgccGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAg                                                                                    <  2:708946/138‑1 (MQ=255)
                                           aTGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAg                                                                                    <  2:571796/140‑1 (MQ=255)
                                                 ccgTTGGTGTTGAAACATAAACCGCGCCCCGACCCCCCATTAGCCATATCGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATc                                                                              <  2:277080/140‑1 (MQ=255)
                                                  gtttGGTGTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATcc                                                                             <  2:441044/138‑1 (MQ=255)
                                                    ttGGTTTTGTAACGAAACCCGCCCCCCGCCCCGCCAGTAGCCATAACGCCCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAg                                                                           <  2:323309/140‑1 (MQ=255)
                                                                    aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCa                                                           >  1:566600/1‑140 (MQ=255)
                                                                            aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCt                                                    >  2:29424/1‑140 (MQ=255)
                                                                                 cACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATc                                               >  1:657558/1‑140 (MQ=255)
                                                                                          tAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAAt                                      >  1:68209/1‑140 (MQ=255)
                                                                                             ccATAACGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGc                                   >  1:464732/1‑140 (MQ=255)
                                                                                                 aaCGACCGCGTTAGCTCGGATCAGCGCCAACGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCagag                               <  1:343231/140‑1 (MQ=255)
                                                                                                 aaCGACCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCagag                               <  1:16752/140‑1 (MQ=255)
                                                                                                     aCCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGt                           <  1:825920/140‑1 (MQ=255)
                                                                                                     aCCGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGt                           <  2:657558/140‑1 (MQ=255)
                                                                                                      ccGCGTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAAAGCTAGAGCTGCGGGAATTGCAGAGAGGTc                          >  1:236871/1‑140 (MQ=255)
                                                                                                          gTTAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGa                      <  2:767912/140‑1 (MQ=255)
                                                                                                           ttAGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAGAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGaa                     >  1:522644/1‑140 (MQ=255)
                                                                                                             aGCACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACa                   <  1:131320/140‑1 (MQ=255)
                                                                                                               cACGGATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGc                 <  1:774170/140‑1 (MQ=255)
                                                                                                                  ggATCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCgtg              <  2:566600/140‑1 (MQ=255)
                                                                                                                    aTCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGtgtg            <  2:642929/140‑1 (MQ=255)
                                                                                                                    aTCTGCACCAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGtgtg            <  2:68209/140‑1 (MQ=255)
                                                                                                                    aTCTGCACCAGCGTACCTTCCACCATGTTCATTGCCACCCGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGtgtg            <  2:264339/140‑1 (MQ=255)
                                                                                                                    aTCTGCACCAGCGTACCCTCCACCATGTTCATTGCCACCCGGCCGCGAGCATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGtgtg            <  2:226655/140‑1 (MQ=255)
                                                                                                                            cAGCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGTGTGCAGCatat    <  1:29424/140‑1 (MQ=255)
                                                                                                                              gCGTACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGTGTGCAGCATATgg  >  1:170358/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                       
AACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGGAATTGCAGAGAGGTCTGGAACAGCGTGTGCAGCATATGG  >  NC_000913/4381653‑4381917

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: