Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,375,433 A→G 30.2% D237G (GAT→GGT)  dcuD → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,375,4330AG30.2% 84.6 / 53.0 63D237G (GAT→GGT) dcuDpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (19/25);  new base G (11/8);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTC  >  W3110S.gb/3375366‑3375503
                                                                   |                                                                      
gTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATa                                                                      <  1:4374672/70‑1 (MQ=255)
gTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATa                                                                      <  1:2964646/70‑1 (MQ=255)
  gCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATg                                                                   >  1:3888497/1‑71 (MQ=255)
   cAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGt                                                                  >  1:122635/1‑71 (MQ=255)
       gcgcTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCccg                                                              >  1:2763013/1‑71 (MQ=255)
            tttGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGct                                                         >  1:2318998/1‑71 (MQ=255)
            tttGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATaa                                                                     <  1:1806704/59‑1 (MQ=255)
            tttGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATaa                                                                     <  1:2411552/59‑1 (MQ=255)
              tGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGctc                                                        <  1:2662332/70‑1 (MQ=255)
              tGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGctc                                                        <  1:3510279/70‑1 (MQ=255)
               gACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGctct                                                       <  1:1373639/70‑1 (MQ=255)
               gACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTa                                                      >  1:4667394/1‑71 (MQ=255)
                aCAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTAt                                                     <  1:2931568/71‑1 (MQ=255)
                  aaaaaaGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAAt                                                                    >  1:1684109/1‑54 (MQ=255)
                     aaaGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGCTCTATTACGc                                                >  1:654485/1‑71 (MQ=255)
                     aaaGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGCTCTAATACGc                                                >  1:417564/1‑71 (MQ=255)
                      aaGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGcc                                               <  1:1399226/71‑1 (MQ=255)
                       aGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCa                                              <  1:1297317/71‑1 (MQ=255)
                       aGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCa                                              <  1:5121554/71‑1 (MQ=255)
                       aGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCa                                              >  1:1863888/1‑71 (MQ=255)
                       aGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCa                                              >  1:712757/1‑71 (MQ=255)
                        gATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCAt                                             <  1:4908609/71‑1 (MQ=255)
                        gATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCccgc                                                             >  1:1845516/1‑55 (MQ=255)
                          tAAAAATATCAATCACGAACAGGCAGAGCAAAACGCTCTCGATAATGTCCCGCCGCTCTATTACGCCAttt                                           >  1:3283646/1‑71 (MQ=255)
                            aaaaTATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTa                                                                      <  1:4175590/42‑1 (MQ=255)
                            aaaaTATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTa                                                                      <  1:5129943/42‑1 (MQ=255)
                                tatCAATCACGAACAGGCAGAGCAAAAAGCTCTCGGTAATGTCCCGCCGctc                                                        <  1:5038946/52‑1 (MQ=255)
                                   caatcaCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCACGCCGCTCTATTACGCCATTTTACCTGTGa                                  >  1:4962070/1‑71 (MQ=255)
                                     atcaCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATg                                <  1:2687800/71‑1 (MQ=255)
                                      tcaCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGctct                                                       <  1:642192/47‑1 (MQ=255)
                                        aCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCAtt                                            >  1:1657923/1‑56 (MQ=255)
                                        aCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCg                             >  1:641502/1‑71 (MQ=255)
                                        aCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCg                             >  1:2137050/1‑71 (MQ=255)
                                              aGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATc                       >  1:1464616/1‑71 (MQ=255)
                                                  agagCAAAAAGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGa                   <  1:1933898/71‑1 (MQ=255)
                                                      cAAAAAGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgct               >  1:3896359/1‑71 (MQ=255)
                                                      cAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgct               <  1:3042368/71‑1 (MQ=255)
                                                       aaaaaGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCATTTTAc                                        >  1:1305759/1‑45 (MQ=255)
                                                        aaaaGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcg             >  1:1792984/1‑71 (MQ=255)
                                                        aaaaGCTCTCGATAATGTCCCGCCGCTCTATTACGCCAtttt                                          <  1:3774727/42‑1 (MQ=255)
                                                        aaaaGCTCTCGATAATGTCCCGCCGCTCTATTACGCCAtttt                                          <  1:806232/42‑1 (MQ=255)
                                                        aaaaGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcg             <  1:1157571/71‑1 (MQ=255)
                                                        aaaaGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcg             >  1:3326921/1‑71 (MQ=255)
                                                         aaaGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCATTTTACCtg                                     >  1:3440307/1‑46 (MQ=255)
                                                         aaaGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATcc                      <  1:3133302/61‑1 (MQ=255)
                                                         aaaGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgc                <  1:2738429/67‑1 (MQ=255)
                                                          aaGCTCTCGGTAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGa                   >  1:5142578/1‑63 (MQ=255)
                                                           aGCTCTCGATAATGTCCCGCCGCTCTATTACTCCATTTTAc                                        >  1:1436819/1‑41 (MQ=255)
                                                           aGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTAc                                        >  1:408999/1‑41 (MQ=255)
                                                           aGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTAc                                        >  1:2371317/1‑41 (MQ=255)
                                                           aGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTAc                                        >  1:2507334/1‑41 (MQ=255)
                                                           aGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctc              >  1:4837081/1‑67 (MQ=255)
                                                             ctctCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           <  1:1009535/52‑1 (MQ=255)
                                                             ctctCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           <  1:3366369/52‑1 (MQ=255)
                                                             ctctCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           <  1:4359920/52‑1 (MQ=255)
                                                                tcGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcgcc           <  1:3241699/65‑3 (MQ=255)
                                                                 cGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGt    >  1:3150674/1‑71 (MQ=255)
                                                                  gATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           <  1:4413723/47‑1 (MQ=255)
                                                                  gATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           <  1:1455873/47‑1 (MQ=255)
                                                                  gATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           >  1:4131943/1‑47 (MQ=255)
                                                                  gATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGtt                           <  1:4059435/47‑1 (MQ=255)
                                                                  gATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGt    <  1:5086156/70‑1 (MQ=255)
                                                                  gATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGt    <  1:3607723/70‑1 (MQ=255)
                                                                   aTAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTc  <  1:2014035/71‑1 (MQ=255)
                                                                   |                                                                      
GTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTC  >  W3110S.gb/3375366‑3375503

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: