Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,168,260 G→A 12.2% G200D (GGT→GAT)  ndh → respiratory NADH dehydrogenase 2/cupric reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,168,2600GA12.2% 137.2 / 8.5 49G200D (GGT→GAT) ndhrespiratory NADH dehydrogenase 2/cupric reductase
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base A (4/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGCGACGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTATTTTGCC  >  W3110S.gb/1168194‑1168326
                                                                  |                                                                  
gcggcgCGACGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTAc                                                                <  1:2558807/71‑1 (MQ=255)
 cggcgCGACGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACa                                                               <  1:2798749/71‑1 (MQ=255)
 cggcgCGACGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACa                                                               <  1:23045/71‑1 (MQ=255)
        aCGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACaaa                                                             >  1:2314145/1‑66 (MQ=255)
        aCGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACaaa                                                             >  1:2386381/1‑66 (MQ=255)
        aCGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACa                                                               <  1:3941255/64‑1 (MQ=255)
        aCGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACa                                                               <  1:374204/64‑1 (MQ=255)
         cGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTg                                                       <  1:123958/71‑1 (MQ=255)
              gtAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCaa                                                  >  1:794677/1‑71 (MQ=255)
               tAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAAc                                                 >  1:3546985/1‑71 (MQ=255)
               tAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGATTACAAAGGCCTGACCAAc                                                 >  1:709604/1‑71 (MQ=255)
               tAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGATTACAAAGGCCTGACCAAc                                                 >  1:2535646/1‑71 (MQ=255)
               tAGAACTCTCCGATGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAAc                                                 >  1:3470143/1‑71 (MQ=255)
                      ctcCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACg                                                >  1:3789335/1‑65 (MQ=255)
                      ctcCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGccc                                          >  1:3268025/1‑71 (MQ=255)
                       tcCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCaa                                                  >  1:688697/1‑62 (MQ=255)
                       tcCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCaa                                                  >  1:1982776/1‑62 (MQ=255)
                       tcCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCaa                                                  >  1:1546599/1‑62 (MQ=255)
                              aaTTGCACAACGCGGTCAAGCAACTGCACAGCTACGATTACaaa                                                             <  1:1885356/44‑1 (MQ=255)
                                  gCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGc                              >  1:3304352/1‑71 (MQ=255)
                                    acaACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCt                             <  1:1618385/70‑1 (MQ=255)
                                    acaACGCGGTCAAGCAACTGCACAGCTACGATTACAAAGGCCTGACCAACGAAGCCCTGAACGt                                   <  1:715743/64‑1 (MQ=255)
                                     caACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTgg                           <  1:615158/71‑1 (MQ=255)
                                      aaCGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGt                                   <  1:1833744/62‑1 (MQ=255)
                                      aaCGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAACCCCTGa                                       <  1:632985/58‑1 (MQ=255)
                                           ggTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAg                     <  1:2551721/71‑1 (MQ=255)
                                            gTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGc                                            <  1:664213/47‑1 (MQ=255)
                                            gTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGc                              <  1:4245221/61‑1 (MQ=255)
                                             tCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGa                       <  1:941664/67‑1 (MQ=255)
                                             tCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGa                       <  1:668785/67‑1 (MQ=255)
                                              cAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTg                            >  1:2120601/1‑61 (MQ=255)
                                                 gCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTg                                        <  1:958775/46‑1 (MQ=255)
                                                 gCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTg                                        <  1:964078/46‑1 (MQ=255)
                                                 gCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTg                                        <  1:3895151/46‑1 (MQ=255)
                                                   aaCTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCgg                  >  1:3437397/1‑66 (MQ=255)
                                                     cTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACg           <  1:3398125/71‑1 (MQ=255)
                                                      tGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACg           >  1:2916793/1‑70 (MQ=255)
                                                       gCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTa         >  1:3050278/1‑71 (MQ=255)
                                                       gCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTa         >  1:2067012/1‑71 (MQ=255)
                                                        cacaGCTACGGTTACAAAGGCCTGACCAACGACCCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTAt        >  1:475975/1‑71 (MQ=255)
                                                        cacaGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGaa                                      >  1:2394186/1‑41 (MQ=255)
                                                        cacaGCTACGATTACAAAGGCCTGACCAACGAAGCCCTGAAc                                     >  1:2917233/1‑42 (MQ=255)
                                                        cacaGCTACGATTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTAt        >  1:1983823/1‑71 (MQ=255)
                                                         acaGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTAtt       <  1:589135/71‑1 (MQ=255)
                                                           atcTACGGGTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGgaga              <  1:4087511/60‑1 (MQ=255)
                                                           aGCTACGGTTACAAAGGCCTGACCAACGAATCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTAtttt     <  1:3053831/71‑1 (MQ=255)
                                                           aGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTTAACGTAACGCTGGTAGAAGCGGGAGAACGTAtttt     <  1:2326635/71‑1 (MQ=255)
                                                             cTACGGTTACAAAGCCCTGACCAACGAAGCCCTGAACGTAACGCTGTTAGAAGCGGGAGAACGTATTTTGc   <  1:2882361/71‑1 (MQ=255)
                                                              tACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTATTTTGcc  >  1:3776959/1‑71 (MQ=255)
                                                                  |                                                                  
GCGGCGCGACGGGTGTAGAACTCTCCGCTGAATTGCACAACGCGGTCAAGCAACTGCACAGCTACGGTTACAAAGGCCTGACCAACGAAGCCCTGAACGTAACGCTGGTAGAAGCGGGAGAACGTATTTTGCC  >  W3110S.gb/1168194‑1168326

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