Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,986,216 A→G 13.3% Y81H (TAT→CAT)  slp ← outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,986,2160AG13.3% 73.3 / 5.2 30Y81H (TAT→CAT) slpouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (14/12);  new base G (2/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CAAGGAAGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAACCTTC  >  W3110S.gb/3986150‑3986286
                                                                  |                                                                      
cAAGGAAGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATGGCTa                                                                    >  1:706317/1‑71 (MQ=255)
cAAGGAAGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTa                                                                    >  1:662216/1‑71 (MQ=255)
    gAAGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATGGCt                                                                     <  1:1597512/66‑1 (MQ=255)
      aGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAAc                                                              >  1:581666/1‑71 (MQ=255)
      aGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAAc                                                              >  1:259561/1‑71 (MQ=255)
          gCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTa                                                          <  1:721170/71‑1 (MQ=255)
           cTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTaa                                                         <  1:136715/71‑1 (MQ=255)
              tGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATa                                                       >  1:1218917/1‑70 (MQ=255)
              tGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATa                                                       <  1:467909/70‑1 (MQ=255)
                    gCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGAt                                                >  1:218175/1‑71 (MQ=255)
                        gCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCt                                            <  1:277543/71‑1 (MQ=255)
                        gCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCt                                            <  1:882709/71‑1 (MQ=255)
                             ctgCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTaacaa                                       >  1:802612/3‑71 (MQ=255)
                                           ggCTTCAATATCAGGCTTCGCATAGCTATCCAACGGt                                                           <  1:1411334/37‑1 (MQ=255)
                                           ggCTTCAATATCAGGCTTCGCATAGCTATCCAACGGt                                                           <  1:5187/37‑1 (MQ=255)
                                           ggCTTCAATATCAGGCTTCGCATAGCTATACAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGc                         >  1:315425/1‑71 (MQ=255)
                                              ttCAATATCAGGCTTCGCATGGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCAt                      >  1:791195/1‑71 (MQ=255)
                                                cAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCAtt                     >  1:967748/1‑70 (MQ=255)
                                                cAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCAtt                     >  1:276708/1‑70 (MQ=255)
                                                 aaTATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTaa                                          >  1:791108/1‑48 (MQ=255)
                                                   tatCAGGCTTCGCATAGCTATCCAACGGTAATACaga                                                   <  1:1053635/37‑1 (MQ=255)
                                                   tatCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGc                         <  1:1416664/63‑1 (MQ=255)
                                                     tCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAAc               >  1:109076/1‑71 (MQ=255)
                                                             cGCATGGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGAt         <  1:959390/69‑1 (MQ=255)
                                                             cGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGAt         <  1:675795/69‑1 (MQ=255)
                                                              gCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATaa       <  1:539126/70‑1 (MQ=255)
                                                              gCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAAc      >  1:1417624/1‑71 (MQ=255)
                                                               cATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAAcc     >  1:555147/1‑71 (MQ=255)
                                                               cATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAAcc     <  1:1223464/71‑1 (MQ=255)
                                                               cATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAAcc     >  1:886109/1‑71 (MQ=255)
                                                                  aGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAACCTTc  >  1:329649/1‑71 (MQ=255)
                                                                  |                                                                      
CAAGGAAGCCGCTTTGTCTGGCGAGCAGTCGGCCCTGATAGTTGGCTTCAATATCAGGCTTCGCATAGCTATCCAACGGTAATACAGAGATTTCTAACAACGTATCCGTTTTGCCATTGATAACGTTGATAACCTTC  >  W3110S.gb/3986150‑3986286

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