Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,871,558 T→C 12.6% intergenic (+1384/+221) yibH → / ← yibG hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,871,5580TC12.6% 252.6 / 25.9 100intergenic (+1384/+221)yibH/yibGhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/57);  new base C (9/3);  total (40/60)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 102 total aligned reads displayed.

AGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACG  >  W3110S.gb/3871490‑3871627
                                                                    |                                                                     
aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCAt                                                                       <  1:3369156/69‑1 (MQ=255)
    aCTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTaa                                                                 <  1:2729565/71‑1 (MQ=255)
    aCTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTa                                                                  <  1:3380013/70‑1 (MQ=255)
    aCTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTa                                                                  <  1:4489906/70‑1 (MQ=255)
        gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt                                                              <  1:3232027/70‑1 (MQ=255)
             cTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCt                                                        >  1:1739430/1‑71 (MQ=255)
             cTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCt                                                        <  1:361493/71‑1 (MQ=255)
             cTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCt                                                        >  1:3857693/1‑71 (MQ=255)
             cTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCt                                                        >  1:2083586/1‑71 (MQ=255)
             cTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCt                                                        >  1:4322973/1‑71 (MQ=255)
             cTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCt                                                        <  1:1247859/71‑1 (MQ=255)
                aGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGta                                                     >  1:1246863/1‑71 (MQ=255)
                aGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGta                                                     <  1:4589469/71‑1 (MQ=255)
                aGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGta                                                     >  1:1671764/1‑71 (MQ=255)
                aGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCACACCTAAAGTTATGCTGta                                                     >  1:2656392/1‑71 (MQ=255)
                aGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCACACCTAAAGTTATGCTGta                                                     >  1:3657913/1‑71 (MQ=255)
                   acaTCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATaa                                                  <  1:1001314/71‑1 (MQ=255)
                    catcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATaaa                                                 <  1:4247565/71‑1 (MQ=255)
                     atcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTg                                                       <  1:2532583/64‑1 (MQ=255)
                     atcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTg                                                       <  1:2634782/64‑1 (MQ=255)
                     atcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:3050307/68‑1 (MQ=255)
                     atcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATaa                                                  <  1:421726/69‑1 (MQ=255)
                     atcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAAc                                                <  1:3183516/71‑1 (MQ=255)
                      tcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACa                                               <  1:3363017/71‑1 (MQ=255)
                      tcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACa                                               <  1:4236132/71‑1 (MQ=255)
                      tcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACa                                               <  1:1619970/71‑1 (MQ=255)
                      tcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACa                                               <  1:3099585/71‑1 (MQ=255)
                      tcatAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACa                                               <  1:936733/71‑1 (MQ=255)
                        atAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACa                                               >  1:286930/1‑69 (MQ=255)
                        atAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAgg                                             >  1:3743246/1‑71 (MQ=255)
                        atAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAg                                              <  1:4001751/70‑1 (MQ=255)
                        atAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAg                                              <  1:1689664/70‑1 (MQ=255)
                         tAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt                                                              <  1:2090896/53‑1 (MQ=255)
                         tAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:2097744/64‑1 (MQ=255)
                          aaaaaTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAgg                                             >  1:3035041/1‑69 (MQ=255)
                           aaaaTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt                                                              <  1:1095385/51‑1 (MQ=255)
                           aaaaTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTgg                                          >  1:1470888/1‑71 (MQ=255)
                           aaaaTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTg                                           <  1:4009885/70‑1 (MQ=255)
                              aTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCt                                        <  1:1895865/70‑1 (MQ=255)
                              aTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTa                                       <  1:4140518/71‑1 (MQ=255)
                              aTTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTa                                       <  1:4345454/71‑1 (MQ=255)
                               ttATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt                                                              <  1:2512040/47‑1 (MQ=255)
                               ttATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTaa                                      <  1:3430118/71‑1 (MQ=255)
                               ttATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTa                                       <  1:759212/70‑1 (MQ=255)
                                 atcatcACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:842729/56‑1 (MQ=255)
                                    atcaCATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCa                                 <  1:1766318/71‑1 (MQ=255)
                                     tcaCATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATccc                                  >  1:2165980/1‑69 (MQ=255)
                                       acaTGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATg                                                          <  1:4214559/43‑1 (MQ=255)
                                       acaTGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGt                                                      >  1:3957004/1‑47 (MQ=255)
                                       acaTGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGt                                                      >  1:3962388/1‑47 (MQ=255)
                                       acaTGGTATATCTGAAACACACGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAgg                              >  1:3979073/1‑71 (MQ=255)
                                       acaTGGTATATCTGAAACACACGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAgg                              >  1:2175511/1‑71 (MQ=255)
                                           ggTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGAt                          >  1:2055577/1‑71 (MQ=255)
                                             tataTCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTaa                                      <  1:2025532/57‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCATACCTAAAGTTATGCTg                                                       <  1:2331928/38‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:2482158/42‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:1472706/42‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:3589486/42‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAgggg                            <  1:3025067/65‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAgggg                            <  1:3640499/65‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGAt                      <  1:785658/71‑1 (MQ=255)
                                               tatCTGAAACACACGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGAt                      <  1:769363/71‑1 (MQ=255)
                                                atCTGAAACACACGAAAGCATAGCTAAAGTTATGCTg                                                       <  1:517041/37‑1 (MQ=255)
                                                atCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATaaa                                                 <  1:38111/43‑1 (MQ=255)
                                                 tatGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTc                    >  1:4014651/3‑71 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAgg                                             >  1:2961353/1‑46 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAg                                              >  1:1583138/1‑45 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTg                                           <  1:1148752/48‑1 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTg                                           <  1:346616/48‑1 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGAt                          >  1:4284029/1‑65 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGa                       >  1:474111/1‑68 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTc                    >  1:4344380/1‑71 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTc                    >  1:1939517/1‑71 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTc                    >  1:2125867/1‑71 (MQ=255)
                                                 tCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTc                    >  1:4260921/1‑71 (MQ=255)
                                                 tCTGAAACACACGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCt                                        <  1:1353579/51‑1 (MQ=255)
                                                  cTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTc                    >  1:4217529/1‑70 (MQ=255)
                                                   tGAAACACACGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCTaa                                      >  1:368379/1‑51 (MQ=255)
                                                     aaaCACACGAAAGCATACCTAAAGTTATGCTGtata                                                   <  1:2656880/36‑1 (MQ=255)
                                                      aaCACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGaa               <  1:3940311/71‑1 (MQ=255)
                                                       acacacGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGaaa              >  1:2518617/1‑71 (MQ=255)
                                                       acacacGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGaaa              >  1:3388447/1‑71 (MQ=255)
                                                       acacacGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGaa               <  1:1344406/70‑1 (MQ=255)
                                                          cacGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCg           <  1:1434150/71‑1 (MQ=255)
                                                           acGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATccc                                  <  1:565123/47‑1 (MQ=255)
                                                            cGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCg                        >  1:3229284/1‑56 (MQ=255)
                                                            cGAAAGCACACCTAAAGTTATGCTGTATAAACAgg                                             >  1:2313201/1‑35 (MQ=255)
                                                            cGAAAGCACACCTAAAGTTATGCTGTATAAACAgg                                             >  1:4357262/1‑35 (MQ=255)
                                                            cGAAAGCACACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTcaa                  >  1:2422127/1‑60 (MQ=255)
                                                            cGAAAGCAAACCTAAAGTTATGCTGTATAAACAgg                                             >  1:3879443/1‑35 (MQ=255)
                                                             gAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCaa                                >  1:4320130/1‑47 (MQ=255)
                                                             gAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGt        >  1:2425295/1‑71 (MQ=255)
                                                               aaGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc      <  1:1307912/71‑1 (MQ=255)
                                                               aaGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc      <  1:2470639/71‑1 (MQ=255)
                                                               aaGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc      <  1:4360081/71‑1 (MQ=255)
                                                                aGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGaaa              >  1:3903060/1‑62 (MQ=255)
                                                                aGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCa     <  1:1195849/71‑1 (MQ=255)
                                                                   aTACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCTGGTTCAACg  <  1:2771321/71‑1 (MQ=255)
                                                                   aTACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACg  <  1:3325209/71‑1 (MQ=255)
                                                                    tACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc      <  1:1685418/66‑1 (MQ=255)
                                                                    |                                                                     
AGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACG  >  W3110S.gb/3871490‑3871627

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: