Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,267,116 G→C 7.0% V236V (GTC→GTG qor ← quinone oxidoreductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,267,1160GC 7.0% 159.9 / 2.0 57V236V (GTC→GTGqorquinone oxidoreductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (20/32);  new base C (2/2);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AACGAGCCTTTTTGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACA  >  W3110S.gb/4267048‑4267186
                                                                    |                                                                      
aaCGAGCCTTTTTGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGAc                                                                      <  1:3900355/71‑1 (MQ=255)
 aCGAGCCTTTTTGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGAc                                                                      >  1:318935/1‑70 (MQ=255)
 aCGAGCCTTTTTGATTGAGAATGCCTAAGTTCACACCAGTAACCGCACCTGATGAGTTGCCAAAACTGAc                                                                      >  1:120167/1‑70 (MQ=255)
    aGCCTTTTTGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCAtt                                                                  <  1:3809543/71‑1 (MQ=255)
      ccTTTTTGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGAcc                                                                     <  1:3164/66‑1 (MQ=255)
          tttGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGc                                                              >  1:856152/1‑69 (MQ=255)
           ttgattgaGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGcc                                                             >  1:2645060/1‑69 (MQ=255)
           ttgattgaGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCcgc                                                           >  1:2988674/1‑71 (MQ=255)
           ttgattgaGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCcgc                                                           >  1:3738911/1‑71 (MQ=255)
           ttgattgaGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCcgc                                                           >  1:1698503/1‑71 (MQ=255)
              attgaGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGAc                                                                      <  1:749342/57‑1 (MQ=255)
                  agaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTg                                                    <  1:497100/71‑1 (MQ=255)
                   gaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTg                                                    <  1:1613341/70‑1 (MQ=255)
                      tGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCcgg                                                           >  1:2222515/1‑59 (MQ=255)
                      tGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGtt                                                     <  1:1661813/66‑1 (MQ=255)
                          tAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGc                                                   <  1:1575356/64‑1 (MQ=255)
                          tAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGc                                                   <  1:1082203/64‑1 (MQ=255)
                           aaGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTTACCATTAAGCCGCGGCGTTGCAGGCAATc                                           >  1:3516420/1‑71 (MQ=255)
                            aGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATc                                           <  1:2977048/70‑1 (MQ=255)
                              ttCACACCGGTAACCGCACCTGATGAGTTCCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAg                                        <  1:1067695/71‑1 (MQ=255)
                                cacaCCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAg                                                 >  1:1753439/1‑60 (MQ=255)
                                cacaCCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGc                                       <  1:2086871/70‑1 (MQ=255)
                                      ggTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAg                                                               >  1:3474551/1‑40 (MQ=255)
                                       gTAACCGCACCTGATGAGTTGCCAAAACTCACCATTAAGCCGCGGCGTTGCAGGCAATCCAGc                                       >  1:1003902/1‑63 (MQ=255)
                                       gTAACCGCACCTGATGAGTTGCCAAAACTCACCATTAAGCCGCGGCGTTGCAGGCAATCCAGc                                       >  1:2634261/1‑63 (MQ=255)
                                        tAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTcc                              <  1:2466235/71‑1 (MQ=255)
                                        tAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTc                               <  1:3351638/70‑1 (MQ=255)
                                         aaCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTcc                              <  1:2471141/70‑1 (MQ=255)
                                         aaCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTcc                              >  1:2555851/1‑70 (MQ=255)
                                          aCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTccc                             >  1:3037898/1‑70 (MQ=255)
                                               aCCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGgtgt                       <  1:1668353/71‑1 (MQ=255)
                                                 cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTg                                                    <  1:1221528/40‑1 (MQ=255)
                                                 cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTg                                                    <  1:2976042/40‑1 (MQ=255)
                                                 cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGt                                 <  1:1169531/59‑1 (MQ=255)
                                                 cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGtct                     <  1:2158857/71‑1 (MQ=255)
                                                 cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGtct                     >  1:1557471/1‑71 (MQ=255)
                                                  tgatgaGTTGCCAAAACTGACCATTAAGCCgcggc                                                        <  1:2641406/35‑1 (MQ=255)
                                                    atgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTg                  >  1:3297683/1‑71 (MQ=255)
                                                     tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGt                                 >  1:1051862/1‑55 (MQ=255)
                                                     tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGt                                 >  1:4347498/1‑55 (MQ=255)
                                                       aGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGcc                <  1:1608131/70‑1 (MQ=255)
                                                        gTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTcc                              <  1:2745450/55‑1 (MQ=255)
                                                        gTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTcc                              <  1:1002847/55‑1 (MQ=255)
                                                            ccAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACgg           <  1:656710/70‑1 (MQ=255)
                                                            ccAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGa          <  1:1688595/71‑1 (MQ=255)
                                                             cAAAACTCACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTGCCAg                           <  1:3988228/53‑1 (MQ=39)
                                                             cAAAACTCACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAg                           <  1:2147906/53‑1 (MQ=255)
                                                              aaaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGgt                         <  1:2780400/54‑1 (MQ=255)
                                                              aaaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGgt                         <  1:3491752/54‑1 (MQ=255)
                                                               aaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGgtgt                       >  1:2000377/1‑55 (MQ=255)
                                                               aaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATc       >  1:3502498/1‑71 (MQ=255)
                                                                aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCa                                         >  1:2538134/1‑36 (MQ=255)
                                                                aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCg      <  1:3700921/71‑1 (MQ=255)
                                                                aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCg      <  1:2368341/71‑1 (MQ=255)
                                                                aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCg      <  1:1505449/71‑1 (MQ=255)
                                                                aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCg      <  1:988593/71‑1 (MQ=255)
                                                                 aCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCg      >  1:733928/1‑70 (MQ=255)
                                                                    gACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTaca  >  1:548198/1‑71 (MQ=255)
                                                                    gACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTaca  >  1:56200/1‑71 (MQ=255)
                                                                    |                                                                      
AACGAGCCTTTTTGATTGAGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACA  >  W3110S.gb/4267048‑4267186

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: