Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,564,985 T→C 29.8% F21L (TTT→CTT)  yihS → predicted glucosamine isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,564,9850TC29.8% 63.5 / 37.4 47F21L (TTT→CTT) yihSpredicted glucosamine isomerase
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base C (9/5);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCA  >  W3110S.gb/3564917‑3565051
                                                                    |                                                                  
tAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCttt                                                                  <  1:4944489/71‑1 (MQ=255)
tAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGAACGCATCttt                                                                  <  1:429655/71‑1 (MQ=255)
 aaaaGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCctt                                                                  >  1:1589011/1‑70 (MQ=255)
 aaaaGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCctt                                                                  >  1:4695779/1‑70 (MQ=255)
    aGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAtt                                                              <  1:2253590/71‑1 (MQ=255)
     ggAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt                                                             <  1:2982501/71‑1 (MQ=255)
          gAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGt                                                         >  1:1979839/1‑70 (MQ=255)
           aaaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTa                                                        <  1:2834635/70‑1 (MQ=255)
           aaaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTa                                                        >  1:4898519/1‑70 (MQ=255)
             aTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa                                                      >  1:3730606/1‑70 (MQ=255)
              tGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaaaa                                                    >  1:3738647/1‑71 (MQ=255)
              tGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTaaaaa                                                    <  1:3400643/71‑1 (MQ=255)
                gTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAAtt                                                  >  1:3520366/1‑71 (MQ=255)
                 tttAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTc                                                 <  1:5102347/71‑1 (MQ=255)
                     acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCg                                               <  1:3348542/69‑1 (MQ=255)
                     acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGt                                              <  1:4186897/70‑1 (MQ=255)
                     acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGt                                              <  1:4819518/70‑1 (MQ=255)
                        cccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCg                                               <  1:3084037/66‑1 (MQ=255)
                         ccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCt                                                                    <  1:1681279/44‑1 (MQ=255)
                          cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAAt                                                   >  1:4126793/1‑60 (MQ=255)
                          cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTa                                             >  1:850788/1‑66 (MQ=255)
                          cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTa                                             >  1:182346/1‑66 (MQ=255)
                          cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGcc                                        >  1:2747598/1‑71 (MQ=255)
                          cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGcc                                        >  1:2902009/1‑71 (MQ=255)
                          cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGcc                                        <  1:2667358/71‑1 (MQ=255)
                           tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAg                                            >  1:597141/1‑66 (MQ=255)
                             aGCCACAACCGTTGGCTGGAACAGGAAACCGCCCGCATCTTTGATTTTGGTAAAAATTCCGTAGTgac                                        >  1:4527379/1‑66 (MQ=255)
                                 acaACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGt                                 >  1:2428076/1‑71 (MQ=255)
                                 acaACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGt                                 >  1:3086710/1‑71 (MQ=255)
                                  caACCGTTGGCTGCAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtt                                >  1:1938676/1‑71 (MQ=255)
                                    aCCGTTGGCTGGAACAGTAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtttt                              >  1:324013/1‑71 (MQ=255)
                                    aCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtttt                              <  1:3218888/71‑1 (MQ=255)
                                    aCCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtttt                              >  1:3583797/1‑71 (MQ=255)
                                    aCCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtttt                              >  1:3665382/1‑71 (MQ=255)
                                     ccGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTg                             >  1:4579502/1‑71 (MQ=255)
                                             tGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGt                       >  1:1819521/1‑69 (MQ=255)
                                                aaCAGGAAACCGACCGCATCTTTGATTTTGGTaaaaa                                                    >  1:99105/1‑37 (MQ=255)
                                                 aCAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCg                                       <  1:4852940/49‑1 (MQ=255)
                                                    ggAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATa              <  1:4782204/71‑1 (MQ=255)
                                                    ggAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATa              <  1:2292615/71‑1 (MQ=255)
                                                    ggAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATa              <  1:2291360/71‑1 (MQ=255)
                                                     gAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATa              <  1:3163796/70‑1 (MQ=255)
                                                         ccgaccgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAAt               <  1:4911379/65‑1 (MQ=255)
                                                         ccgaccgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAAt               <  1:450566/65‑1 (MQ=255)
                                                             ccgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaa             <  1:678707/63‑1 (MQ=255)
                                                             ccgCATCTTTGATTTTGGTAAAAATTACGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaa             <  1:1126025/63‑1 (MQ=255)
                                                               gCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAg           >  1:2561613/1‑63 (MQ=255)
                                                                cATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCa  <  1:3193835/71‑1 (MQ=255)
                                                                cATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCa  <  1:667402/71‑1 (MQ=255)
                                                                cATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATc   >  1:2391001/1‑70 (MQ=255)
                                                                 aTCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGTCAAATCa  <  1:2067105/70‑1 (MQ=255)
                                                                 aTCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAg           >  1:705403/1‑61 (MQ=255)
                                                                    tttGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTa                     >  1:170701/1‑48 (MQ=255)
                                                                    |                                                                  
TAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCA  >  W3110S.gb/3564917‑3565051

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: