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.90.6% 367.9 / 44.6 96coding (1460/1527 nt)abgTp‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (7/2);  new base . (44/43);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-04

ATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGAC  >  NC_000913/1400266‑1400354
                                                |                                        
aTACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGGAACATcag                                         <  2:3848194/50‑1 (MQ=255)
    ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:258689/1‑50 (MQ=255)
    ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   <  2:3339179/50‑1 (MQ=255)
    ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   <  1:445862/50‑1 (MQ=255)
    ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   <  1:4325640/50‑1 (MQ=255)
    ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:580038/1‑50 (MQ=255)
    ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                   >  2:3895545/1‑50 (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                                >  2:1475188/1‑50 (MQ=255)
       gACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                                <  1:4247199/50‑1 (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                             >  1:2550098‑M1/1‑47 (MQ=255)
          cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             >  2:3514038‑M1/1‑47 (MQ=255)
          cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             <  1:1836900‑M1/50‑4 (MQ=255)
          cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             <  2:1488310‑M1/50‑4 (MQ=255)
          cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                             >  1:2740692‑M1/1‑47 (MQ=255)
           tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  2:725151‑M1/1‑46 (MQ=255)
           tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  2:3663864‑M1/1‑46 (MQ=255)
           tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  1:1513384‑M1/1‑46 (MQ=255)
            aTCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  2:524827‑M1/1‑45 (MQ=255)
            aTCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            >  1:4468845‑M1/1‑45 (MQ=255)
             tCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATaccac                             >  1:1998121/1‑49 (MQ=255)
             tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgca                          <  2:2992384‑M1/50‑7 (MQ=255)
             tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            <  2:4024954‑M1/48‑5 (MQ=255)
             tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                            <  2:1116333‑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                      >  2:3942373‑M1/1‑39 (MQ=255)
                  aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                      >  1:3447053‑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)
                                                |                                        
ATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGAC  >  NC_000913/1400266‑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: