Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,185,626 Δ1 bp intergenic (‑32/‑218) potA ← / → pepT spermidine preferential ABC transporter ATP binding subunit/peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,185,6230G.92.3% 90.9 / 2.2 26intergenic (‑29/‑221)potA/pepTspermidine preferential ABC transporter ATP binding subunit/peptidase T
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (8/16);  total (8/18)
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

CGGTGAAAGCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTGAATCGCCAATC  >  NC_000913/1185547‑1185698
                                                                            |                                                                           
cGGTGAAAGCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCgggtt                                                                         <  2:236217/81‑3 (MQ=255)
 ggTGAAAGCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGttt                                                                        <  2:236218/81‑1 (MQ=255)
  gTGAAAGCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGtttg                                                                       <  2:236219/81‑2 (MQ=255)
     aaaGCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTgg                                                                   >  2:236255/1‑81 (MQ=255)
        gCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGtt                                                                 >  1:236265/1‑80 (MQ=255)
             cTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAaccac                                                           <  2:236220/81‑1 (MQ=255)
                            aaTTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTa                                             <  1:236221/80‑1 (MQ=255)
                              ttttttACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATg                                          <  1:236222/81‑1 (MQ=255)
                              ttttttACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATg                                          <  1:236223/81‑1 (MQ=255)
                                          tCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTa                                 <  2:236224/78‑1 (MQ=255)
                                             cATGTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTa                           <  1:236225/81‑1 (MQ=255)
                                              aTGTAAACGCAACGGAAGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTAc                          <  1:236226/81‑1 (MQ=255)
                                                gTAAACGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACgg                        >  2:236263/1‑81 (MQ=255)
                                                   aaCGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACggta                      >  2:236260/1‑78 (MQ=255)
                                                   aaCGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCaa                     <  2:236228/81‑1 (MQ=255)
                                                   aaCGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCa                      <  1:236227/80‑1 (MQ=255)
                                                   aaCGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCa                      >  2:236253/1‑80 (MQ=255)
                                                   aaCGCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCa                      >  1:236257/1‑80 (MQ=255)
                                                      gCAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAAcc                  <  1:236229/81‑1 (MQ=255)
                                                       cAACGGATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCg                 <  2:236230/81‑1 (MQ=255)
                                                           ggATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTg             <  1:236231/81‑1 (MQ=255)
                                                           ggATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTg             <  2:236232/81‑1 (MQ=255)
                                                           ggATGGCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTg             <  2:236233/81‑1 (MQ=255)
                                                                gCTTACCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTGAATCg        <  2:236234/81‑1 (MQ=255)
                                                                    aCCGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTGAATCGcc      <  1:236235/79‑1 (MQ=255)
                                                                      cGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTGAATCGCCAATc  >  1:236269/1‑81 (MQ=255)
                                                                      cGATGCGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTGAATCGCCAATc  >  1:236258/1‑81 (MQ=255)
                                                                            |                                                                           
CGGTGAAAGCGAACTCGGTTGTTTATTCAATTTTTTACTCTGTCCCATGTAAACGCAACGGATGGCTTACCGATGCGGGGTTTGTGGTTAACCACCTTGGTGACTCTTAATGAGGGCGGTAATTCTACGGCAAACCGCTTGAATCGCCAATC  >  NC_000913/1185547‑1185698

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