Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,769,899 C→T 8.3% S167S (TCG→TCA ydiJ ← predicted FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,769,8990CT 8.3% 187.1 / 5.9 60S167S (TCG→TCAydiJpredicted FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base T (2/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCG  >  W3110S.gb/1769830‑1769969
                                                                     |                                                                      
gagaATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGc                                                                        >  1:2673460/1‑70 (MQ=255)
  gaATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGAt                                                                     <  1:805728/71‑1 (MQ=255)
   aaTATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGTGATc                                                                    >  1:1071559/1‑71 (MQ=255)
    atatCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATcc                                                                   <  1:2180251/71‑1 (MQ=255)
       tCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTg                                                                <  1:3324663/71‑1 (MQ=255)
       tCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTg                                                                <  1:449169/71‑1 (MQ=255)
        cGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTg                                                                <  1:3266797/70‑1 (MQ=255)
           cccccAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCg                                                            <  1:1867459/71‑1 (MQ=255)
                aacaacACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGa                                                          >  1:2891455/1‑68 (MQ=255)
                 acaacaCCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCa                                                      >  1:1239586/1‑71 (MQ=255)
                 acaacaCCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCa                                                      >  1:1597163/1‑71 (MQ=255)
                   aacaCCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATc                                                    >  1:692450/1‑71 (MQ=255)
                   aacaCCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATc                                                    >  1:1206861/1‑71 (MQ=255)
                   aacaCCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATc                                                    >  1:1931190/1‑71 (MQ=255)
                    acacCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCg                                                   >  1:340991/1‑71 (MQ=255)
                      acCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGAc                                                              >  1:469979/1‑58 (MQ=255)
                      acCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGAc                                                              >  1:2782025/1‑58 (MQ=255)
                      acCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGAc                                                              >  1:1510415/1‑58 (MQ=255)
                      acCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGAc                                                              >  1:1993242/1‑58 (MQ=255)
                      acCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGAc                                                              >  1:1086987/1‑58 (MQ=255)
                       ccGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGAt                                                                     <  1:2757816/50‑1 (MQ=255)
                       ccGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGAt                                                                     <  1:3023179/50‑1 (MQ=255)
                         gcgcgTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGAt                                                                     <  1:4016797/48‑1 (MQ=255)
                          cgcgTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTg                                             <  1:3327081/71‑1 (MQ=255)
                          cgcgTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTg                                             <  1:1472656/71‑1 (MQ=255)
                                 gCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGTGATCCCTg                                                                <  1:1584721/45‑1 (MQ=255)
                                 gCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGTGATCCCTg                                                                <  1:1149063/45‑1 (MQ=255)
                                 gCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATcc                                      <  1:3613971/71‑1 (MQ=255)
                                   cAAGTACGTGATCTGACGTTTTGCCATAGACCAGTGATCCCTGAc                                                              >  1:1630611/1‑45 (MQ=255)
                                    aaGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGa                                                          >  1:3121490/1‑48 (MQ=255)
                                    aaGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGa                                                          >  1:1080479/1‑48 (MQ=255)
                                    aaGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGa                                                          >  1:2220472/1‑48 (MQ=255)
                                     aGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATg                                                        >  1:1572399/1‑49 (MQ=255)
                                       tACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCgg                                                  <  1:545912/53‑1 (MQ=255)
                                         cGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCgg                                                           >  1:1399215/1‑42 (MQ=255)
                                         cGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAg                               <  1:3465267/70‑1 (MQ=255)
                                         cGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAg                               >  1:815759/1‑70 (MQ=255)
                                         cGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGc                              <  1:3452457/71‑1 (MQ=255)
                                         cGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGc                              >  1:3833355/1‑71 (MQ=255)
                                         cGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGc                              >  1:3986326/1‑71 (MQ=255)
                                             aTCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTg                                             <  1:3396608/52‑1 (MQ=255)
                                             aTCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTg                                             <  1:669968/52‑1 (MQ=255)
                                              tCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTg                                             <  1:3813886/51‑1 (MQ=255)
                                              tCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTg                                             <  1:3938935/51‑1 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCg                                                            >  1:50115/1‑35 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCg                                                            >  1:360868/1‑35 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCg                                                            >  1:1996248/1‑35 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCg                                                            >  1:1419103/1‑35 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATcc                                      <  1:2218947/57‑1 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATcc                                      <  1:2034726/57‑1 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCAc                                   >  1:823946/1‑60 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGcc                        >  1:536498/1‑71 (MQ=255)
                                               cTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATAATCCCAc                                   >  1:586756/1‑60 (MQ=255)
                                                  aCGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGc                         <  1:3364065/67‑1 (MQ=255)
                                                     ttttGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGat                                           <  1:3304958/46‑1 (MQ=255)
                                                       ttGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGat                                           >  1:1229992/1‑44 (MQ=255)
                                                       ttGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGat                                           >  1:856798/1‑44 (MQ=255)
                                                               gACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGa          >  1:903882/1‑69 (MQ=255)
                                                                 ccAGTGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGt      <  1:2044919/71‑1 (MQ=255)
                                                                   aGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCa                                    >  1:4082640/1‑39 (MQ=255)
                                                                   aGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCa                                    >  1:1359582/1‑39 (MQ=255)
                                                                     cGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCg  >  1:519026/1‑71 (MQ=255)
                                                                     cGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCg  >  1:2161744/1‑71 (MQ=255)
                                                                     |                                                                      
GAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTTTTGCCATAGACCAGCGATCCCTGACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCG  >  W3110S.gb/1769830‑1769969

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: