Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,041,961 C→T 6.5% intergenic (+99/+89) livJ → / ← yhhK leucine/isoleucine/valine transporter subunit/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,041,9610CT 6.5% 183.2 / 2.1 62intergenic (+99/+89)livJ/yhhKleucine/isoleucine/valine transporter subunit/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base T (2/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

AAATTGCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTT  >  W3110S.gb/4041892‑4042031
                                                                     |                                                                      
aaaTTGCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATcc                                                                       <  1:2765304/71‑1 (MQ=255)
aaaTTGCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATcc                                                                       <  1:2740653/71‑1 (MQ=255)
 aaTTGCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCa                                                                      <  1:3431640/71‑1 (MQ=255)
   ttGCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCaaa                                                                    <  1:480475/71‑1 (MQ=255)
     gCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCaaa                                                                    <  1:3540448/69‑1 (MQ=255)
     gCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAg                                                                   <  1:2447318/70‑1 (MQ=255)
         atatCGAGTACGTTGCTTCATGCCGGATTCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAAt                                                               >  1:962036/1‑70 (MQ=255)
         atatCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAAt                                                               <  1:1554038/70‑1 (MQ=255)
              gAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTa                                                          >  1:3042304/1‑70 (MQ=255)
              gAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTa                                                          <  1:3467818/70‑1 (MQ=255)
              gAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTa                                                          <  1:703198/70‑1 (MQ=255)
              gAGTACGTTGCTTCAAGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTa                                                          <  1:2255029/70‑1 (MQ=255)
                       gCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAaca                                                >  1:563416/1‑71 (MQ=255)
                       gCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAaca                                                >  1:1911924/1‑71 (MQ=255)
                               ccGGATGCGGCGTAAACGCCTTATCCGGCGTACAAGATCCAAAGAAATCAGTATATTGCAACACACATTg                                         <  1:1739856/70‑1 (MQ=37)
                               ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATTCAAAGAAATCAGTAAATTGCAGCACACATTg                                         <  1:3488744/70‑1 (MQ=255)
                               ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATTCAAAGAAATCAGTAAATTGCAGCACACATTg                                         <  1:549171/70‑1 (MQ=255)
                               ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACAtt                                          >  1:2218351/1‑69 (MQ=255)
                               ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTg                                         >  1:1417776/1‑70 (MQ=255)
                                     gcggcgTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCt                                  >  1:1069109/1‑71 (MQ=255)
                                       ggcgTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAacac                                               <  1:1109863/56‑1 (MQ=255)
                                        gcgTAAACGCCTTATCCGGCCTACAAGATTCAAAGAAATCAGTAAATTGCAGCACACATTGTa                                       >  1:2380531/1‑63 (MQ=37)
                                        gcgTAAACGCCTTATCCGGCCTACAAGATTCAAAGAAATCAGTAAATTGCAGCACACATTGTa                                       >  1:3669056/1‑63 (MQ=37)
                                           tAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGt                                                           <  1:1815746/40‑1 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTaaa                                                        >  1:716220/1‑41 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTaaa                                                        >  1:1321158/1‑41 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTaaa                                                        >  1:1413622/1‑41 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCa                                                   >  1:1800547/1‑46 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCa                                                   >  1:3664427/1‑46 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                          >  1:3537195/1‑71 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                          <  1:3619841/71‑1 (MQ=255)
                                             aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                          >  1:3780058/1‑71 (MQ=255)
                                              aCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAg                                                            <  1:915647/36‑1 (MQ=255)
                                              aCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAg                                                            <  1:1483299/36‑1 (MQ=255)
                                              aCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGt                         <  1:271361/71‑1 (MQ=255)
                                               cGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAAcacaca                                            >  1:1712232/1‑51 (MQ=255)
                                                  cTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAgcg                     >  1:3195828/1‑71 (MQ=255)
                                                   ttATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAgcg                     >  1:3064916/1‑70 (MQ=255)
                                                     aTCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAAcacaca                                            >  1:686924/1‑45 (MQ=255)
                                                     aTCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:702827/1‑71 (MQ=255)
                                                      tCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTg                                         <  1:3776046/47‑1 (MQ=255)
                                                      tCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCAt                  >  1:2593400/1‑70 (MQ=255)
                                                        cGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTa                                       >  1:897941/1‑47 (MQ=255)
                                                          gCCTACAAGATGCAAAGAAATCAGTAAATTGCAACACACAt                                           <  1:1917016/41‑1 (MQ=255)
                                                          gCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGc             <  1:1090155/71‑1 (MQ=255)
                                                           ccTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCa            <  1:2195639/71‑1 (MQ=255)
                                                            cTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa           <  1:2793776/71‑1 (MQ=255)
                                                             tACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATaa                             <  1:2574454/52‑1 (MQ=255)
                                                             tACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATaa                             <  1:2831936/52‑1 (MQ=255)
                                                              aCAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAg                                      <  1:124760/42‑1 (MQ=255)
                                                              aCAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                          <  1:979619/54‑1 (MQ=255)
                                                              aCAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                          <  1:1066536/54‑1 (MQ=255)
                                                               cAAGATCCAAAGAAATCAGTAAATTGCAACACACAt                                           >  1:1535589/1‑36 (MQ=255)
                                                               cAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAatat        >  1:2459931/1‑71 (MQ=255)
                                                               cAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAatat        >  1:102832/1‑71 (MQ=255)
                                                                aaGATCCAAAGAAATCAGTAAATTGCAACACACATTg                                         <  1:3787279/37‑1 (MQ=255)
                                                                aaGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAg                       >  1:3564275/1‑55 (MQ=255)
                                                                aaGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAatat        >  1:634801/1‑70 (MQ=255)
                                                                aaGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAata         >  1:1698581/1‑69 (MQ=255)
                                                                 aGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCa      >  1:923871/1‑71 (MQ=255)
                                                                    tCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACtt   >  1:2985848/1‑71 (MQ=255)
                                                                    tCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACtt   >  1:2981288/1‑71 (MQ=255)
                                                                    tCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACtt   >  1:2234704/1‑71 (MQ=255)
                                                                    tCCAAAGAAATCAGTAAATTGCAACACACATTTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACttt  >  1:2905439/1‑71 (MQ=255)
                                                                     |                                                                      
AAATTGCTAATATCGAGTACGTTGCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTT  >  W3110S.gb/4041892‑4042031

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: