Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,551,633 A→G intergenic (‑337/‑3) gntP ← / → uxuA fructuronate transporter/mannonate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,551,6330AG97.3% 127.8 / NA 37intergenic (‑337/‑3)gntP/uxuAfructuronate transporter/mannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/16);  minor base C (1/0);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCACGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAATCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACCGCGCTGC  >  NC_000913/4551536‑4551729
                                                                                                 |                                                                                                
aCGCCACGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTc                                                                                                <  1:578878/100‑1 (MQ=255)
 cGCCACGGCTGCGGTCTGGTTTTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCa                                                                                               <  2:949374/100‑1 (MQ=255)
 cGCCACGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCa                                                                                               <  2:727117/100‑1 (MQ=255)
    cACGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATgg                                                                                            <  2:852610/100‑1 (MQ=255)
    cACGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATgg                                                                                            <  1:170708/100‑1 (MQ=255)
      cGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGaa                                                                                          >  1:924172/1‑100 (MQ=255)
       ggCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTTTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAAc                                                                                         <  2:696283/100‑1 (MQ=255)
           gCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGAc                                                                                     <  1:1144382/100‑1 (MQ=255)
            cGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGAcc                                                                                    <  1:237295/100‑1 (MQ=255)
              gTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTg                                                                                  >  1:469239/1‑100 (MQ=255)
                 tGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGgcg                                                                               <  1:66108/100‑1 (MQ=255)
                 tGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGgcg                                                                               <  1:536033/100‑1 (MQ=255)
                   gttgttCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCt                                                                             >  1:1092728/1‑100 (MQ=255)
                       ttCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTa                                                                         <  2:1198085/100‑1 (MQ=255)
                            ccGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGcc                                                                    >  2:202411/1‑100 (MQ=255)
                            ccGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAACTCATGGAACAGACCTGGCGCTGGTACGGcc                                                                    >  2:519513/1‑100 (MQ=255)
                                 tACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAAc                                                               >  2:409769/1‑100 (MQ=255)
                                 tACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAAc                                                               >  2:897480/1‑100 (MQ=255)
                                    cTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGAt                                                            >  2:1262280/1‑100 (MQ=255)
                                      aaaCAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATcc                                                          >  2:1119292/1‑100 (MQ=255)
                                            cTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGATCGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTc                                                    <  2:21070/100‑1 (MQ=255)
                                               cAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCttt                                                 <  2:469235/100‑1 (MQ=255)
                                                  ggTTCCGAGTCGCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGc                                              <  2:1261129/100‑1 (MQ=255)
                                                   gTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCt                                             >  2:1162072/1‑100 (MQ=255)
                                                           tctcTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCg                                     >  1:196448/1‑100 (MQ=255)
                                                             tctTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTc                                   >  1:1116322/1‑100 (MQ=255)
                                                             tctTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTc                                   <  1:202413/100‑1 (MQ=255)
                                                                    tgtgGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAggcggg                            >  1:110467/1‑100 (MQ=255)
                                                                      tgGAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAggcgggcg                          <  2:1043951/100‑1 (MQ=255)
                                                                       ggAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGgcgc                         >  2:886856/1‑100 (MQ=255)
                                                                       ggAACCCACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGgcgc                         >  1:1189315/1‑100 (MQ=255)
                                                                             cACTATGTGAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTgg                   >  2:195671/1‑100 (MQ=255)
                                                                                   gtgAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGt             >  2:581372/1‑100 (MQ=255)
                                                                                    tgAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGtt            >  1:619196/1‑100 (MQ=255)
                                                                                     gAAAGAGGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTa           >  2:1065971/1‑100 (MQ=255)
                                                                                        agagGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACcg        >  1:273294/1‑100 (MQ=255)
                                                                                        agagGAAAAGTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACcg        >  1:273423/1‑100 (MQ=255)
                                                                                              aaagTCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACCGCGCTGc  >  2:692723/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACGCCACGGCTGCGGTCTGGTTGTTCATCCGGATACCTAAACAACTCCAGGGTTCCGCGTCTCTTTGCTGTGGAACCCACTATGTGAAAGAGGAAAAATCATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACCGCGCTGC  >  NC_000913/4551536‑4551729

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