Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,397,191 A→G 80.0% F174S (TTC→TCC)  nuoM ← NADH:ubiquinone oxidoreductase, membrane subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,397,1910AG80.0% 100.1 / 17.9 45F174S (TTC→TCC) nuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base G (21/15);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAG  >  W3110S.gb/2397129‑2397255
                                                              |                                                                
cAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTg                                                             <  1:2180530/68‑1 (MQ=255)
cAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTg                                                             <  1:2076227/68‑1 (MQ=255)
 aaaaaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTg                                                             <  1:1493729/67‑1 (MQ=255)
  aaaaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGTGTGTAAATGAAGGACTTGGtt                                                          <  1:1316724/69‑1 (MQ=255)
  aaaaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGtt                                                          <  1:656216/69‑1 (MQ=255)
  aaaaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGtt                                                          <  1:263865/69‑1 (MQ=255)
  aaaaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGtt                                                          <  1:1151422/69‑1 (MQ=255)
       cAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTgccgc                                                     <  1:2020872/69‑1 (MQ=255)
          cgcCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGt                                                  <  1:1036954/69‑1 (MQ=255)
                      aaTCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGtttt                                      <  1:2212485/69‑1 (MQ=255)
                       aTCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTa                                     <  1:3006632/69‑1 (MQ=255)
                         cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGtttt                                      >  1:444942/1‑66 (MQ=255)
                         cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGttt                                       >  1:3504842/1‑65 (MQ=255)
                         cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc                                   >  1:376512/1‑69 (MQ=255)
                         cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc                                   >  1:148004/1‑69 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTgccg                                                      <  1:374039/49‑1 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTgccg                                                      <  1:167426/49‑1 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTgccg                                                      <  1:1300409/49‑1 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATa                                              >  1:1633587/1‑57 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATa                                              >  1:1971838/1‑57 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTa                                     >  1:600197/1‑66 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTa                                     >  1:109545/1‑66 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTAcc                                   >  1:3095826/1‑68 (MQ=255)
                          aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACATGGTTGCCGCCGTGATACGCGTTTTAcc                                   >  1:2393661/1‑68 (MQ=255)
                            cATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGccgccg                                                   >  1:2927847/1‑50 (MQ=255)
                            cATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGt                                                  >  1:3533963/1‑51 (MQ=255)
                            cATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGt                                                  >  1:3052179/1‑51 (MQ=255)
                            cATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGt                                                  >  1:196906/1‑51 (MQ=255)
                                  cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg                           <  1:128760/68‑1 (MQ=255)
                                    ggCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        <  1:1453862/69‑1 (MQ=255)
                                    ggCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        <  1:1486186/69‑1 (MQ=255)
                                    ggCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        <  1:2337171/69‑1 (MQ=255)
                                    ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        <  1:531764/69‑1 (MQ=255)
                                     gCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        >  1:2817813/1‑68 (MQ=255)
                                      ccACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGttt                                       <  1:1193441/52‑1 (MQ=255)
                                       cACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTAt                     >  1:210833/1‑69 (MQ=255)
                                        aCTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGttt                                       >  1:1308507/1‑50 (MQ=255)
                                           cGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTa                                     <  1:2701247/49‑1 (MQ=255)
                                                       aaTGAAGGACTTGGTTGCCGCCGTGCTACGCGTTTTACCGTCa                               >  1:562235/1‑43 (MQ=38)
                                                       aaTGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCttt                       >  1:3658118/1‑51 (MQ=255)
                                                       aaTGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        >  1:1304449/1‑50 (MQ=255)
                                                       aaTGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        >  1:1840797/1‑50 (MQ=255)
                                                       aaTGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                        >  1:1158463/1‑50 (MQ=255)
                                                       aaTGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCg       >  1:2558731/1‑67 (MQ=255)
                                                       aaTGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCg       >  1:2166051/1‑67 (MQ=255)
                                                          gaaggaCTTGGTTGCCGCCGTGATACGCGTTTTACCGTCATAGGCTTTATGCCCCCACAGTGCGATCAg  <  1:2663127/69‑1 (MQ=255)
                                                              |                                                                
CAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAG  >  W3110S.gb/2397129‑2397255

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