Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,152,234 G→A 39.5% D6N (GAT→AAT)  yghB → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,152,2340GA39.5% 28.2 / 40.9 38D6N (GAT→AAT) yghBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (9/14);  new base A (7/8);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTT  >  W3110S.gb/3152168‑3152293
                                                                  |                                                           
aGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGatat                                                         <  1:2575030/71‑1 (MQ=255)
       cAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                                   <  1:1361824/70‑1 (MQ=255)
       cAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCt                                                  <  1:996862/71‑1 (MQ=255)
            ggCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGct                                             >  1:558894/1‑71 (MQ=255)
               aGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAatat                                                         <  1:1079378/56‑1 (MQ=255)
                   aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                      >  1:858547/1‑71 (MQ=255)
                   aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                      >  1:2109974/1‑71 (MQ=255)
                    aaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                      >  1:1302444/1‑70 (MQ=255)
                      ggCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                                   >  1:1131563/1‑71 (MQ=255)
                      ggCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAaca                                    >  1:1152228/1‑70 (MQ=255)
                       gCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                        <  1:2069144/65‑1 (MQ=255)
                        cTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                 <  1:1727199/71‑1 (MQ=255)
                        cTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGa                                 <  1:212908/71‑1 (MQ=255)
                         ttAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcatc                                                     <  1:1747502/50‑1 (MQ=255)
                           aaaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                      >  1:1934416/1‑63 (MQ=255)
                            aaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGAc                                <  1:2257818/68‑1 (MQ=255)
                               gCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGcc                          <  1:928029/71‑1 (MQ=255)
                                 tGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAac                                     <  1:1497261/58‑1 (MQ=255)
                                 tGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                        >  1:2616371/1‑71 (MQ=255)
                                     ccACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAcca                                    >  1:1401147/1‑55 (MQ=255)
                                     ccACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTg                    >  1:1417768/1‑71 (MQ=255)
                                      cacaGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                      <  1:2628794/52‑1 (MQ=255)
                                      cacaGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg                   <  1:845577/71‑1 (MQ=255)
                                        caGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCg                 >  1:2569219/1‑71 (MQ=255)
                                         aGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCACCACGACTTTGCCGCGCTGGCg                 >  1:555257/1‑70 (MQ=255)
                                          ggAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa               <  1:1072953/71‑1 (MQ=255)
                                              aGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGAt              <  1:991034/68‑1 (MQ=255)
                                                tCCATGGCTGTTATTCAAAATATCATCGCTGCGct                                             <  1:1288192/35‑1 (MQ=255)
                                                 ccATGGCTGTTATTCAAAATATCATCGCTGCGCTCTAGCAACACGACTTTGCCGCGCTGGCGGATCCTCa         <  1:1210438/70‑1 (MQ=255)
                                                  cATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCata       >  1:909996/1‑71 (MQ=255)
                                                   aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCatat      <  1:633580/71‑1 (MQ=255)
                                                   aTGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                      >  1:202191/1‑39 (MQ=255)
                                                    tGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa               >  1:1835292/1‑61 (MQ=255)
                                                      gCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGt   <  1:670715/71‑1 (MQ=255)
                                                       cTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGtt  <  1:2120559/71‑1 (MQ=255)
                                                       cTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGtt  <  1:566567/71‑1 (MQ=255)
                                                          ttATTCAAAATATCATCGCTGCGCTCTGGCAACACGAc                                <  1:420227/38‑1 (MQ=255)
                                                            aTTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                      >  1:1660812/1‑46 (MQ=255)
                                                                  |                                                           
AGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTT  >  W3110S.gb/3152168‑3152293

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