Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,397,191 A→G 55.1% 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,1910AG55.1% 15.2 / 63.9 49F174S (TTC→TCC) nuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base A (10/12);  new base G (16/11);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTG  >  W3110S.gb/2397124‑2397248
                                                                   |                                                         
gTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACt                                                        >  1:998876/1‑71 (MQ=255)
gTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACt                                                        >  1:1901536/1‑71 (MQ=255)
gTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACt                                                        >  1:148258/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:1057865/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:818920/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:731948/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:851721/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:1886693/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:1298922/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:1446638/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTgg                                                     >  1:137392/1‑71 (MQ=255)
   aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTg                                                      >  1:667011/1‑70 (MQ=255)
             aGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGt                                           <  1:1160631/71‑1 (MQ=255)
             aGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGt                                           <  1:1454528/71‑1 (MQ=255)
                       tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATATGCGtt                                 >  1:1550071/1‑71 (MQ=255)
                       tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGtt                                 >  1:1733607/1‑71 (MQ=255)
                       tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGtt                                 >  1:1582083/1‑71 (MQ=255)
                       tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGtt                                 >  1:1551756/1‑71 (MQ=255)
                                     aCCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg                    >  1:1347354/1‑70 (MQ=255)
                                     aCCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGc                   >  1:53851/1‑71 (MQ=255)
                                     aCCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGc                   >  1:360614/1‑71 (MQ=255)
                                     aCCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGc                   >  1:281513/1‑71 (MQ=255)
                                     aCCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGc                   >  1:1945691/1‑71 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1394857/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1231019/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1454936/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1108402/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1131254/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:984262/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1418752/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:783122/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1872206/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1347664/71‑1 (MQ=255)
                                       cAGGCCACTCGCCTGGGGGTAAATGAAGAACTTGGTTGCCGCCGt                                           <  1:1067686/45‑1 (MQ=255)
                                          gCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTAt              <  1:875035/71‑1 (MQ=255)
                                          gCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTAt              <  1:564587/71‑1 (MQ=255)
                                          gCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTAt              <  1:370001/71‑1 (MQ=255)
                                              cTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCTGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:850430/64‑1 (MQ=255)
                                              cTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1078145/64‑1 (MQ=255)
                                              cTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1787425/64‑1 (MQ=255)
                                                 gCCTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGcccc         <  1:1971653/69‑1 (MQ=255)
                                                  ccTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCcac      <  1:1207333/71‑1 (MQ=255)
                                                   cTGGGTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCcac      <  1:1102828/70‑1 (MQ=255)
                                                     gggTGTAAATGAAGGACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                 <  1:1830876/57‑1 (MQ=255)
                                                      ggTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTg  >  1:1552273/1‑71 (MQ=255)
                                                      ggTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTg  >  1:632580/1‑71 (MQ=255)
                                                      ggTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTg  >  1:1031516/1‑71 (MQ=255)
                                                      ggTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTg  >  1:1420998/1‑71 (MQ=255)
                                                      ggTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTg  >  1:1038472/1‑71 (MQ=255)
                                                                   |                                                         
GTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTG  >  W3110S.gb/2397124‑2397248

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