Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,773,435 T→C 22.5% L488L (TTG→CTG)  yfjW → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,773,4350TC22.5% 77.0 / 20.4 40L488L (TTG→CTG) yfjWpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (7/24);  new base C (3/6);  total (10/30)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGT  >  W3110S.gb/2773374‑2773499
                                                             |                                                                
tATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTCTg                                                                <  1:649404/64‑1 (MQ=255)
 aTCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTCTg                                                                <  1:1836685/63‑1 (MQ=255)
  tCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:802109/62‑1 (MQ=255)
  tCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:218262/62‑1 (MQ=255)
  tCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGAc                                                         >  1:691845/1‑69 (MQ=255)
  tCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGAc                                                         >  1:1466892/1‑69 (MQ=255)
  tCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGAc                                                         >  1:1629733/1‑69 (MQ=255)
    aaCGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTa                                                              >  1:1325782/1‑62 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:433601/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:295766/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:237974/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:316529/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:189871/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:1892757/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:1690732/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:1674581/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:394894/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:146792/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTg                                                                <  1:1037002/54‑1 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTCTGTATGGACTAGTtaata                                                >  1:1414169/1‑70 (MQ=255)
          aaaCATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTCTGTATGGACTAGTtaa                                                  >  1:1218456/1‑68 (MQ=255)
                        aaGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTa                                                              >  1:1644610/1‑42 (MQ=255)
                        aaGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTa                                                              >  1:537547/1‑42 (MQ=255)
                        aaGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTa                                                              >  1:311855/1‑42 (MQ=255)
                        aaGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTCTGTATGGACTAGTTAATAAAAGCAAAAAAATc                                  >  1:51360/1‑70 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:1092121/70‑1 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:735485/70‑1 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:1924957/70‑1 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:140075/70‑1 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:1385612/70‑1 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:1357010/70‑1 (MQ=255)
                                    ggAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGc                      <  1:135245/70‑1 (MQ=255)
                                           atttttTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCAcagcag               <  1:1625826/70‑1 (MQ=255)
                                           attattTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCAcagcag               <  1:1951277/70‑1 (MQ=255)
                                           attattTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCAcagcag               <  1:472720/70‑1 (MQ=255)
                                           attattTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCAcagcag               <  1:116473/70‑1 (MQ=255)
                                                       gggAGTCTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTttgt  <  1:1493858/71‑1 (MQ=255)
                                                       gggAGTCTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTttgt  <  1:448305/71‑1 (MQ=255)
                                                       gggAGTCTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTttgt  <  1:1132602/71‑1 (MQ=255)
                                                       gggAGTCTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTttgt  <  1:833515/71‑1 (MQ=255)
                                                             |                                                                
TATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGT  >  W3110S.gb/2773374‑2773499

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: