Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,185,245 (T)6→5 coding (350/1137 nt) potA ← spermidine preferential ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,185,2400T.96.2% 99.3 / ‑1.3 26coding (355/1137 nt)potAspermidine preferential ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (9/16);  minor base G (1/0);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAA  >  NC_000913/1185162‑1185314
                                                                              |                                                                          
gTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGttt                                                                          <  2:236141/81‑1 (MQ=255)
       ttGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATg                                                                  >  2:236193/1‑81 (MQ=255)
         gAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATgc                                                                 <  1:236143/80‑1 (MQ=255)
               aGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAAccc                                                          <  2:236144/81‑1 (MQ=255)
                      cAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAATGcc                                                   <  2:236145/81‑1 (MQ=255)
                         cTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCcaca                                                <  2:236146/81‑1 (MQ=255)
                            cACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACAttt                                             <  2:236147/81‑1 (MQ=255)
                             aCCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACAtttt                                            >  1:236192/1‑81 (MQ=255)
                                aTCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCg                                          >  2:236205/1‑80 (MQ=255)
                                   cGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAaca                                      <  2:236148/81‑1 (MQ=255)
                                    gCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAacac                                     >  2:236187/1‑81 (MQ=255)
                                    gCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAacac                                     <  2:236149/81‑1 (MQ=255)
                                     cAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACg                                    <  1:236150/81‑1 (MQ=255)
                                       gggCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGt                                  <  1:236151/81‑1 (MQ=255)
                                          cTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCAt                               <  1:236152/81‑1 (MQ=255)
                                                      ggggCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAg                   <  2:236153/81‑1 (MQ=255)
                                                       gggCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGt                  >  1:236191/1‑81 (MQ=255)
                                                          cGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCg               >  2:236182/1‑81 (MQ=255)
                                                                 tCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCtt         <  1:236155/80‑1 (MQ=255)
                                                                 tCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCtt         <  1:236154/80‑1 (MQ=255)
                                                                  cagcagCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTg       >  2:236197/1‑81 (MQ=255)
                                                                   agcagcGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTgg      <  1:236156/81‑1 (MQ=255)
                                                                   agcagcGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTg       <  1:236157/80‑1 (MQ=255)
                                                                     cagcGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaa    <  1:460615/81‑1 (MQ=255)
                                                                      agcGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaa    >  1:236213/1‑80 (MQ=255)
                                                                      agcGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaa    >  2:236208/1‑80 (MQ=255)
                                                                        gcggggGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaaaa  >  1:236186/3‑81 (MQ=255)
                                                                              |                                                                          
GTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAA  >  NC_000913/1185162‑1185314

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