Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 897,586 G→A R168H (CGC→CAC)  potI → putrescine ABC transporter membrane subunit PotI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,5860GA97.3% 126.6 / ‑3.7 37R168H (CGC→CAC) potIputrescine ABC transporter membrane subunit PotI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/24);  minor base T (0/1);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

ATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCGCTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTACCGATGATCA  >  NC_000913/897510‑897663
                                                                            |                                                                             
aTCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCg                                                                           <  2:186058/81‑1 (MQ=255)
    ggctggcGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGa                                                                          >  2:186108/1‑78 (MQ=255)
     gctggcGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGa                                                                          <  2:186061/77‑1 (MQ=255)
     gctggcGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATaga                                                                      <  2:186060/81‑1 (MQ=255)
     gctggcGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATa                                                                        <  1:186059/79‑1 (MQ=255)
          cgcATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGaagaag                                                                 <  1:186062/81‑1 (MQ=255)
           gcATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAagc                                                                <  2:186064/81‑1 (MQ=255)
           gcATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAagc                                                                <  2:186063/81‑1 (MQ=255)
            cATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAagca                                                               <  2:186065/81‑1 (MQ=255)
             aTGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAagcag                                                              <  1:186066/81‑1 (MQ=255)
               gTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAagcagc                                                             >  2:186114/1‑80 (MQ=255)
                  aCTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATg                                                         >  2:186113/1‑81 (MQ=255)
                   cTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCGAGATTGCGGAAACTGGATCACTCGATAGAAGAAACAGCGATgg                                                        <  2:186067/81‑1 (MQ=255)
                      tttGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGa                                                       <  2:186068/79‑1 (MQ=255)
                       ttGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGAtct                                                    >  2:186111/1‑81 (MQ=255)
                             ggCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGc                                              <  1:186070/81‑1 (MQ=255)
                             ggCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGc                                              <  2:186069/81‑1 (MQ=255)
                                  aTGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGAc                                           <  1:186071/79‑1 (MQ=255)
                                  aTGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCGCGGTGCGAcg                                          <  2:186072/80‑1 (MQ=255)
                                      ggCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCt                                     >  2:186122/1‑81 (MQ=255)
                                       gCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCt                                     >  2:186131/1‑80 (MQ=255)
                                       gCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTg                                    <  1:186073/81‑1 (MQ=255)
                                              ttATTTCTTCAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAgtgt                             <  1:186074/81‑1 (MQ=255)
                                                  ttcttcAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGttttt                         >  1:186116/1‑81 (MQ=255)
                                                    cttcAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCg                       >  1:186123/1‑81 (MQ=255)
                                                       cAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGt                      <  2:186076/79‑1 (MQ=255)
                                                       cAAGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTc                     <  2:186075/80‑1 (MQ=255)
                                                        aaGATTGCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCa                    <  1:186077/80‑1 (MQ=255)
                                                              gCGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGc              >  1:186115/1‑80 (MQ=255)
                                                               cGGGAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTa            >  2:186125/1‑81 (MQ=255)
                                                                gggAACTGGATCTCTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTAc           <  2:186078/81‑1 (MQ=255)
                                                                gggAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTAc           >  2:186118/1‑81 (MQ=255)
                                                                gggAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTAc           >  1:186117/1‑81 (MQ=255)
                                                                 ggAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTTATTACGCTAcc          <  1:186079/81‑1 (MQ=255)
                                                                 ggAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTAcc          <  2:186080/81‑1 (MQ=255)
                                                                  gAACTGGATCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTACCg         <  2:186081/81‑1 (MQ=255)
                                                                         aTCACTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTACCGATGATCa  <  2:186082/81‑1 (MQ=255)
                                                                            |                                                                             
ATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCGGTCGTTATTTCTTCAAGATTGCGGGAACTGGATCGCTCGATAGAAGAAGCAGCGATGGATCTCGGTGCGACGCCGCTGAAAGTGTTTTTCGTCATTACGCTACCGATGATCA  >  NC_000913/897510‑897663

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