Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 542,820 C→T 37.8% A151T (GCC→ACC)  ylbA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb542,8200CT37.8% 28.0 / 35.8 37A151T (GCC→ACC) ylbAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (10/13);  new base T (6/8);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

TTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTG  >  W3110S.gb/542752‑542889
                                                                    |                                                                     
tttGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGc                                                                       <  1:1435223/69‑1 (MQ=255)
          aaaTCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAAc                                                           >  1:903125/1‑71 (MQ=255)
          aaaTCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAAc                                                           <  1:1843399/71‑1 (MQ=255)
             tCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAg                                                        <  1:2660452/71‑1 (MQ=255)
                agcagAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCa                                                     <  1:2439829/71‑1 (MQ=255)
                 gcagAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCAc                                                    <  1:3649338/71‑1 (MQ=255)
                   agAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACgg                                                  >  1:1203271/1‑71 (MQ=255)
                    gAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGc                                                                       <  1:3663510/49‑1 (MQ=255)
                    gAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACgg                                                  >  1:270445/1‑70 (MQ=255)
                       tAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCa                                                                      >  1:644559/1‑47 (MQ=255)
                       tAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAAcc                                                          <  1:302299/59‑1 (MQ=255)
                            tCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCAc                                                    <  1:1394654/60‑1 (MQ=255)
                            tCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACg                                                   <  1:4616714/61‑1 (MQ=255)
                             cGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGa                                                              <  1:553855/49‑1 (MQ=255)
                              gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCtt                                       >  1:4127261/1‑71 (MQ=255)
                                 cATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCt                                     <  1:1405264/70‑1 (MQ=255)
                                 cATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCa                                                     >  1:4229663/1‑54 (MQ=255)
                                    gCCTTCATAATGAATGCGTTCCAGTTCGCTGGTAt                                                                     >  1:2057090/1‑35 (MQ=255)
                                      cTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGa                               <  1:4272108/71‑1 (MQ=255)
                                       ttCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGa                                                              <  1:1639118/39‑1 (MQ=255)
                                       ttCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGt                                                 <  1:3022160/52‑1 (MQ=255)
                                        tCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAAc                                                           <  1:1714735/41‑1 (MQ=255)
                                        tCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAAc                                                           <  1:4206777/41‑1 (MQ=255)
                                        tCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCa                                                     <  1:4364957/47‑1 (MQ=255)
                                        tCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa                                    >  1:547830/1‑64 (MQ=255)
                                         cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGaaa                                                            >  1:762415/1‑39 (MQ=255)
                                         cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGaaa                                                            >  1:1879565/1‑39 (MQ=255)
                                         cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa                                    >  1:627842/1‑63 (MQ=255)
                                          aTAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACAt                           <  1:781388/71‑1 (MQ=255)
                                          aTAATGAATGCGTTCCAGTTCGCTGGCATTGCCATAAACCAGCCACGGTGCATAGCCTTCTACCGGAACa                            <  1:558517/70‑1 (MQ=255)
                                          aTAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACa                            >  1:4339035/1‑70 (MQ=255)
                                             atgaatgCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa                                    >  1:2517252/1‑59 (MQ=255)
                                              tgaatgCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATa                                            >  1:937225/1‑50 (MQ=255)
                                                tatgCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTAc                                   <  1:984020/56‑1 (MQ=255)
                                                  tgCGTTCCAGTTCGCTGGCATTGCCATAAACCAGCCACGTTGCATAGCCTTCTACCGGAACATAGCGgcgc                   <  1:674065/71‑1 (MQ=255)
                                                              cGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGc                                          >  1:3064999/1‑36 (MQ=255)
                                                              cGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATaaa          >  1:4056693/1‑68 (MQ=255)
                                                                 tGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACAt                           <  1:2196164/48‑1 (MQ=255)
                                                                   gCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTg  <  1:2204043/71‑1 (MQ=255)
                                                                    |                                                                     
TTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTG  >  W3110S.gb/542752‑542889

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: