Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,343,669 G→T 7.1% E101* (GAA→TAA)  melB → melibiose:sodium symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,343,6690GT 7.1% 182.6 / 3.9 56E101* (GAA→TAA) melBmelibiose:sodium symporter
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base T (2/2);  total (26/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 = 3.71e-01

AGGATCTGGGATGCTATTAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCTTATCCGCG  >  NC_000913/4343522‑4343814
                                                                                                                                                   |                                                                                                                                                 
aGGATCTGGGATGCTATTAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTg                                                                                                                                                   <  5:195388/148‑1 (MQ=255)
        ggATGCTATTAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAgg                                                                                                                                               <  6:278684/144‑1 (MQ=255)
        ggATGCTATTAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGtac                                                                                                                                            >  1:303498/1‑147 (MQ=255)
              tATTAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGAGCATCTGTTTGAAGGTACTACTCAg                                                                                                                                     <  7:187069/148‑1 (MQ=255)
                 tAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGAtt                                                                                                                                  >  3:92244/1‑148 (MQ=255)
                  aaCGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGaa                                                                                                                                                 >  1:315853/1‑132 (MQ=255)
                  aaCGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGaa                                                                                                                                                 <  2:315853/132‑1 (MQ=255)
                   aCGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGt                                                                                                                                >  4:277132/1‑148 (MQ=255)
                      aTCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCtt                                                                                                                             >  3:107916/1‑148 (MQ=255)
                       tCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGt                                                                                                                                              <  4:17466/130‑1 (MQ=255)
                       tCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGt                                                                                                                                              >  3:17466/1‑130 (MQ=255)
                              aTGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGc                                                                                                                      <  2:303498/147‑1 (MQ=255)
                              aTGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGc                                                                                                                      <  8:23504/147‑1 (MQ=255)
                                 ggatggatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGa                                                                                                                  <  1:88073/148‑1 (MQ=255)
                                   atggatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGAcc                                                                                                                >  3:124139/1‑148 (MQ=255)
                                     ggatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGAt                                                                                                                                   <  4:65587/127‑1 (MQ=255)
                                     ggatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGAt                                                                                                                                   >  3:65587/1‑127 (MQ=255)
                                      gatTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTa                                                                                                              >  3:19814/1‑147 (MQ=255)
                                                   aCGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCa                                                                                                >  3:110916/1‑148 (MQ=255)
                                                      cgatcgcgatGGTGTAAGGTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGa                                                                                             <  5:279193/148‑1 (MQ=255)
                                                           gcgatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTac                                                                                        <  4:92244/148‑1 (MQ=255)
                                                             gatGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGTAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTacac                                                                                      <  5:155002/148‑1 (MQ=255)
                                                              atGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTTGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTacac                                                                                      >  6:65697/1‑147 (MQ=255)
                                                                 gggTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCAt                                                                                   <  2:4617/147‑1 (MQ=255)
                                                                   gTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTAt                                                                                >  8:296210/1‑148 (MQ=255)
                                                                     aaGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATg                                                                               >  3:109845/1‑147 (MQ=255)
                                                                     aaGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATg                                                                               >  4:107205/1‑147 (MQ=255)
                                                                      aGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGa                                                                             <  4:109845/148‑1 (MQ=255)
                                                                          cAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCctc                                                                                                       >  1:70332/1‑118 (MQ=255)
                                                                          cAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCctc                                                                                                       <  2:70332/118‑1 (MQ=255)
                                                                               ccTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCtt                                                                    <  3:277132/148‑1 (MQ=255)
                                                                               ccTGGATCCTGATCGGTAAGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACTCCATTATGGATATTCCCtt                                                                    <  8:114312/148‑1 (MQ=255)
                                                                                     tCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTc                                                              <  3:107205/148‑1 (MQ=255)
                                                                                       cTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACAATTATGGATATTCCCTTCTGGTCGc                                                            >  7:186684/1‑148 (MQ=255)
                                                                                             ggTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTc                                                      <  4:107916/148‑1 (MQ=255)
                                                                                             ggTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTc                                                      <  4:19814/148‑1 (MQ=255)
                                                                                               tACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTcc                                                     >  1:107575/1‑147 (MQ=255)
                                                                                                   ttGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCa                                                <  2:132272/148‑1 (MQ=255)
                                                                                                                tAATCTTATTTATCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATa                                    <  4:124139/147‑1 (MQ=255)
                                                                                                                    cTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGc                                                                                                 <  5:245299/82‑1 (MQ=255)
                                                                                                                    cTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGc                                                                                                 >  6:245299/1‑82 (MQ=255)
                                                                                                                       aTTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACg                                <  4:110916/144‑1 (MQ=255)
                                                                                                                                  ttAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACt                 >  4:28690/1‑148 (MQ=255)
                                                                                                                                   tAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGaacaa                   <  2:107575/145‑1 (MQ=255)
                                                                                                                                         gcATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCt          <  7:256199/148‑1 (MQ=255)
                                                                                                                                           aTCTGTTTTAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGa                       <  6:23772/133‑1 (MQ=255)
                                                                                                                                           aTCTGTTTTAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGa                       <  6:19341/133‑1 (MQ=255)
                                                                                                                                           aTCTGTTTTAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGa                       >  5:23772/1‑133 (MQ=255)
                                                                                                                                           aTCTGTTTTAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGa                       >  5:19341/1‑133 (MQ=255)
                                                                                                                                             cTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATa                                    <  2:244137/118‑1 (MQ=255)
                                                                                                                                             cTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATa                                    >  1:244137/1‑118 (MQ=255)
                                                                                                                                             cTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCTTATc      >  2:215097/1‑148 (MQ=255)
                                                                                                                                                tttGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACg                                           <  3:138915/108‑1 (MQ=255)
                                                                                                                                                tttGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACg                                           >  4:138915/1‑108 (MQ=255)
                                                                                                                                                 ttGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCTTAt       >  2:193850/1‑143 (MQ=255)
                                                                                                                                                 ttGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCTTAt       <  1:193850/143‑1 (MQ=255)
                                                                                                                                                 ttGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCTTATCcgcg  >  8:203446/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGATCTGGGATGCTATTAACGATCCGATTATGGGATGGATTGTAAATGCTACGCGATCGCGATGGGGTAAGTTCAAACCCTGGATCCTGATCGGTACGTTGGCAAACTCTGTAATCTTATTTCTCCTCTTTAGTGCGCATCTGTTTGAAGGTACTACTCAGATTGTCTTTGTTTGCGTGACCTACATCCTCTGGGGCATGACTTACACCATTATGGATATTCCCTTCTGGTCGCTGGTTCCAACCATCACGCTCGATAAACGTGAGCGCGAACAACTGGTTCCTTATCCGCG  >  NC_000913/4343522‑4343814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑