Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,353,056	T→A	P191P (CCA→CCT)	sapD ←	putrescine ABC exporter ATP binding protein SapD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,353,056	0	T	A	100.0%	148.0 / NA	47	P191P (CCA→CCT)	sapD	putrescine ABC exporter ATP binding protein SapD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (24/23);  total (24/23)

TTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGT  >  NC_000913/1352921‑1353186                                                                                                                                        |                                                                                                                                   ttGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCt                                                                                                                                >  2:279891/1‑140 (MQ=255)          gTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGt                                                                                                                       >  1:47086/1‑140 (MQ=255)                cACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCg                                                                                                                 >  2:171978/1‑140 (MQ=255)                cACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCg                                                                                                                 >  2:621320/1‑140 (MQ=255)                   gTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCgg                                                                                                                >  1:677722/1‑138 (MQ=255)                   gTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCgg                                                                                                                <  2:677722/138‑1 (MQ=255)                   gTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGtt                                                                                                              >  1:803390/1‑140 (MQ=255)                    ttAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTc                                                                                                             >  1:305492/1‑140 (MQ=255)                                     ttGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGAc                                                                                            >  1:839234/1‑140 (MQ=255)                                     ttGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGAc                                                                                            >  1:235662/1‑140 (MQ=255)                                            aGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTg                                                                                     >  1:696498/1‑140 (MQ=255)                                            aGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTg                                                                                     >  2:329856/1‑140 (MQ=255)                                                  tGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGc                                                                               <  1:171978/140‑1 (MQ=255)                                                  tGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGc                                                                               <  1:621320/140‑1 (MQ=255)                                                  tGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGc                                                                               >  1:753030/1‑140 (MQ=255)                                                  tGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGc                                                                               <  2:305492/140‑1 (MQ=255)                                                  tGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGc                                                                               >  2:147330/1‑140 (MQ=255)                                                           tGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAAt                                                                      >  2:836582/1‑140 (MQ=255)                                                           tGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAAt                                                                      >  1:845288/1‑140 (MQ=255)                                                                 aTAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTAGCCAg                                                                            <  1:485585/128‑1 (MQ=255)                                                                 aTAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTAGCCAg                                                                            >  2:485585/1‑128 (MQ=255)                                                                   aaGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCa                                                                        <  1:136446/130‑1 (MQ=255)                                                                   aaGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCa                                                                        >  2:136446/1‑130 (MQ=255)                                                                      cAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTAtcatca                                                           >  1:573251/1‑140 (MQ=255)                                                                          aTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATAACttt                                                       <  2:696498/140‑1 (MQ=255)                                                                           tAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTc                                                      <  2:573251/140‑1 (MQ=255)                                                                           tAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTc                                                      <  2:839234/140‑1 (MQ=255)                                                                           tAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTc                                                      <  2:235662/140‑1 (MQ=255)                                                                                gTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTatcat                                                             <  2:87862/128‑1 (MQ=255)                                                                                gTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTatcat                                                             >  1:87862/1‑128 (MQ=255)                                                                                   ttGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTc                                              >  1:834050/1‑140 (MQ=255)                                                                                      cTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCAcc                                           >  1:278882/1‑140 (MQ=255)                                                                                       tGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCt                                          <  1:329856/140‑1 (MQ=255)                                                                                        gTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCtt                                         >  1:831147/1‑140 (MQ=255)                                                                                          tttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCg                                       <  2:47086/140‑1 (MQ=255)                                                                                           ttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCgg                                      >  2:540704/1‑140 (MQ=255)                                                                                             ggTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTc                                    <  2:845288/140‑1 (MQ=255)                                                                                                tGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAAc                                 <  1:671444/140‑1 (MQ=255)                                                                                                  agaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTc                               <  1:819918/140‑1 (MQ=255)                                                                                                  agaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTc                               <  2:831147/140‑1 (MQ=255)                                                                                                         tCAGCAGGCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTTCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGa                        >  2:31972/1‑140 (MQ=255)                                                                                                               ggCGAAAGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTg                  <  1:342477/140‑1 (MQ=255)                                                                                                                   aaaGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgca              <  1:147330/140‑1 (MQ=255)                                                                                                                     aGATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatc            <  2:753030/140‑1 (MQ=255)                                                                                                                      gATTTGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatcg           <  1:540704/140‑1 (MQ=255)                                                                                                                        tttGCGCCTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatcgca         <  2:834050/140‑1 (MQ=255)                                                                                                                               cTGGGTTGTAGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTtgt  <  1:836582/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   TTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGT  >  NC_000913/1352921‑1353186

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