Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,400,314 Δ1 bp coding (1460/1527 nt) abgT ← p‑aminobenzoyl glutamate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,400,3140A.99.0% 475.0 / ‑2.0 96coding (1460/1527 nt)abgTp‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (49/46);  minor base C (0/1);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400270‑1400353
                                            |                                       
ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  <  1:5612766/50‑1 (MQ=255)
ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  <  1:8362303/50‑1 (MQ=255)
ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  <  1:7061090/50‑1 (MQ=255)
ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  <  1:4460027/50‑1 (MQ=255)
 cGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  >  1:359003/1‑49 (MQ=38)
  ggACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATa                                <  1:6921358/50‑1 (MQ=255)
  ggACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATa                                <  1:6955792/50‑1 (MQ=255)
  ggACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATa                                <  1:4192937/50‑1 (MQ=255)
  ggACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATa                                >  1:8439757/1‑50 (MQ=255)
  ggACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATa                                >  1:4507532/1‑50 (MQ=255)
   gcccTATCGGCAGACCGACAAGATGCCACGCCAGCACATCGCAGCCATac                               <  1:1306733/48‑1 (MQ=25)
   gACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                               <  1:3685955/50‑1 (MQ=255)
   gACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                               >  1:2989511/1‑50 (MQ=255)
    acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  1:6570235‑M1/49‑2 (MQ=255)
    acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  1:2374586‑M1/49‑2 (MQ=255)
    acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  1:6378805‑M1/49‑2 (MQ=255)
    aCCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                               >  1:8354460/1‑49 (MQ=38)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  1:6683874‑M1/50‑3 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  1:7871804‑M1/50‑3 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  1:7640571‑M1/50‑3 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             >  1:2823167‑M1/1‑48 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  1:2353524‑M1/50‑3 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  1:1195807‑M1/50‑3 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  1:1478150‑M1/50‑3 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             >  1:1717903‑M1/1‑48 (MQ=255)
     ccTATCGGCAGACCGACAAGATACCACGCCAGCACACCGCAGCCATACta                             <  1:8371207‑M1/50‑3 (MQ=255)
      cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                            <  1:800820‑M1/50‑4 (MQ=255)
      cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                            >  1:1709163‑M1/1‑47 (MQ=255)
       tATCTGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:564403‑M1/1‑46 (MQ=255)
       tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:7903895‑M1/1‑46 (MQ=255)
       tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:7292492‑M1/1‑46 (MQ=255)
       tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:955271‑M1/1‑46 (MQ=255)
       tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:3455098‑M1/1‑46 (MQ=255)
       tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCACCCATACtatg                           >  1:4872757‑M1/1‑46 (MQ=255)
         tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgca                         <  1:5573707‑M1/50‑7 (MQ=255)
         tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgc                          <  1:3466558‑M1/49‑6 (MQ=255)
         tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgc                          <  1:5157362‑M1/49‑6 (MQ=255)
         tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           <  1:4096634‑M1/48‑5 (MQ=255)
         tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:8303136‑M1/1‑44 (MQ=255)
          cGGCAGACCGACAAGATACCACGCTAGCACATCGCAGCCATACtatg                           >  1:5112820‑M1/1‑43 (MQ=255)
          cGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:3694191‑M1/1‑43 (MQ=255)
           ggCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcaga                       <  1:7153568‑M1/50‑9 (MQ=255)
           ggCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcaga                       <  1:527902‑M1/50‑9 (MQ=255)
           ggCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcag                        >  1:7712694‑M1/1‑42 (MQ=255)
            gCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaa                      <  1:5964704‑M1/50‑10 (MQ=255)
             cAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     >  1:7938721‑M1/1‑40 (MQ=255)
             cAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     <  1:7208236‑M1/50‑11 (MQ=255)
              aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     >  1:5437435‑M1/1‑39 (MQ=255)
              aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     >  1:1808052‑M1/1‑39 (MQ=255)
               gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaaca                   <  1:1132381‑M1/50‑13 (MQ=255)
               gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaaca                   <  1:1764253‑M1/50‑13 (MQ=255)
               gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaac                    >  1:6176613‑M1/1‑38 (MQ=255)
                aCCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaa                  >  1:5501641‑M1/1‑37 (MQ=255)
                aCCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaa                  <  1:4539608‑M1/50‑14 (MQ=255)
                aCCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaa                  >  1:6371714‑M1/1‑37 (MQ=255)
                aCCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaa                  <  1:357476‑M1/50‑14 (MQ=255)
                 ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                 >  1:8139247‑M1/1‑36 (MQ=255)
                 ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                 <  1:1236618‑M1/50‑15 (MQ=255)
                  cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                >  1:7294064‑M1/1‑35 (MQ=255)
                  cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                >  1:4242184‑M1/1‑35 (MQ=255)
                  cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                <  1:6461665‑M1/50‑16 (MQ=255)
                   gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc               >  1:977811‑M1/1‑34 (MQ=255)
                   gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc               >  1:5064979‑M1/1‑34 (MQ=255)
                   gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc               <  1:6457879‑M1/50‑17 (MQ=255)
                    aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              <  1:7717683‑M1/50‑18 (MQ=255)
                    aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              <  1:8556905‑M1/50‑18 (MQ=255)
                    aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              >  1:3791966‑M1/1‑33 (MQ=255)
                     cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg             >  1:7211206‑M1/1‑32 (MQ=255)
                      aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            >  1:1932386‑M1/1‑31 (MQ=255)
                      aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            >  1:8046834‑M1/1‑31 (MQ=255)
                      aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            >  1:8275288‑M1/1‑31 (MQ=255)
                      aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            >  1:7044678‑M1/1‑31 (MQ=255)
                      aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            >  1:6795176‑M1/1‑31 (MQ=255)
                       aGATACCACGCCAGCACATCGCAGTCATACtatgcagaaacaatgctgtc           <  1:7583968‑M1/50‑21 (MQ=255)
                       aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  1:2673315‑M1/1‑30 (MQ=255)
                       aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           <  1:1791748‑M1/50‑21 (MQ=255)
                       aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  1:1660721‑M1/1‑30 (MQ=255)
                        gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcg          <  1:658536‑M1/50‑22 (MQ=255)
                        gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  1:1171619‑M1/1‑29 (MQ=255)
                        gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  1:5446946‑M1/1‑29 (MQ=255)
                         aTACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcga         <  1:4889959‑M1/50‑23 (MQ=255)
                          tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat        <  1:3768618‑M1/50‑24 (MQ=255)
                          tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat        >  1:8265604‑M1/1‑27 (MQ=255)
                            ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:5568064‑M1/50‑26 (MQ=255)
                            ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:5239990‑M1/1‑25 (MQ=255)
                             cACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:4758792‑M1/1‑24 (MQ=255)
                              aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:3809948‑M1/50‑28 (MQ=255)
                              aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:4290804‑M1/50‑28 (MQ=255)
                              aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:5804664‑M1/1‑23 (MQ=255)
                              aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:6862613‑M1/1‑23 (MQ=255)
                              aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:1460340‑M1/1‑23 (MQ=255)
                               cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:7199055‑M1/50‑29 (MQ=255)
                               cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:5990133‑M1/1‑22 (MQ=255)
                               cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:6043228‑M1/1‑22 (MQ=255)
                               cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgagggctg   <  1:6570321‑M1/50‑29 (MQ=255)
                                gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:3829003‑M1/1‑21 (MQ=255)
                                            |                                       
CCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400270‑1400353

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