Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,320,162 T→C 5.5% Y436C (TAT→TGT)  yhbX ← predicted hydrolase, inner membrane

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,320,1620TC 5.5% 212.1 / 3.1 73Y436C (TAT→TGT) yhbXpredicted hydrolase, inner membrane
Reads supporting (aligned to +/- strand):  ref base T (41/28);  new base C (2/2);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCG  >  W3110S.gb/3320093‑3320229
                                                                     |                                                                   
aCGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATa                                                                    <  1:2753244/71‑1 (MQ=255)
       tataCCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGttc                                                             >  1:1276558/1‑71 (MQ=255)
         taCCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGttctt                                                           <  1:3915013/71‑1 (MQ=255)
         taCCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAACAGACGttctt                                                           >  1:955902/1‑71 (MQ=255)
         taCCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAACAGACGttct                                                            >  1:571866/1‑70 (MQ=255)
            cAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGttcttctt                                                        >  1:3790409/1‑71 (MQ=255)
               aTAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGt                                                     <  1:292854/71‑1 (MQ=255)
               aTAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGt                                                     <  1:707029/71‑1 (MQ=255)
                  aaCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAg                                                                   <  1:3866863/54‑1 (MQ=255)
                   aCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAg                                                                   <  1:5069272/53‑1 (MQ=255)
                   aCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTa                                                    <  1:2663765/68‑1 (MQ=255)
                    cATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGt                                                >  1:5073348/1‑71 (MQ=255)
                     aTCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAgg                                                  >  1:4991760/1‑68 (MQ=255)
                     aTCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGt                                                >  1:1683142/1‑70 (MQ=255)
                      tCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCa                                              <  1:699899/71‑1 (MQ=255)
                       cGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCAc                                             >  1:3042558/1‑71 (MQ=255)
                        gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACg                                            >  1:4655900/1‑71 (MQ=255)
                            cATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATACACGTTCTTCTTCGTAGGGTCACGTTcc                                        <  1:3642876/71‑1 (MQ=255)
                             aTGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATCAAAATAGACGTTCTTCTTCGTAGGGTCACGTTc                                         >  1:3442427/1‑69 (MQ=255)
                               gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGttc                                                             >  1:5075936/1‑47 (MQ=255)
                               gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCAc                                             <  1:666266/63‑1 (MQ=255)
                               gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAACAGACGtt                                                              <  1:923689/46‑1 (MQ=255)
                                atatGCCTTCTGGCTAGCCTCCCTGCCTCCATGAAAACAGACGtt                                                              <  1:1806879/45‑1 (MQ=255)
                                atatGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGc                                    >  1:5044481/1‑71 (MQ=255)
                                   tGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATa                                                                    <  1:647665/36‑1 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGa                                                                  >  1:572489/1‑36 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTa                                                    >  1:3941065/1‑50 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTgg                               >  1:3144901/1‑71 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTgg                               >  1:4995752/1‑71 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTgg                               >  1:4841760/1‑71 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTg                                >  1:4167270/1‑70 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTg                                >  1:1248734/1‑70 (MQ=255)
                                     ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTg                                >  1:4980478/1‑70 (MQ=255)
                                       tgctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTa                                                    <  1:3101780/48‑1 (MQ=255)
                                       tgctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTa                                                    <  1:3219025/48‑1 (MQ=255)
                                        gctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGt                                                               <  1:841313/36‑1 (MQ=255)
                                         ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCgg                           >  1:4362502/1‑71 (MQ=255)
                                         ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCgg                           >  1:2052587/1‑71 (MQ=255)
                                         ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCgg                           >  1:1451015/1‑71 (MQ=255)
                                            gCTAGCCTCCCTGCCTCCATGAAAATAGGCGTTCTTCTTCGTAGGGGCACGTTCCAGGCCGTGGTCGGCaa                        <  1:589306/71‑1 (MQ=255)
                                              tAGCCTCCCTGCCTCCATGAAAATAGACGttcttcttc                                                       >  1:4843630/1‑38 (MQ=255)
                                              tAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCg                                                      >  1:1150452/1‑39 (MQ=255)
                                                    cccTGCCTCCATGAAAATAGACGTTCTTCTTCGTa                                                    <  1:102176/35‑1 (MQ=255)
                                                     ccTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGt                              >  1:4314657/1‑56 (MQ=255)
                                                       tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGt                                           <  1:332140/41‑1 (MQ=255)
                                                       tGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGt                                           <  1:2768813/41‑1 (MQ=255)
                                                         ccTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGtt                                          >  1:2621389/1‑40 (MQ=255)
                                                         ccTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTc                                         >  1:3538140/1‑41 (MQ=255)
                                                          cTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGt                                           <  1:1424193/38‑1 (MQ=255)
                                                           tCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCaaaa                      >  1:2614507/1‑58 (MQ=255)
                                                            ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAg                                      >  1:1454838/1‑41 (MQ=255)
                                                            ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTgg                               <  1:3466198/48‑1 (MQ=255)
                                                            ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc        >  1:3470568/1‑71 (MQ=255)
                                                             cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTgg                               <  1:1862263/47‑1 (MQ=255)
                                                             cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcg       >  1:1291186/1‑71 (MQ=255)
                                                             cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcg       >  1:4954883/1‑71 (MQ=255)
                                                              aTGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAg                                      >  1:1516235/1‑39 (MQ=255)
                                                              aTGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAg                                      >  1:1941925/1‑39 (MQ=255)
                                                               tGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGt                              <  1:3764636/46‑1 (MQ=255)
                                                               tGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGt                              <  1:2849782/46‑1 (MQ=255)
                                                               tGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACa                  >  1:3245199/1‑58 (MQ=255)
                                                               tGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACa                  >  1:3925063/1‑58 (MQ=255)
                                                               tGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGcgcg     >  1:4948956/1‑71 (MQ=255)
                                                                gAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGt                              <  1:3638958/45‑1 (MQ=255)
                                                                gAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACAt                 <  1:5048605/58‑1 (MQ=255)
                                                                gAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACAt                 <  1:3053552/58‑1 (MQ=255)
                                                                gAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcgg    >  1:448349/1‑71 (MQ=255)
                                                                 aaaaTAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc   >  1:2711822/1‑71 (MQ=255)
                                                                 aaaaTAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc   >  1:3107781/1‑71 (MQ=255)
                                                                  aaaTAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGa               >  1:2069673/1‑58 (MQ=255)
                                                                  aaaTAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggcg  <  1:4520330/71‑1 (MQ=255)
                                                                  aaaTAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggcg  <  1:1221912/71‑1 (MQ=255)
                                                                   aaTAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGt                              >  1:1047622/1‑42 (MQ=255)
                                                                     |                                                                   
ACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCG  >  W3110S.gb/3320093‑3320229

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: