Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,521	C→T	A343V (GCG→GTG)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,521	0	C	T	97.7%	140.6 / NA	44	A343V (GCG→GTG)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/26);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAA  >  NC_000913/3220387‑3220654                                                                                                                                       |                                                                                                                                      tGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTgg                                                                                                                                  >  1:330144/1‑140 (MQ=255)  gCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGc                                                                                                                                 <  2:382276/140‑1 (MQ=255)  gCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGc                                                                                                                                 >  1:314829/1‑140 (MQ=255)  gCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGc                                                                                                                                 <  2:458376/140‑1 (MQ=255)        gcgcTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCAt                                                                                                                           <  1:308272/140‑1 (MQ=255)           cTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCaa                                                                                                                        >  2:116554/1‑140 (MQ=255)             cggcggcggCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAAtg                                                                                                                      >  2:566635/1‑140 (MQ=255)                   cggcgTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGc                                                                                                                <  1:566635/140‑1 (MQ=255)                          aTGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCa                                                                                                         <  2:519552/140‑1 (MQ=255)                                 ttGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTa                                                                                                  <  2:330144/140‑1 (MQ=255)                                 ttGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTa                                                                                                  <  2:364197/140‑1 (MQ=255)                                 ttGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTa                                                                                                  <  1:595879/140‑1 (MQ=255)                                    gcgcGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCg                                                                                               >  2:121326/1‑140 (MQ=255)                                     cgcgACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCgg                                                                                              >  2:496580/1‑140 (MQ=255)                                     cgcgACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCgg                                                                                              >  2:214312/1‑140 (MQ=255)                                                      aaGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTCACCGGCTCAGGCTGGGCaaaaa                                                                             >  2:417305/1‑140 (MQ=255)                                                       agagGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAg                                                                            >  2:289940/1‑140 (MQ=255)                                                          ggTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCg                                                                         <  2:828893/140‑1 (MQ=255)                                                           gtgtTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGAACATCAATTTGTTGCTGGAGCAGAACTTATCGGCTCAGGCTGGGCAAAAAGGCGa                                                                        >  2:611042/1‑140 (MQ=255)                                                                       cTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGa                                                            <  1:496580/140‑1 (MQ=255)                                                                       cTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGa                                                            <  2:646888/140‑1 (MQ=255)                                                                       cTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGa                                                            <  1:214312/140‑1 (MQ=255)                                                                          ttCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACgg                                                         >  1:540254/1‑140 (MQ=255)                                                                          ttCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACgg                                                         >  1:733754/1‑140 (MQ=255)                                                                           tCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGt                                                        <  2:524713/140‑1 (MQ=255)                                                                                    cATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAACGTGTTGCTGGAGCAGAACTTACCGGGTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCaa                                               >  1:574555/1‑140 (MQ=255)                                                                                       tCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTg                                            <  1:116554/140‑1 (MQ=255)                                                                                       tCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTg                                            >  1:253573/1‑140 (MQ=255)                                                                                        ccTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAACGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTgg                                           <  2:574555/140‑1 (MQ=255)                                                                                             aTACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCgg                                      >  2:468825/1‑140 (MQ=255)                                                                                              tACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCggg                                     <  2:37915/140‑1 (MQ=255)                                                                                                 caccaccTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAat                                  >  1:101487/1‑140 (MQ=255)                                                                                                 caccaccTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAat                                  >  1:410348/1‑140 (MQ=255)                                                                                                   ccaccTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                >  1:372215/1‑140 (MQ=255)                                                                                                   ccaccTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCCGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                >  1:569869/1‑140 (MQ=255)                                                                                                     accTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATcc                              <  2:372215/140‑1 (MQ=255)                                                                                                     accTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATcc                              <  2:253573/140‑1 (MQ=255)                                                                                                       cTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCg                            <  2:733754/140‑1 (MQ=255)                                                                                                         ttGGCGGCAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGAt                          <  2:540254/140‑1 (MQ=255)                                                                                                           ggcggcAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCt                        <  1:611042/140‑1 (MQ=255)                                                                                                           ggcggcAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCt                        <  1:34237/140‑1 (MQ=255)                                                                                                           ggcggcAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCt                        <  2:314829/140‑1 (MQ=255)                                                                                                             cggcAACCCGCTGGCCTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTgg                      <  2:410348/140‑1 (MQ=255)                                                                                                                            cTGTGCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTg       <  2:750685/140‑1 (MQ=255)                                                                                                                              gtgCGGCGGTGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGt     <  2:101487/140‑1 (MQ=255)                                                                                                                                 cggcggtgCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTaaa  <  1:468825/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAA  >  NC_000913/3220387‑3220654

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