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.98.8% 410.5 / ‑1.5 85coding (1460/1527 nt)abgTp‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (30/54);  minor base C (0/1);  total (30/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353
                                             |                                       
cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGcc                                   <  2:1570057/50‑1 (MQ=255)
 ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  <  2:4195954/50‑1 (MQ=255)
 ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  >  1:4340260/1‑50 (MQ=255)
 ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  >  1:2228030/1‑50 (MQ=255)
 ccGGAACTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCa                                  <  2:3048486/50‑1 (MQ=255)
   ggACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATa                                <  1:132118/50‑1 (MQ=255)
    gACCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                               <  2:504414/50‑1 (MQ=255)
     acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  1:2941107‑M1/49‑2 (MQ=255)
     acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  1:354‑M1/49‑2 (MQ=255)
     acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  1:87353‑M1/49‑2 (MQ=255)
     acctATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACt                              <  2:2959780‑M1/49‑2 (MQ=255)
     acctATCGGCAGACCGACAAGATACAACGCCAGCACATCGCAGCCATACt                              <  1:248676‑M1/49‑2 (MQ=255)
     aCCTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATac                               >  1:4185098/1‑49 (MQ=38)
      ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  2:1331607‑M1/50‑3 (MQ=255)
      ccTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACta                             <  2:1552850‑M1/50‑3 (MQ=255)
       cTATCGGCAGACCGACAAGATACCGCGCCAGCACATCGCAGCCATACtat                            <  1:594315‑M1/50‑4 (MQ=255)
       cTATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtat                            >  2:417947‑M1/1‑47 (MQ=255)
        tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  2:1138793‑M1/1‑46 (MQ=255)
        tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:3890439‑M1/1‑46 (MQ=255)
        tATCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:1358286‑M1/1‑46 (MQ=255)
          tCGGCAGACCGACAAGATGCCACGCCAGCACATCGCAGCCATACtatgc                          <  1:3044776‑M1/49‑6 (MQ=255)
          tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCGTACtatg                           <  2:880130‑M1/48‑5 (MQ=255)
          tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATGCtatgc                          <  2:1590181‑M1/49‑6 (MQ=255)
          tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgca                         <  1:1924387‑M1/50‑7 (MQ=255)
          tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgca                         <  1:2478189‑M1/50‑7 (MQ=255)
          tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:416712‑M1/1‑44 (MQ=255)
          tCGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           <  1:3304904‑M1/48‑5 (MQ=255)
           cGGCAGACCGACAATATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:3638254‑M1/1‑43 (MQ=255)
           cGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:185135‑M1/1‑43 (MQ=255)
           cGGCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatg                           >  1:3519961‑M1/1‑43 (MQ=255)
            ggCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcag                        <  2:4426809‑M1/49‑8 (MQ=255)
            ggCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcag                        <  2:244872‑M1/49‑8 (MQ=255)
             gCAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaa                      <  1:2418923‑M1/50‑10 (MQ=255)
              cAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     <  2:1721302‑M1/50‑11 (MQ=255)
              cAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     <  2:1389636‑M1/50‑11 (MQ=255)
              cAGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     <  2:4161687‑M1/50‑11 (MQ=255)
               aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     >  2:3381778‑M1/1‑39 (MQ=255)
               aGACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaa                     >  1:2376078‑M1/1‑39 (MQ=255)
                gACCGACAAGATACCACGCCAGCGCATCGCAGCCATACtatgcagaaaca                   <  2:3506634‑M1/50‑13 (MQ=255)
                gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaaca                   <  2:546576‑M1/50‑13 (MQ=255)
                gACCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaaca                   <  2:626872‑M1/50‑13 (MQ=255)
                 aCCGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaa                  >  2:4364328‑M1/1‑37 (MQ=255)
                 aCCGACAAGATACCACGCCAGCACACCGCAGCCATACtatgcagaaacaa                  <  1:2948622‑M1/50‑14 (MQ=255)
                  ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                 >  1:152967‑M1/1‑36 (MQ=255)
                  ccGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaat                 <  2:4221842‑M1/50‑15 (MQ=255)
                   cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                >  1:511022‑M1/1‑35 (MQ=255)
                   cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                <  1:2374560‑M1/50‑16 (MQ=255)
                   cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                <  1:2711154‑M1/50‑16 (MQ=255)
                   cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                <  2:337110‑M1/50‑16 (MQ=255)
                   cGACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatg                >  2:25787‑M1/1‑35 (MQ=255)
                    gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc               <  1:1587107‑M1/50‑17 (MQ=255)
                    gACAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgc               >  2:2509037‑M1/1‑34 (MQ=255)
                     aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              >  2:4028067‑M1/1‑33 (MQ=255)
                     aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              <  1:206947‑M1/50‑18 (MQ=255)
                     aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              <  1:1147814‑M1/50‑18 (MQ=255)
                     aCAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgct              >  2:1855835‑M1/1‑33 (MQ=255)
                      cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg             >  2:2786458‑M1/1‑32 (MQ=255)
                      cAAGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctg             >  1:3181168‑M1/1‑32 (MQ=255)
                       aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            <  2:2033987‑M1/50‑20 (MQ=255)
                       aaGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgt            <  1:3551874‑M1/50‑20 (MQ=255)
                        aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           <  1:4012367‑M1/50‑21 (MQ=255)
                        aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           <  2:4404872‑M1/50‑21 (MQ=255)
                        aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           <  2:2050049‑M1/50‑21 (MQ=255)
                        aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  2:3677660‑M1/1‑30 (MQ=255)
                        aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           <  1:836138‑M1/50‑21 (MQ=255)
                        aGATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           <  2:394245‑M1/50‑21 (MQ=255)
                         gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcg          <  1:502430‑M1/50‑22 (MQ=255)
                         gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  1:1553208‑M1/1‑29 (MQ=255)
                         gATACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtc           >  2:3075773‑M1/1‑29 (MQ=255)
                          aTACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcga         <  2:2744679‑M1/50‑23 (MQ=255)
                           tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat        <  2:3864368‑M1/50‑24 (MQ=255)
                           tACCACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgat        <  1:2553783‑M1/50‑24 (MQ=255)
                             ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg      >  2:3729358‑M1/1‑25 (MQ=255)
                             ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:3061147‑M1/1‑25 (MQ=255)
                             ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatgg      >  2:330503‑M1/1‑25 (MQ=255)
                             ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:228679‑M1/49‑25 (MQ=255)
                             ccACGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatg       <  2:1752962‑M1/49‑25 (MQ=255)
                               aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:1729930‑M1/50‑28 (MQ=255)
                               aCGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:1581281‑M1/50‑28 (MQ=255)
                                cGCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:1023820‑M1/1‑22 (MQ=255)
                                 gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:953704‑M1/50‑30 (MQ=255)
                                 gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  2:4359788‑M1/50‑30 (MQ=255)
                                 gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:1373583‑M1/50‑30 (MQ=255)
                                 gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:665080‑M1/50‑30 (MQ=255)
                                 gCCAGCACATCGCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  2:2338422‑M1/1‑21 (MQ=255)
                                             |                                       
CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353

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