Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,356,628	C→T	A310T (GCC→ACC)	ycjJ ←	putrescine importer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,356,628	0	C	T	94.4%	125.1 / ‑4.2	36	A310T (GCC→ACC)	ycjJ	putrescine importer
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (12/22);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGCATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGAACAGCTTACC  >  REL606/1356555‑1356700                                                                          |                                                                         tGCACATAGCCAAACACGCGCCCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTAt                                                                        >  1:949813/1‑76 (MQ=255)    acaTAGCCAAACACGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTAtgtg                                                                     <  2:494163/76‑1 (MQ=255)    acaTAGCCAAACACGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTAtgtg                                                                     >  2:111964/1‑76 (MQ=255)       tAGCCAAACACGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGAGACACGCTGGCATGTGagg                                                                  >  2:1347335/1‑76 (MQ=255)           cAAACACGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGaggccag                                                              <  2:1336745/76‑1 (MQ=255)              acacGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGaggccaggcc                                                           <  2:969913/76‑1 (MQ=255)              acacGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGaggccaggcc                                                           <  2:396843/76‑1 (MQ=255)              acacGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGaggccaggcc                                                           >  2:826920/1‑76 (MQ=255)                acGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGgccaggccag                                                         <  2:706049/76‑1 (MQ=255)                   cgcTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACg                                                      >  1:1392298/1‑76 (MQ=255)                    gcTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGc                                                     >  1:587269/1‑76 (MQ=255)                     cTCCGAAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCt                                                    <  2:635358/76‑1 (MQ=255)                              cacaTTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTc                                           >  1:873932/1‑76 (MQ=255)                              cacaTGGTCGCGCCCCATCACATACAGCAGACGTGCCACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTc                                           <  2:992781/76‑1 (MQ=255)                                caTTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTcac                                         <  2:702639/76‑1 (MQ=255)                                 aTTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTcaca                                        <  1:1054705/76‑1 (MQ=255)                                   tGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACaaa                                      >  1:1489559/1‑76 (MQ=255)                                   tGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACaaa                                      <  2:1249203/76‑1 (MQ=255)                                            ccATCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGcn                             <  1:184104/76‑2 (MQ=255)                                              aTCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCaaa                           >  2:321537/1‑76 (MQ=255)                                              aTCACATACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCaaa                           <  1:1347335/76‑1 (MQ=255)                                                  caTACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGaa                       >  2:1280364/1‑76 (MQ=255)                                                   aTACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGaaa                      <  1:185979/76‑1 (MQ=255)                                                   aTACAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGaaa                      <  2:756238/76‑1 (MQ=255)                                                     aCAGCGGGCGTGACCCGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAAt                    <  1:272491/76‑1 (MQ=255)                                                     aCAGCAGACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAAt                    <  2:839179/76‑1 (MQ=255)                                                       agcagACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCg                  <  1:1322855/76‑1 (MQ=255)                                                        gcagACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGa                 <  1:691518/76‑1 (MQ=255)                                                         cagACGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGAc                <  2:1267169/76‑1 (MQ=255)                                                           gACGTGACACGCTGGCATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTg              <  1:1362036/76‑1 (MQ=255)                                                            aCGTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTgg             >  2:1145781/1‑76 (MQ=255)                                                              gTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGaa           <  2:781245/76‑1 (MQ=255)                                                              gTGACACGCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGaa           <  1:1318002/76‑1 (MQ=255)                                                                     gCTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGAACAGCTTa    <  2:315715/76‑1 (MQ=255)                                                                      cTGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGAACAGCTTAc   >  1:1347102/1‑76 (MQ=255)                                                                       tGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGAACAGCTTAcc  >  1:637227/1‑76 (MQ=255)                                                                       tGGTATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGAACAGCTTAcc  >  1:131553/1‑76 (MQ=255)                                                                          |                                                                         TGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCGCGCCCCATCACATACAGCAGACGTGACACGCTGGCATGTGAGGCCAGGCCAGACGCTAACGTGTTCACAAACGTGGTGCAAAGGAAAATCGACTGGAACAGCTTACC  >  REL606/1356555‑1356700

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