Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,751 A→G 20.4% E611E (GAA→GAG 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,7510AG20.4% 48.1 / 7.3 34E611E (GAA→GAGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (7/20);  new base G (5/2);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 6.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAT  >  NC_000913/1708603‑1708894
                                                                                                                                                    |                                                                                                                                               
caaTCTCCCGCGCCAAACCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                                                                <  1:1559656/148‑1 (MQ=255)
 tATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAcc                                                                                                                                               <  1:1728567/150‑1 (MQ=255)
   tCGACCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCag                                                                                                                                             <  2:1947099/150‑1 (MQ=255)
      cccGCGCCCAACCCCGCAACCAGGAACACCAACAGGCAAATCCGGAACCAGAACCCCAGGCCGATCCCCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                                                                                                                                          <  2:596522/150‑1 (MQ=14)
         gcgactaaGCGCGCAATCAGGAACATCATCAGGATAATCAGCAACCAGCACAACAGGTCGATCCTCGTAAACCTGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAc                                                                                                                                       <  1:1059263/144‑1 (MQ=14)
          ctccAACGCGCGCAAGCTGGAACAGCATCAGGCTAATGAGGAACCATATCAAGATGTCGATCCGCTAAAAGCCTCCGTCGAAGCAGCTATAGCACTTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                                                                      <  1:676839/148‑1 (MQ=14)
          cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                                                                      <  1:1028296/150‑1 (MQ=255)
          cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                                                                      <  1:1806748/150‑1 (MQ=255)
            ccaactcgcgcAAGCTGGAACACCAACAGGATAATCAGGAACCAGAACCACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                                                                    <  1:161097/144‑1 (MQ=14)
                       atgCTGTAAGAGCTACAGGCTAATCAGGAAGCAGCACCAGAGGTGGATCCTCGCAAAGCCGCCGTGGAAGCAGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCAATCCgcgc                                                                                                                         <  2:1774532/148‑1 (MQ=14)
                              aagatcAACAGGCTAATTCGGAACCAGAACAAGAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                                                                  <  2:1950816/145‑1 (MQ=255)
                              aagatcAACAGGATACTCCCGAACCAGAACACGAGGTCGATCCTCGTAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                                                                  <  2:503254/145‑1 (MQ=14)
                              aaCACCAACAGGATAATTAGGAACCAGCACAAGTGGTCGATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                                                                  <  1:1166187/150‑1 (MQ=14)
                                              tGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                                                         >  2:2072134/1‑111 (MQ=255)
                                                          acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGCAGCCGCTATTGCCCGAGCCaa                                                                                      >  1:1969601/1‑150 (MQ=255)
                                                           caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCGAAGCCGCGCGCGAAGACTCTTTTGCGCGTggaaag                                                                                     >  2:1232782/1‑145 (MQ=11)
                                                             agaGGTGGATCCTCGCAAAGCGGCAGTCGATGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAGAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGc                                                                                   <  1:701101/148‑1 (MQ=14)
                                                              caGGTTGAGCCTCGCAAATCCGCAGTAGAAGCCGCTATTGCCCTTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCa                                                                                  <  1:1026686/150‑1 (MQ=14)
                                                                                 ccgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga                                                                                                                                            >  1:2065235/1‑73 (MQ=21)
                                                                                  cgccgcCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGCACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACa                                                              <  1:708174/150‑1 (MQ=255)
                                                                                           aGCCGCCATTGCCCCGGCCCAACCGCGCAAGCTGGAACACCAACAGGCAAATGCGGAACCAGCACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGc                                                     <  1:20133/150‑1 (MQ=11)
                                                                                                            ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAGGCCGCCGTCGACGCCGCTTTTGACCGTGCCCAAGCACGCAAGCTGTCACAGCAACAGGGTAATGCAGAGCAAGAACAGCAgg                                    >  1:1641112/1‑150 (MQ=14)
                                                                                                             cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCgcgc                                                                                                                         <  1:1372695/64‑1 (MQ=255)
                                                                                                             cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCgcgc                                                                                                                         >  2:1372695/1‑64 (MQ=25)
                                                                                                                  cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                                                                                     <  2:2087703/79‑1 (MQ=11)
                                                                                                                  cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCGGc                                                                                                     >  1:2087703/1‑79 (MQ=11)
                                                                                                                      cAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg                          >  2:2012936/1‑150 (MQ=255)
                                                                                                                                 aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg               <  1:2012936/150‑1 (MQ=255)
                                                                                                                                    aaCGGGCTAATGCGTAAGCACAAGAACTTTTCTATCTTAGCATAGCCGACGTTGAAGCCTCTTTTGCGCGCGCCATTGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg            <  2:1924067/150‑1 (MQ=2)
                                                                                                                                           tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGcc                                                                                        <  2:1961450/67‑1 (MQ=25)
                                                                                                                                           tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGcc                                                                                        >  1:1961450/1‑67 (MQ=25)
                                                                                                                                             aTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCCAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa   >  1:751097/1‑150 (MQ=255)
                                                                                                                                              tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGACGAGCCAGCACACCAGGTCCATCCCCGCAACGCAGCCGTCGCAGACGCTAt  >  1:1651674/1‑150 (MQ=14)
                                                                                                                                              tGCGGAACCAGAAGAAAAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGTGCCAACGCACCCAAGCTGGAACAGCAGCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt  >  1:1990412/1‑150 (MQ=14)
                                                                                                                                                    |                                                                                                                                               
CTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAT  >  NC_000913/1708603‑1708894

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