Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,041,983 A→G 7.1% intergenic (+121/+67) 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,9830AG 7.1% 145.6 / 2.7 56intergenic (+121/+67)livJ/yhhKleucine/isoleucine/valine transporter subunit/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/20);  new base G (2/2);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGATTAA  >  W3110S.gb/4041915‑4042053
                                                                    |                                                                      
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)
                      aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                                                >  1:3780058/1‑71 (MQ=255)
                      aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                                                >  1:3537195/1‑71 (MQ=255)
                      aaCGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                                                <  1:3619841/71‑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:1066536/54‑1 (MQ=255)
                                       aCAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCg                                                <  1:979619/54‑1 (MQ=255)
                                        cAAGATCCAAAGAAATCAGTAAATTGCAACACACAt                                                                 >  1:1535589/1‑36 (MQ=255)
                                        cAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAatat                              >  1:102832/1‑71 (MQ=255)
                                        cAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAatat                              >  1:2459931/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:2981288/1‑71 (MQ=255)
                                             tCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACtt                         >  1:2234704/1‑71 (MQ=255)
                                             tCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACtt                         >  1:2985848/1‑71 (MQ=255)
                                             tCCAAAGAAATCAGTAAATTGCAACACACATTTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACttt                        >  1:2905439/1‑71 (MQ=255)
                                                 aaGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCa                                      >  1:3766434/1‑54 (MQ=255)
                                                 aaGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                 >  1:2282761/1‑59 (MQ=255)
                                                 aaGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                 >  1:3112238/1‑59 (MQ=255)
                                                 aaGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                 >  1:2573606/1‑59 (MQ=255)
                                                  aGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAata                               >  1:146085/1‑60 (MQ=255)
                                                  aGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGa                     >  1:3030660/1‑70 (MQ=255)
                                                         aGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAAc             >  1:1235421/1‑71 (MQ=255)
                                                           tAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                 <  1:140610/49‑1 (MQ=255)
                                                           tAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACtt           <  1:28679/71‑1 (MQ=255)
                                                             aaTTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGa         >  1:3722345/1‑71 (MQ=255)
                                                             aaTTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGa         >  1:3391030/1‑71 (MQ=255)
                                                              aTTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGAc        <  1:729529/71‑1 (MQ=255)
                                                                  cAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGATGAtt    <  1:1429910/71‑1 (MQ=255)
                                                                  cAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGAtt    <  1:3251134/71‑1 (MQ=255)
                                                                  cAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGAtt    <  1:659661/71‑1 (MQ=255)
                                                                    acacacATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGATTaa  <  1:2486113/71‑1 (MQ=255)
                                                                    acacacATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGATTaa  <  1:705471/71‑1 (MQ=255)
                                                                    acacacATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGATTaa  <  1:2557014/71‑1 (MQ=255)
                                                                    |                                                                      
GCTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAGATCCAAAGAAATCAGTAAATTGCAACACACATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATCACTTTTGAAATCAAACTTGACGATTAA  >  W3110S.gb/4041915‑4042053

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: