Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,400,309 Δ1 bp coding (1465/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,3080A.89.7% 365.7 / 48.3 97coding (1466/1527 nt)abgTp‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (7/3);  new base . (44/43);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

GGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGAC  >  NC_000913/1400261‑1400354
                                               |                                              
ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                              >  2:1575559/1‑48 (MQ=255)
ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaa                                               >  1:567944/1‑49 (MQ=255)
     aTACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGGAACATcag                                         <  2:3848194/50‑1 (MQ=255)
         ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:3895545/1‑50 (MQ=255)
         ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:258689/1‑50 (MQ=255)
         ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   <  2:3339179/50‑1 (MQ=255)
         ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:580038/1‑50 (MQ=255)
         ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   <  1:445862/50‑1 (MQ=255)
         ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   <  1:4325640/50‑1 (MQ=255)
          cGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGcc                                    <  2:5385607/50‑1 (MQ=255)
          cGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:356908/1‑49 (MQ=38)
          cGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCAt                                  <  2:909699/50‑1 (MQ=255)
          cGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCAt                                  <  2:371024/50‑1 (MQ=255)
          cGGACCTATCGGCAGACCGAAAAGATCCCACGCCAGCACATCGCAGCCAt                                  <  2:3517880/50‑1 (MQ=25)
            gACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCa                                   >  1:1477675/1‑49 (MQ=255)
            gACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                                <  1:4247199/50‑1 (MQ=255)
            gACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                                >  2:1475188/1‑50 (MQ=255)
             aCCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                                >  1:393225/1‑49 (MQ=38)
             aCCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                                >  2:3835844/1‑49 (MQ=38)
              ccTATCGGCGGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtg                              <  2:5104904‑M1/50‑3 (MQ=255)
               cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCCGCCATACtat                             >  2:5215530‑M1/1‑47 (MQ=255)
               cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             <  2:1488310‑M1/50‑4 (MQ=255)
               cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             >  2:3514038‑M1/1‑47 (MQ=255)
               cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             >  1:2550098‑M1/1‑47 (MQ=255)
               cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             >  1:2740692‑M1/1‑47 (MQ=255)
               cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             <  1:1836900‑M1/50‑4 (MQ=255)
                tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  2:3663864‑M1/1‑46 (MQ=255)
                tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  1:1513384‑M1/1‑46 (MQ=255)
                tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  2:725151‑M1/1‑46 (MQ=255)
                 aTCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  1:4468845‑M1/1‑45 (MQ=255)
                 aTCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  2:524827‑M1/1‑45 (MQ=255)
                  tCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATaccac                             >  1:1998121/1‑49 (MQ=255)
                  tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgca                          <  2:2992384‑M1/50‑7 (MQ=255)
                  tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            <  2:1116333‑M1/48‑5 (MQ=255)
                  tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            <  2:4024954‑M1/48‑5 (MQ=255)
                  tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            <  2:4152716‑M1/48‑5 (MQ=255)
                  tCGGCAGACCGACAAGATACCACGCCAGCACAGCGCAGCCATACtatgca                          <  2:2276548‑M1/50‑7 (MQ=255)
                   cGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  1:4490945‑M1/1‑43 (MQ=255)
                    ggCAGACCGACAAGGTACCACGCCAGCACATCGCAGCCATACtatgcag                         <  2:3603924‑M1/49‑8 (MQ=255)
                    ggCAGACCGACAAGATACCGCGCCAGCAACATCAGCAGCCATAcaacca                           <  1:1745060/49‑1 (MQ=39)
                     gCAGACCGACAAGATACCTCGCCAGCACATCGCAGCCATACtatgcagaa                       <  1:757698‑M1/50‑10 (MQ=255)
                     gCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaa                       <  2:1817558‑M1/50‑10 (MQ=255)
                      cAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                      <  1:3520184‑M1/50‑11 (MQ=255)
                       aGACTGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                      >  2:235366‑M1/1‑39 (MQ=255)
                       aGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCaaa                         >  1:4240752/1‑48 (MQ=255)
                       aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaac                     <  2:2484021‑M1/50‑12 (MQ=255)
                       aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                      >  1:3447053‑M1/1‑39 (MQ=255)
                       aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                      >  2:3942373‑M1/1‑39 (MQ=255)
                        gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaaca                    <  2:1317150‑M1/50‑13 (MQ=255)
                        gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaac                     >  2:1419719‑M1/1‑38 (MQ=255)
                        gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaac                     >  1:2300411‑M1/1‑38 (MQ=255)
                         aCCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCaaa                         >  2:3698422/1‑46 (MQ=255)
                          ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                  <  1:1017874‑M1/50‑15 (MQ=255)
                          ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                  <  1:3608696‑M1/50‑15 (MQ=255)
                          ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                  >  2:2247009‑M1/1‑36 (MQ=255)
                           cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                 <  1:3375119‑M1/50‑16 (MQ=255)
                           cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                 >  1:656151‑M1/1‑35 (MQ=255)
                            gACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGAt                    >  1:4864632/1‑48 (MQ=255)
                            gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc                >  1:4183500‑M1/1‑34 (MQ=255)
                            gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc                >  1:675691‑M1/1‑34 (MQ=255)
                            gACAAGATACCACGCCAGCACATCGCAGCCAAACtatgcagaaacaatgc                >  2:4174205‑M1/1‑34 (MQ=255)
                             aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct               <  1:545740‑M1/50‑18 (MQ=255)
                             aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct               >  1:1956669‑M1/1‑33 (MQ=255)
                             aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct               <  2:1592199‑M1/50‑18 (MQ=255)
                              cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg              <  1:906676‑M1/50‑19 (MQ=255)
                              cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg              <  2:1114170‑M1/50‑19 (MQ=255)
                              cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg              <  1:4520028‑M1/50‑19 (MQ=255)
                              cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg              >  2:3568443‑M1/1‑32 (MQ=255)
                              cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg              >  2:2874351‑M1/1‑32 (MQ=255)
                               aaGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTg               >  1:5141539/1‑50 (MQ=255)
                               aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt             >  2:3373286‑M1/1‑31 (MQ=255)
                               aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt             <  2:2350279‑M1/50‑20 (MQ=255)
                               aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt             >  2:3673665‑M1/1‑31 (MQ=255)
                                aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc            <  2:4453298‑M1/50‑21 (MQ=255)
                                aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc            <  2:2952782‑M1/50‑21 (MQ=255)
                                aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc            <  1:4089028‑M1/50‑21 (MQ=255)
                                aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc            >  1:5384981‑M1/1‑30 (MQ=255)
                                 gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcg           <  1:2395287‑M1/50‑22 (MQ=255)
                                  aTACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcga          >  2:3374714‑M1/1‑28 (MQ=255)
                                  aTACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcga          <  1:4388648‑M1/50‑23 (MQ=255)
                                   tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat         >  2:547668‑M1/1‑27 (MQ=255)
                                   tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat         <  1:2648631‑M1/50‑24 (MQ=255)
                                   tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat         <  1:4700600‑M1/50‑24 (MQ=255)
                                   tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat         <  1:3940079‑M1/50‑24 (MQ=255)
                                     ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg       >  1:3735321‑M1/1‑25 (MQ=255)
                                     ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg       >  2:2878171‑M1/1‑25 (MQ=255)
                                     ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatg        <  1:5012823‑M1/49‑25 (MQ=255)
                                     ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatg        <  1:2046190‑M1/49‑25 (MQ=255)
                                      cACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggc      >  1:3963888‑M1/1‑24 (MQ=255)
                                       aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgctggcc     >  2:621842‑M1/1‑23 (MQ=255)
                                       aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct     >  1:1985658‑M1/1‑23 (MQ=255)
                                       aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct     >  2:1316140‑M1/1‑23 (MQ=255)
                                       aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct     <  2:797130‑M1/50‑28 (MQ=255)
                                        cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctg    <  2:4273631‑M1/50‑29 (MQ=255)
                                        cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctg    >  2:3815355‑M1/1‑22 (MQ=255)
                                         gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc   <  2:2652381‑M1/50‑30 (MQ=255)
                                         gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc   >  2:212463‑M1/1‑21 (MQ=255)
                                         gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc   >  1:1066729‑M1/1‑21 (MQ=255)
                                             gCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGAc  >  1:4312849/1‑49 (MQ=255)
                                               |                                              
GGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGAC  >  NC_000913/1400261‑1400354

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