Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,844 C→G 20.7% A642A (GCC→GCG rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,8440CG20.7% 27.7 / 7.2 33A642A (GCC→GCGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (20/6);  new base G (4/3);  total (24/9)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTG  >  NC_000913/1708700‑1708992
                                                                                                                                                |                                                                                                                                                    
tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTAGACGCCGCTATTGCCCGCGCCCACGCACAAAACCACGAACAGCAGAAGGCTAAAAACGAGcc                                                                                                                                                 >  2:904526/1‑150 (MQ=14)
   ttaccGTGCCAAATCGCGCAAGCTGGAACAGCAACAGCCTACTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                                                                                                                                              <  1:1132256/147‑1 (MQ=11)
         ggCCACAGGGCGCAAGCTGGAACAGCATCAGGCTAATGCGGAACCAGAACCCCTGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                                                                                                                                        <  2:1941206/149‑1 (MQ=11)
             acatctcgcAAGCTGGAAGAGCAACAGCCTAATTCGGAACCAGAAGAACATGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                                                                                                                                    <  1:620786/144‑1 (MQ=14)
             aaaTCGCTCAAGATGGTACAGTAACAGTCTTTTGCTGAACCAGCATATCAGGTTGATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                                                                                                                                    <  2:1541432/150‑1 (MQ=11)
                 cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGCGCACACGCCCGCAAGCGGGCCCCGCACCAGGGTAATGCCGAGCCAGAACTCCAGGTCGATc                                                                                                                                >  2:1823845/1‑150 (MQ=14)
                                 gCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCTCAATTGCCCGTGCCAAAGCACACAACCTGGAACAGCACCAGGCTAAAGCCGAGCCAGAACATCAGGTACATAAGCGCATAGcggccgg                                                                                                                >  1:635439/1‑149 (MQ=14)
                                                ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGGACCGCAACAGCCAACGGCGGAGCCCGAGGCACGGGGCGATCCGCGCAAAGCTGACGTCGAAGCAGCTAATGc                                                                                                 >  1:1054732/1‑150 (MQ=11)
                                                           aaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTGGATCCGCGCAAAGCCGACGTCGAACCCGCTAGTGCGCGAGCGAAAg                                                                                      >  1:1901651/1‑150 (MQ=11)
                                                                   gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACCCAAGCTCGACCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTATAACCGCGCAAATCCGCAGTTTCCGCCGCTATTGACCGCGCCAACGCGCGaaaa                                                                              >  2:1645052/1‑150 (MQ=14)
                                                                        gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGACCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGCTACGCGCAAAGCGGCCGTCGCAGCGGCTATTGCCCGCGCCACAGCCCGAAAGCGCGa                                                                         >  1:1329427/1‑150 (MQ=11)
                                                                                             cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTGGATCAGCGCAAAGCCGCCGTCGAAGCGGCTTTTGCGCGAGCCAAAGCGCTCAACCGGGAACAGCAAAAGGCTAATGCggg                                                    >  1:1791836/1‑149 (MQ=14)
                                                                                              gCTATTGCCCGTGCCAAAGCACGCAATCTGGAACAGCATCAGCCTAATCCCGACCCAGAACAACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg                                                   <  1:663115/150‑1 (MQ=17)
                                                                                                   tGCCCGTGCCACAGCACGCAATCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGCTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                                              <  1:1518968/150‑1 (MQ=255)
                                                                                                    gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGTGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTCGACGCGGCTATTGCCCGAGCACAAGCGAGCATACGGGCACAGCAGACGGGTTTTGCGGAGCCagaa                                             >  2:721071/1‑150 (MQ=17)
                                                                                                    gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCCAACGCGCGCAAACGGGACCAGCACCCGGCTCATGCCGAGCCcgaa                                             >  1:1489903/1‑150 (MQ=14)
                                                                                                     cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTCGACGCCGCTATTGCCCGAGCGCAGGCGCGCAAACCGGAACAGCACCCGGCTTATGCGGAGCCcgaag                                            >  1:273328/1‑150 (MQ=17)
                                                                                                                     cAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGAAAACGGGAACAGCAACCGGATAATGCGGAGCCAGAAGAACAGGTAGAACcgcc                            >  1:1775355/1‑149 (MQ=14)
                                                                                                                       aGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAACGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGATCAGGTTGATCCGCGCa                          >  1:670743/1‑150 (MQ=255)
                                                                                                                         cTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCCAAGCGCGCAAACGGGAACAGCAGCAGGCTAATGCCGAGCCAGAAGAGCAGGTGGCTCCCCGCaaa                        >  1:1944500/1‑150 (MQ=14)
                                                                                                                          tGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCAAAAGCCAAAAAACGCGAACAGCACCAGGATAAACCCGAGCCCTCATAACAGGATCATCCTCTCAAAg                       >  2:909527/1‑150 (MQ=2)
                                                                                                                             aaCAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                       <  2:138230/99‑1 (MQ=255)
                                                                                                                             aaCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                       >  1:138230/1‑99 (MQ=16)
                                                                                                                             aaCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCCAAAGCGCGAAAACGCGCACAGCACCAGGATAATGCCGAGCAAGAAGAGCAGTTTGATACGCGGAAATCTg                    >  2:1816360/1‑150 (MQ=14)
                                                                                                                              aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCAAAGCCAGCAAACGGGAACAGCAACCGGCTAATACGGAGCCAGAAAAACAGGTAAATCAGCGAAACGCAGc                   >  1:1360393/1‑150 (MQ=14)
                                                                                                                               cAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGCCGCCGCTATTGCCCGCGCCAAAGCCAAAAAACGGGACCAGCAACCGGATGTTGCGGAGCCAGAAGACCATGTTGAACCGCGCAACGcgccc                  >  1:1297887/1‑146 (MQ=14)
                                                                                                                                 gCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCAAAAGCGCGCAAACAGGAACAGCAACCGGCAAATGCGGAGCCAGAAGAACAGCTTGCTCCGCGCAAAGCTGCCGt                >  1:1502335/1‑150 (MQ=14)
                                                                                                                                 gCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCAGTAGACGCCGCTATTGCCCGAGCCAAAGCCCGCAAACGCGAACAGCAACAGGGAAAGGCGGAGCCAGAAGATCAGTTTGATACGAGCAAAGcgtcagg                >  2:1873980/1‑144 (MQ=14)
                                                                                                                                    aCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCCCGCAAACGCGAACAGCAACAGGCTAATGCGGTGCCAGAAGAACAGGTTGCTCCGCGCAACGCTGCAGTCgc             >  1:700751/1‑149 (MQ=14)
                                                                                                                                       ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                       >  1:794144/1‑89 (MQ=255)
                                                                                                                                       ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                       <  2:794143/89‑1 (MQ=255)
                                                                                                                                       ggCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCTCAAACGGGAACAGCAACAGGATAATCCGGAGCAAGAAGAAAAGGTTGATCCGCGAAAAGCTGCAGtagccgc          >  2:1008035/1‑144 (MQ=14)
                                                                                                                                          aaaTGCGGAAGCAGCACACCAGGCCGATCCGCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAATCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc       <  2:1817469/149‑1 (MQ=11)
                                                                                                                                               ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTg  >  1:1145999/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTG  >  NC_000913/1708700‑1708992

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