Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,180,003 T→C 12.7% D369D (GAT→GAC lolE → outer membrane‑specific lipoprotein transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,180,0030TC12.7% 136.7 / 10.4 55D369D (GAT→GAClolEouter membrane‑specific lipoprotein transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (32/16);  new base C (2/5);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATTATTGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCATTGTTGCTGAGTCTTT  >  W3110S.gb/1179935‑1180073
                                                                    |                                                                      
ccGATTATTGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGAttt                                                                      <  1:1735908/71‑1 (MQ=255)
ccGATTATTGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGAtt                                                                       >  1:2856908/1‑70 (MQ=255)
     tattGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGc                                                                 >  1:1399970/1‑71 (MQ=255)
     tattGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGc                                                                 >  1:2614801/1‑71 (MQ=255)
      attGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGc                                                                 >  1:2120215/1‑70 (MQ=255)
      attGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGc                                                                 <  1:2549534/70‑1 (MQ=255)
             ggATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:1765300/1‑71 (MQ=255)
                ttGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGGGATATCTATTTTATTGATTTCCTGCCATCGGAATTg                                                      <  1:579339/71‑1 (MQ=255)
                    aaaGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCAct                                                  >  1:2366105/1‑71 (MQ=255)
                        ttGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGAttt                                                                      <  1:404870/47‑1 (MQ=255)
                        ttGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGAttt                                                                      <  1:2828008/47‑1 (MQ=255)
                        ttGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGAttt                                                                      <  1:1919619/47‑1 (MQ=255)
                        ttGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCgg                                                           <  1:1368625/58‑1 (MQ=255)
                          gATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTg                                                                  >  1:3315843/1‑49 (MQ=255)
                          gATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCActggctg                                             >  1:2112716/1‑70 (MQ=255)
                                tcatcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCAatggctgg                                            <  1:3766986/65‑1 (MQ=255)
                                tcatcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtct                                      >  1:1770283/1‑71 (MQ=255)
                                tcatcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtct                                      >  1:1685502/1‑71 (MQ=255)
                                  atgaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    >  1:740760/4‑71 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCActggc                                               >  1:2191472/1‑60 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    >  1:3745098/1‑71 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    >  1:2811975/1‑71 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    >  1:2823371/1‑71 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    <  1:2610725/71‑1 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    <  1:3315942/71‑1 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    <  1:2397375/71‑1 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    <  1:1983618/71‑1 (MQ=255)
                                  atcaGTTCCTCTCCAGCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    <  1:1974488/71‑1 (MQ=255)
                                    caGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:3830404/1‑48 (MQ=255)
                                    caGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:554669/1‑48 (MQ=255)
                                     aGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:3025475/1‑47 (MQ=255)
                                     aGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:319719/1‑47 (MQ=255)
                                     aGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:601981/1‑47 (MQ=255)
                                     aGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCCGCCATCGGaa                                                         >  1:3746125/1‑47 (MQ=255)
                                      gTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACg                                <  1:2515343/71‑1 (MQ=255)
                                      gTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACg                                <  1:695462/71‑1 (MQ=255)
                                      gTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACg                                <  1:3092063/71‑1 (MQ=255)
                                       ttCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGc                                                     <  1:829231/49‑1 (MQ=255)
                                         ccTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTAc                             >  1:605346/1‑71 (MQ=255)
                                         ccTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTAc                             >  1:1285839/1‑71 (MQ=255)
                                         ccTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTAc                             >  1:2769624/1‑71 (MQ=255)
                                             tcCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGt                         >  1:2536646/1‑71 (MQ=255)
                                             tcCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGt                         >  1:1508160/1‑71 (MQ=255)
                                             tcCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGt                         >  1:1797312/1‑71 (MQ=255)
                                             tcCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCATCTACGTACTGGt                         >  1:3757653/1‑71 (MQ=255)
                                               cAGCGATATCTATTTTATTGATTTCCTGCCATCGGaa                                                         >  1:2106581/1‑37 (MQ=255)
                                               cAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTa                              >  1:1722639/1‑64 (MQ=255)
                                                aGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTcac                      <  1:936204/71‑1 (MQ=255)
                                                aGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTcac                      <  1:345063/71‑1 (MQ=255)
                                                 gCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCa                                                    >  1:2589978/1‑40 (MQ=37)
                                                 gCGATATCTATTTTATTGACTTCTTGCCATCGGAATTGCa                                                    >  1:608689/1‑40 (MQ=37)
                                                  cGATATCTATTTTATTGATTTCCTGCCATAGGAATTGCACTGGCTGGACGTCTTCTaag                                <  1:186727/59‑3 (MQ=255)
                                                     tatCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCAt                 <  1:1258762/71‑1 (MQ=255)
                                                       tCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTAc                                 >  1:481371/1‑53 (MQ=255)
                                                       tCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACg                                <  1:3061476/54‑1 (MQ=255)
                                                       tCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACg                                <  1:2205723/54‑1 (MQ=255)
                                                                ttGATTTCCTGCCATCGGAATTGCACTGGCTGGACGtcttc                                    <  1:1295989/41‑1 (MQ=255)
                                                                 tGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCa                  >  1:2373375/1‑58 (MQ=255)
                                                                 tGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCa                  >  1:3122599/1‑58 (MQ=255)
                                                                   aTTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCATTGTTGCTGAGTCt    <  1:940211/70‑1 (MQ=255)
                                                                    tttCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCATTGTTGCTGAGTCttt  >  1:2139018/1‑71 (MQ=255)
                                                                    |                                                                      
CCGATTATTGAGTGGATTGAAAAGTTGATCGGTCATCAGTTCCTCTCCAGCGATATCTATTTTATTGATTTCCTGCCATCGGAATTGCACTGGCTGGACGTCTTCTACGTACTGGTCACAGCATTGTTGCTGAGTCTTT  >  W3110S.gb/1179935‑1180073

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: