Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 717,029 A→G Y413H (TAC→CAC)  potE ← putrescine:H(+) symporter/putrescine:ornithine antiporter PotE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913717,0290AG100.0% 156.7 / NA 43Y413H (TAC→CAC) potEputrescine:H(+) symporter/putrescine:ornithine antiporter PotE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/18);  total (25/18)

GATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCGACCTTCTGAATAATGACTAATG  >  NC_000913/716913‑717157
                                                                                                                    |                                                                                                                                
gATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAAACAAGGAAGTTCAGGGTGAAACCGTGAAGCATGGCTTCTTcc                                                                                                                  >  2:306815/1‑133 (MQ=255)
 aTCATTCATTCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTccc                                                                                                                 <  1:185500/133‑1 (MQ=255)
 aTCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTccc                                                                                                                 >  1:95671/1‑133 (MQ=255)
           cGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTa                                                                                                       <  2:2849/133‑1 (MQ=255)
                   tttACTTTGTTTTATCAACCGTGTTTGTTTTTGAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCAt                                                                                               <  2:276032/133‑1 (MQ=255)
                     tACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaa                                                                                             <  1:479444/133‑1 (MQ=255)
                     tACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaa                                                                                             <  2:158718/133‑1 (MQ=255)
                     tACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaa                                                                                             <  2:415495/133‑1 (MQ=255)
                     tACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaa                                                                                             <  1:384559/133‑1 (MQ=255)
                       cTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaaaa                                                                                           <  2:95671/133‑1 (MQ=255)
                          tGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGAAGCATGGCATCTACCCAGGATGAGTACGGCGCATAAAAGCt                                                                                        >  2:197861/1‑133 (MQ=255)
                            ttttATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCtat                                                                                      >  2:104588/1‑133 (MQ=255)
                                      cGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGcc                                                                            >  2:443734/1‑133 (MQ=255)
                                          tttATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCcaaca                                                                        >  1:319798/1‑133 (MQ=255)
                                           ttATTTTTCAGTTCATAGCGTGGTGATACCGGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCcaacaa                                                                       <  1:306815/133‑1 (MQ=255)
                                                        cAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGtt                                                          >  2:20287/1‑133 (MQ=255)
                                                                ggtgATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACtt                                                  >  2:345419/1‑133 (MQ=255)
                                                                ggtgATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACAACGCGCCAGCAAAAGCACCAAAGTTTGCAACtt                                                  >  2:232383/1‑133 (MQ=255)
                                                                              aTACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAgggg                                    >  2:32274/1‑133 (MQ=255)
                                                                                aCAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAgggggg                                  >  2:517506/1‑133 (MQ=255)
                                                                                aCAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAgggggg                                  >  1:96300/1‑133 (MQ=255)
                                                                                 cAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGc                                 <  1:20287/133‑1 (MQ=255)
                                                                                    tgtCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGcaca                              >  1:624871/1‑133 (MQ=255)
                                                                                     gtCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACAt                             >  1:112510/1‑133 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAggg                                     <  2:147313/124‑1 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAggg                                     >  1:147313/1‑124 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACAtt                            >  1:143980/1‑133 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACAtt                            >  2:4386/1‑133 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACAtt                            >  2:456773/1‑133 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACAtt                            >  1:41642/1‑133 (MQ=255)
                                                                                      tCCAACCAAGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACAtt                            >  1:316153/1‑133 (MQ=255)
                                                                                       ccaaccaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTa                           >  1:596551/1‑133 (MQ=255)
                                                                                       ccaaccaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTa                           >  2:420618/1‑133 (MQ=255)
                                                                                       ccaaccaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACGTTCGCTTTAGAGGGAGGCACAtaa                           >  2:304160/1‑131 (MQ=255)
                                                                                       ccaaccaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATAGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTa                           >  2:474648/1‑133 (MQ=255)
                                                                                         aaccaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGc                         >  2:483296/1‑133 (MQ=255)
                                                                                          accaaGGAGGGTCACTATGGAATCGTGCAGGATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCg                        <  1:232383/133‑1 (MQ=255)
                                                                                          accaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCg                        <  1:443734/133‑1 (MQ=255)
                                                                                          accaaGGAAGGTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCg                        <  1:4386/133‑1 (MQ=255)
                                                                                               ggaaggTCACGATGGAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCGACCtt                   <  1:517506/133‑1 (MQ=255)
                                                                                                            ggAACCGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCGACCTTCTGAATAATGACt      <  1:319581/133‑1 (MQ=255)
                                                                                                                ccGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAAGAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCGACCTTCTGAATAATGACTAATg  <  2:96300/133‑1 (MQ=255)
                                                                                                                ccGTGCAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCGACCTTCTGAATAATGACTAATg  <  1:480242/133‑1 (MQ=255)
                                                                                                                    |                                                                                                                                
GATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAACAAAAGCAACAAAGTTTGCAACTTTCGCTTTTGAGGGGGGCACATTAGCGACCTTCTGAATAATGACTAATG  >  NC_000913/716913‑717157

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