Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,567,579 T→C 15.8% N313D (AAT→GAT)  yddV ← predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,567,5790TC15.8% 89.4 / 10.7 38N313D (AAT→GAT) yddVpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (14/18);  new base C (4/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGTTTT  >  W3110S.gb/1567510‑1567648
                                                                     |                                                                     
cTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTa                                                                      <  1:1377233/71‑1 (MQ=255)
                   cAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACa                                                                  >  1:3017287/1‑56 (MQ=255)
                    agagCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATgg                                                   >  1:242950/1‑70 (MQ=255)
                     gagCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCATGATGACCATTAAACATGGCGGCACCGATGGaa                                                 <  1:1043281/71‑1 (MQ=255)
                      agCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGaaa                                                <  1:3948016/71‑1 (MQ=255)
                        cTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATCAAACa                                                                  >  1:512055/1‑51 (MQ=255)
                          tCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGt                                              >  1:2223850/1‑69 (MQ=255)
                           cATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTga                                            <  1:543629/70‑1 (MQ=255)
                           cATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTga                                            <  1:4002596/70‑1 (MQ=255)
                            aTCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTgag                                           <  1:1446262/70‑1 (MQ=255)
                             tCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCg                                                       <  1:138912/57‑1 (MQ=255)
                             tCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCg                                                       <  1:313096/57‑1 (MQ=255)
                             tCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCg                                                       <  1:3834313/57‑1 (MQ=255)
                                 cTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATggcgg                                                            <  1:3993317/48‑1 (MQ=255)
                                    tttGAATGAGGCGCTCATAGTCAGGATGACCATTaaa                                                                    <  1:1176621/37‑1 (MQ=255)
                                      tGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGaaa                                                >  1:2929403/1‑55 (MQ=255)
                                         aTGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTc                             >  1:1906908/1‑71 (MQ=255)
                                              gcgcTCATAGTCAGGATGACCATCAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGt                        <  1:1366001/71‑1 (MQ=255)
                                                 cTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTg                     >  1:1065866/1‑71 (MQ=255)
                                                 cTCATAGTCAGGATGACCATCAAACATGGCGGCACCGATGGAAAGTg                                             <  1:1175870/47‑1 (MQ=255)
                                                   cATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGaaa                                                >  1:4512264/1‑42 (MQ=255)
                                                   cATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGaaa                                                >  1:4637944/1‑42 (MQ=255)
                                                   cATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCg                   >  1:2636707/1‑71 (MQ=255)
                                                   cATAGTCAGGATGACCATCAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCg                   >  1:906014/1‑71 (MQ=255)
                                                   cATAGTCAGGATGACCATCAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCg                   >  1:3722107/1‑71 (MQ=255)
                                                   cATAGTCAGGATGACCATCAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCg                   >  1:3636421/1‑71 (MQ=255)
                                                    aTAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCg                   >  1:3023484/1‑70 (MQ=255)
                                                      aGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGagag                                         >  1:3896053/1‑46 (MQ=255)
                                                       gTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTgag                                           <  1:2151116/43‑1 (MQ=255)
                                                       gTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCtt               >  1:2355273/1‑71 (MQ=255)
                                                       gTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCtt               >  1:4580841/1‑71 (MQ=255)
                                                       gTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCtt               <  1:2825600/71‑1 (MQ=255)
                                                          aGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCTCCGTTTGCGGCttt              >  1:2577430/1‑69 (MQ=255)
                                                           ggATGACCATTAAACATGGCGGCACCGATGGAAAGt                                              <  1:2604553/36‑1 (MQ=255)
                                                             aTGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGtt         <  1:1439636/71‑1 (MQ=255)
                                                              tGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGttt        <  1:4457280/71‑1 (MQ=255)
                                                                 ccATTAAACATGGCGGCACCGATGGATAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAg           <  1:2149046/65‑1 (MQ=255)
                                                                 ccATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGt     <  1:4313249/71‑1 (MQ=255)
                                                                    ttAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGt          <  1:1410423/63‑1 (MQ=255)
                                                                    ttAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGtttt  <  1:2474604/71‑1 (MQ=255)
                                                                    ttAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGtttt  <  1:1564633/71‑1 (MQ=255)
                                                                     |                                                                     
CTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGTTTT  >  W3110S.gb/1567510‑1567648

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: