Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,152,234 G→A 26.9% 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,2340GA26.9% 139.1 / 55.3 78D6N (GAT→AAT) yghBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (9/12);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

TTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCGTTGTTT  >  W3110S.gb/3152165‑3152303
                                                                     |                                                                     
ttCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                                      <  1:2864333/71‑1 (MQ=255)
  cAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                                     <  1:3790381/70‑1 (MQ=255)
  cAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                                     <  1:2596218/70‑1 (MQ=255)
  cAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                                     <  1:1656952/70‑1 (MQ=255)
      aGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcat                                                                >  1:2382956/1‑71 (MQ=255)
      aGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcat                                                                >  1:1391355/1‑71 (MQ=255)
          cAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGc                                                             >  1:3199593/1‑70 (MQ=255)
           aaTAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                                      <  1:3163703/60‑1 (MQ=255)
            aTAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgc                                                          >  1:1491770/1‑71 (MQ=255)
             tAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatca                                                                 >  1:3123607/1‑63 (MQ=255)
             tAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcg                                                         >  1:1035728/1‑71 (MQ=255)
               ggCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGct                                                       >  1:2796520/1‑71 (MQ=255)
                    aTAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcat                                                                <  1:3025371/57‑1 (MQ=255)
                    aTAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                                  <  1:3847470/55‑1 (MQ=255)
                    aTAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGc                                                  <  1:181614/71‑1 (MQ=255)
                    aTAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGc                                                  <  1:2143954/71‑1 (MQ=255)
                     tAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGc                                                  <  1:4769018/70‑1 (MQ=255)
                     tAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGc                                                  >  1:2162360/1‑70 (MQ=255)
                      aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                                             >  1:2519501/1‑58 (MQ=255)
                      aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                                             >  1:1291089/1‑58 (MQ=255)
                      aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:1125350/1‑71 (MQ=255)
                      aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:938842/1‑71 (MQ=255)
                      aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:1140102/1‑71 (MQ=255)
                         ggCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                                             <  1:1129887/71‑1 (MQ=255)
                          gCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                                >  1:4005755/1‑67 (MQ=255)
                            ttAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                                      >  1:249594/1‑43 (MQ=255)
                            ttAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                                      >  1:665353/1‑43 (MQ=255)
                            ttAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcatc                                                               >  1:2348490/1‑50 (MQ=255)
                            ttAAACGCTGTTACACAGGAAAGTCCATGGATGTTATTCAAGATatcatc                                                               >  1:2689817/1‑50 (MQ=255)
                             tAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                                  <  1:3361487/46‑1 (MQ=255)
                               aaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                                             <  1:4499372/65‑1 (MQ=255)
                               aaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACttt                                       <  1:3565152/71‑1 (MQ=255)
                               aaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACttt                                       <  1:3120010/71‑1 (MQ=255)
                               aaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACttt                                       <  1:898450/71‑1 (MQ=255)
                                aCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACttt                                       >  1:4426438/1‑70 (MQ=255)
                                    tGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCg                                                              <  1:3146970/43‑1 (MQ=255)
                                       tCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               >  1:174312/1‑71 (MQ=255)
                                        ccACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               <  1:3568291/70‑1 (MQ=255)
                                        ccACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               <  1:2569805/70‑1 (MQ=255)
                                            aGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCgg                          <  1:2320292/71‑1 (MQ=255)
                                             ggAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa                         <  1:808884/71‑1 (MQ=255)
                                             ggAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa                         <  1:4491625/71‑1 (MQ=255)
                                               aaaGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATc                       >  1:2705516/1‑71 (MQ=255)
                                               aaaGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATc                       >  1:2195572/1‑71 (MQ=255)
                                               aaaGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATc                       >  1:4164781/1‑71 (MQ=255)
                                               aaaGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATc                       >  1:244733/1‑71 (MQ=255)
                                                 aGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCt                     <  1:944159/71‑1 (MQ=255)
                                                  gTCCATGGCTGTTATTCAAAATATCATCGCTGCGctct                                                     <  1:3798289/38‑1 (MQ=255)
                                                   tCCATGGCTGTTATTCAAGATATCATCGCTGCGctc                                                      <  1:3578029/36‑1 (MQ=255)
                                                   tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                                >  1:241989/1‑58 (MQ=255)
                                                   tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTc                    >  1:3564742/1‑70 (MQ=255)
                                                   tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa                   <  1:1134164/71‑1 (MQ=255)
                                                   tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa                   >  1:4068058/1‑71 (MQ=255)
                                                   tCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACttt                                       <  1:3370988/51‑1 (MQ=255)
                                                   tCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACttt                                       <  1:3885475/51‑1 (MQ=255)
                                                     cATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                           <  1:3873819/45‑1 (MQ=255)
                                                     cATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCata                 <  1:1077682/71‑1 (MQ=255)
                                                     cATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACttt                                       >  1:284845/1‑49 (MQ=255)
                                                     cATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACttt                                       >  1:3715270/1‑49 (MQ=255)
                                                      aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                           <  1:2500380/44‑1 (MQ=255)
                                                      aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                                  <  1:4486063/53‑1 (MQ=255)
                                                      aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                                  <  1:2290036/53‑1 (MQ=255)
                                                        ggCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTg              >  1:753445/1‑71 (MQ=255)
                                                         gCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGt             <  1:85237/71‑1 (MQ=255)
                                                          cTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGtt            <  1:4398843/71‑1 (MQ=255)
                                                            gTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAg          >  1:4074318/1‑71 (MQ=255)
                                                             ttATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGc         <  1:1022955/71‑1 (MQ=255)
                                                              tATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCg        >  1:4651141/1‑71 (MQ=255)
                                                              tATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCg        >  1:2481341/1‑71 (MQ=255)
                                                               aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                                >  1:861173/1‑46 (MQ=255)
                                                               aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgt       >  1:2641138/1‑71 (MQ=255)
                                                               aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgt       >  1:4786224/1‑71 (MQ=255)
                                                                ttCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGtt            <  1:1586317/65‑1 (MQ=255)
                                                                ttCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgtt      <  1:4072629/71‑1 (MQ=255)
                                                                ttCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgtt      <  1:1308122/71‑1 (MQ=255)
                                                                  cAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAg          >  1:1631181/1‑65 (MQ=255)
                                                                  cAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgttg     <  1:4807942/70‑1 (MQ=255)
                                                                    aGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCGTTGttt  <  1:807490/71‑1 (MQ=255)
                                                                     |                                                                     
TTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCGTTGTTT  >  W3110S.gb/3152165‑3152303

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: