Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,343,621 T→C L85L (TTG→CTG)  melB → melibiose:H(+)/Na(+)/Li(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,343,6210TC97.1% 111.7 / ‑3.9 34L85L (TTG→CTG) melBmelibiose:H(+)/Na(+)/Li(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (10/23);  total (10/24)
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.56e-01

GATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTT  >  NC_000913/4343543‑4343696
                                                                              |                                                                           
gATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTg                                                                           >  2:824867/1‑81 (MQ=255)
  tCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGc                                                                         <  1:824829/81‑1 (MQ=255)
    cGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCaa                                                                       <  2:824829/81‑1 (MQ=255)
          tGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAActct                                                                  <  2:824830/80‑1 (MQ=255)
           gggATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTg                                                                 >  1:824875/1‑80 (MQ=255)
             gatggatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTaa                                                              <  1:824832/81‑1 (MQ=255)
                 gatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCtt                                                          <  1:824833/81‑1 (MQ=255)
                      tAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTtc                                                     <  1:824834/81‑1 (MQ=255)
                               cgcgatcgcgatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCttt                                             <  2:824834/80‑1 (MQ=255)
                               cgcgatcgcgatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCttt                                             <  1:824836/80‑1 (MQ=255)
                                   atcgcgatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTgc                                        >  1:824903/1‑81 (MQ=255)
                                    tcgcgatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTgcg                                       <  2:824836/81‑1 (MQ=255)
                                        gatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATc                                   >  2:824880/1‑81 (MQ=255)
                                        gatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATc                                   <  1:824838/81‑1 (MQ=255)
                                          tGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTg                                 <  2:824838/81‑1 (MQ=255)
                                            gggTAAGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGtt                               <  1:824840/81‑1 (MQ=255)
                                             ggTAAGTTCAAACCCTTGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGttt                              <  2:824840/81‑1 (MQ=255)
                                             ggTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGttt                              <  2:824841/81‑1 (MQ=255)
                                                 aGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAg                          <  1:824843/81‑1 (MQ=255)
                                                 aGTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAg                          <  2:824843/81‑1 (MQ=255)
                                                  gTTCAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAgg                         >  2:824877/1‑81 (MQ=255)
                                                     cAAACCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCTCATCTGTTTGAAGGta                       <  2:824844/80‑1 (MQ=255)
                                                        aCCCTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGtactac                   <  1:824846/81‑1 (MQ=255)
                                                         cccTGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTactact                  <  2:824846/81‑1 (MQ=255)
                                                            tGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTactact                  >  1:824883/1‑78 (MQ=255)
                                                            tGGATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTactact                  >  2:824866/1‑78 (MQ=255)
                                                             ggATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTactact                  >  2:824887/1‑77 (MQ=255)
                                                              gATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGAt             <  2:824847/81‑1 (MQ=255)
                                                              gATCCTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGAt             <  1:824849/81‑1 (MQ=255)
                                                                  cTGATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTc         <  1:824850/81‑1 (MQ=255)
                                                                    gATCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCtt       >  2:824888/1‑81 (MQ=255)
                                                                     aTCGGTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCt        <  2:824850/79‑1 (MQ=255)
                                                                        ggTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCtttgtt   >  1:824887/1‑81 (MQ=255)
                                                                         gTACGCTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCtttgttt  <  2:824851/81‑1 (MQ=255)
                                                                              |                                                                           
GATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTT  >  NC_000913/4343543‑4343696

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: