Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,193,567 G→A 9.0% R226R (CGC→CGT yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,193,5670GA 9.0% 126.3 / 2.1 44R226R (CGC→CGTyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (2/2);  total (23/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 = 1.00e+00

ACAACCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCTG  >  W3110S.gb/2193503‑2193634
                                                                |                                                                   
acaacCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGc                                                               >  1:4322756/1‑71 (MQ=255)
acaacCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGc                                                               >  1:2040894/1‑71 (MQ=255)
           ccTCTGAAGCGAACAGAATCAAAAATATCGGTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAg                                                    >  1:2994282/1‑71 (MQ=255)
               tGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGtt                                                  <  1:2133050/69‑1 (MQ=255)
               tGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGtt                                                  <  1:359876/69‑1 (MQ=255)
               tGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTgg                                                >  1:3807243/1‑71 (MQ=255)
               tGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTgg                                                >  1:3552432/1‑71 (MQ=255)
                   gCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCaa                                            <  1:2846179/71‑1 (MQ=255)
                      aaCAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATc                                         <  1:2225548/71‑1 (MQ=255)
                        cAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCa                                       >  1:893708/1‑71 (MQ=255)
                        cAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCa                                       <  1:3851322/71‑1 (MQ=255)
                        cAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCa                                       >  1:1095168/1‑71 (MQ=255)
                          gAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCaa                                            <  1:1978056/64‑1 (MQ=255)
                          gAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAAt                                          >  1:1415107/1‑66 (MQ=255)
                          gAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACg                                     >  1:2203150/1‑71 (MQ=255)
                             tCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTc                                  >  1:3698127/1‑71 (MQ=255)
                               aaaaaTATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCa                                       >  1:2334106/1‑64 (MQ=255)
                                  aaTATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGAt                             <  1:2682629/71‑1 (MQ=255)
                                   atatCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATcc                                        <  1:1837626/59‑1 (MQ=255)
                                   atatCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATa                            <  1:3027937/71‑1 (MQ=255)
                                        gCTTGATGTGTACATATCACCCACGCGAAGCGTGCCgaga                                                      <  1:1342807/40‑1 (MQ=255)
                                            gATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTa                     <  1:3030276/69‑1 (MQ=255)
                                            gATGTGTACATATCACACACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTAt                    <  1:1025453/70‑1 (MQ=255)
                                             aTGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTg                  >  1:1026299/1‑71 (MQ=255)
                                             aTGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTg                  >  1:2643627/1‑71 (MQ=255)
                                             aTGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTg                  >  1:3625353/1‑71 (MQ=255)
                                                tgtACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCa                                <  1:3624642/54‑1 (MQ=255)
                                                  tACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCa                                <  1:3021163/52‑1 (MQ=255)
                                                  tACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCa                                <  1:3762931/52‑1 (MQ=255)
                                                  tACATATCACCCACGCGAAGCGTGACGAGAAGTTGGGCAAATCCACGTTCCAGAt                             <  1:2586553/55‑1 (MQ=255)
                                                    cATATCACCCACGCGAAGCGTGCCTAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTc            <  1:398224/70‑1 (MQ=255)
                                                    cATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTc            >  1:1886528/1‑70 (MQ=255)
                                                    cATATCACCCACACGAAGCGTGCCGAGAAGTTGGGCAAAt                                          >  1:1453370/1‑40 (MQ=255)
                                                     atatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCca          >  1:562496/1‑71 (MQ=255)
                                                      tatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAg                               >  1:4396283/1‑49 (MQ=255)
                                                      tatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGa                              >  1:198182/1‑50 (MQ=255)
                                                      tatCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGa                              >  1:715892/1‑50 (MQ=255)
                                                      tatCACCCACACGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTAt                    <  1:2491543/60‑1 (MQ=255)
                                                      tatCACCCACACGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTAt                    <  1:2802266/60‑1 (MQ=255)
                                                        tCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCcaca        >  1:2423509/1‑70 (MQ=255)
                                                        tCACCCACACGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCcaca        >  1:2119397/1‑70 (MQ=255)
                                                            ccACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAgg                        <  1:382105/50‑1 (MQ=255)
                                                             cACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCTg  <  1:4286558/71‑1 (MQ=255)
                                                             cACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCTg  >  1:184440/1‑71 (MQ=255)
                                                                gcgAAGCGTGCCGAGAAGTTGGGCAAATCCTCGTTCCAGATACAGGGTATTGCTTTTCCACAcc      <  1:1645123/64‑1 (MQ=255)
                                                                gcgAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACAcc      <  1:2607858/64‑1 (MQ=255)
                                                                |                                                                   
ACAACCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTGATGTGTACATATCACCCACGCGAAGCGTGCCGAGAAGTTGGGCAAATCCACGTTCCAGATACAGGGTATTGCTTTTCCACACCCCTG  >  W3110S.gb/2193503‑2193634

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