Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 772,842 G→A 13.9% G321G (GGG→GGA cydA → cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb772,8420GA13.9% 195.4 / 19.5 72G321G (GGG→GGAcydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base A (6/4);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTC  >  W3110S.gb/772775‑772911
                                                                   |                                                                     
ccGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGAt                                                                     >  1:4506432/1‑70 (MQ=255)
 cGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATg                                                                    <  1:3231393/70‑1 (MQ=255)
 cGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATg                                                                    <  1:486327/70‑1 (MQ=255)
  gTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGaa                                                                  <  1:3921716/71‑1 (MQ=255)
      aTACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCg                                                              <  1:1829540/71‑1 (MQ=255)
        aCCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTa                                                            >  1:2590190/1‑71 (MQ=255)
           ccGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActc                                                         >  1:1036348/1‑71 (MQ=255)
            cGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActct                                                        >  1:1500584/1‑71 (MQ=255)
             ggTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActctc                                                       >  1:4452400/1‑71 (MQ=255)
             ggTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActctc                                                       >  1:3844527/1‑71 (MQ=255)
              gTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActctct                                                      >  1:1731364/1‑71 (MQ=255)
              gTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActctct                                                      >  1:639878/1‑71 (MQ=255)
              gTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTActctc                                                       >  1:2837938/1‑70 (MQ=255)
                tATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAGCGCATTCGTAACGGGATGAAGGCGTACTCTCTGc                                                    >  1:1850529/1‑71 (MQ=255)
                tATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGc                                                    >  1:816677/1‑71 (MQ=255)
                 aTCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGaa                                                                  >  1:4121288/1‑56 (MQ=255)
                  tCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTa                                                            >  1:3972731/1‑61 (MQ=255)
                    ggCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCg                                                 <  1:3193464/70‑1 (MQ=255)
                        tGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTa                                                            >  1:2293831/1‑55 (MQ=255)
                          aaaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTATTCTCTGCTCGAACAACt                                          <  1:1405965/71‑1 (MQ=255)
                          aaaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGaacaa                                            <  1:3018290/69‑1 (MQ=255)
                           aaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTa                                                            <  1:2179208/52‑1 (MQ=255)
                           aaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTa                                                            <  1:1614872/52‑1 (MQ=255)
                           aaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTg                                         >  1:4375141/1‑71 (MQ=255)
                           aaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTg                                         <  1:2258906/71‑1 (MQ=255)
                           aaGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTg                                         <  1:2805585/71‑1 (MQ=255)
                             gagCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCg                                       >  1:479594/1‑71 (MQ=255)
                              agCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGt                                      >  1:2127138/1‑71 (MQ=255)
                                cTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATg                                                                    <  1:2608098/39‑1 (MQ=255)
                                cTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTc                                    <  1:3494819/71‑1 (MQ=255)
                                 tGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTc                                    <  1:1000521/70‑1 (MQ=255)
                                  gATGGTGCAGCATGAAGAACGCATTCGTAACGGGATg                                                                    <  1:3184336/37‑1 (MQ=255)
                                  gATGGTGCAGCATGAAGAACGCATTCGTAACGGGATg                                                                    <  1:2746842/37‑1 (MQ=255)
                                    tGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGt                                >  1:40583/1‑71 (MQ=255)
                                          agcaTGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAAGAACTGCGTTCTGGTTCTAcc                          <  1:1786297/71‑1 (MQ=255)
                                          agcaTGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTAcc                          <  1:2152436/71‑1 (MQ=255)
                                           gcaTGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCg                         <  1:2299476/71‑1 (MQ=255)
                                            caTGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTc                                    >  1:3304847/1‑59 (MQ=255)
                                            caTGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGa                        >  1:1606913/1‑71 (MQ=255)
                                             aTGAAGAACGCATTCGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGAc                       >  1:1725201/1‑71 (MQ=255)
                                             aTGAAGAACGCATTCGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGAc                       >  1:3659125/1‑71 (MQ=255)
                                              tGAAGAACGCATTCGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGAc                       >  1:3638903/1‑70 (MQ=255)
                                               gaagaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGa                        >  1:4130551/1‑68 (MQ=255)
                                                aagaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAg                    >  1:4461582/1‑71 (MQ=255)
                                                aagaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAg                    >  1:3529272/1‑71 (MQ=255)
                                                  gaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTAcc                          <  1:3835193/63‑1 (MQ=255)
                                                  gaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGAc                       >  1:2342336/1‑66 (MQ=255)
                                                  gaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAga                   >  1:2054973/1‑69 (MQ=255)
                                                  gaaCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGc                  <  1:2649124/71‑1 (MQ=255)
                                                   aaCGCATTCGTAACTGGATGAAGGCGTACTCTCTGCTCGaaca                                             <  1:789045/43‑1 (MQ=255)
                                                    aCGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGa                        >  1:4401553/1‑63 (MQ=255)
                                                       cATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCg                         <  1:3177711/59‑1 (MQ=255)
                                                        aTTCGTAACGGAATGAAGGCGTACTCTCTGCTCGa                                                <  1:553619/35‑1 (MQ=255)
                                                         ttCGTAACGGGATGAAGGCGTACTCTCTGCTCGaac                                              <  1:4515159/36‑1 (MQ=255)
                                                         ttCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCa                     <  1:981492/61‑1 (MQ=255)
                                                         ttCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCa                     <  1:2409122/61‑1 (MQ=255)
                                                         ttCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTc             <  1:2717951/69‑1 (MQ=255)
                                                         ttCGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGt           <  1:2550478/71‑1 (MQ=255)
                                                          tCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGt                                >  1:626760/1‑49 (MQ=255)
                                                           cGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCa                     >  1:335174/1‑59 (MQ=255)
                                                           cGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGa         >  1:4392651/1‑71 (MQ=255)
                                                           cGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACAAGGCGGTTCGTGa         >  1:1630366/1‑71 (MQ=255)
                                                           cGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGa         >  1:966120/1‑71 (MQ=255)
                                                           cGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGa         >  1:3065650/1‑71 (MQ=255)
                                                           cGTAACGGAATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGa         >  1:659234/1‑71 (MQ=255)
                                                             tAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGAcc       <  1:3251044/71‑1 (MQ=255)
                                                               aCGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTa                            >  1:1721756/1‑48 (MQ=255)
                                                               aCGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTAc                           >  1:742801/1‑49 (MQ=255)
                                                               aCGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAg     <  1:2111975/71‑1 (MQ=255)
                                                               aCGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAg     >  1:1858144/1‑71 (MQ=255)
                                                               aCGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAg     <  1:2925619/71‑1 (MQ=255)
                                                                cGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGAc                       <  1:42494/52‑1 (MQ=255)
                                                                  ggATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTc  <  1:474843/71‑1 (MQ=255)
                                                                  ggATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTc  <  1:3881896/71‑1 (MQ=255)
                                                                   |                                                                     
CCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTC  >  W3110S.gb/772775‑772911

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: