Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,236,524	G→A	D105N (GAT→AAT)	sanA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,236,524	0	G	A	98.9%	257.4 / ‑4.3	88	D105N (GAT→AAT)	sanA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (58/29);  total (59/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TAAATTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAGATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTG  >  W3110S.gb/2236456‑2236589                                                                     |                                                                  taaATTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAAt                                                                 >  1:2164576/1‑71 (MQ=255) taaATTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAAt                                                                 >  1:2500734/1‑71 (MQ=255) taaATTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAAt                                                                 >  1:1467684/1‑71 (MQ=255)   aaTTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAAtt                                                                >  1:2486035/1‑70 (MQ=255)    aTTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTa                                                              <  1:366876/71‑1 (MQ=255)    aTTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTa                                                              <  1:4674734/71‑1 (MQ=255)    aTTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTa                                                              <  1:3997191/71‑1 (MQ=255)      tATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAAt                                                            >  1:2746609/1‑71 (MQ=255)        tCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATc                                                           <  1:3000620/70‑1 (MQ=255)            ttACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACGATGCGCAAAAATTTAATCgctgc                                                      >  1:4421237/1‑71 (MQ=255)              aCTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCgctgct                                                     <  1:4649169/70‑1 (MQ=255)               cTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCgctgct                                                     <  1:3920150/69‑1 (MQ=255)                tGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAAttt                                                               <  1:3408887/58‑1 (MQ=255)                tGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAAttt                                                               <  1:4678644/58‑1 (MQ=255)                     ggcgATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAAt                                                            >  1:1996806/1‑56 (MQ=255)                       cgATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGAcc                                           <  1:2086713/71‑1 (MQ=255)                         aTAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCa                                         <  1:1221620/71‑1 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGCCCCAt                                        >  1:592612/1‑70 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCAt                                        >  1:3033732/1‑70 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCAt                                        >  1:2764957/1‑70 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCAt                                        >  1:2880656/1‑70 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCAt                                        >  1:2082380/1‑70 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCAt                                        >  1:4042110/1‑70 (MQ=255)                           aaCGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCAt                                        >  1:1336476/1‑70 (MQ=255)                                 ttGCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCgctgct                                                     <  1:755774/51‑1 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGGTGGTGTCGACCCATCAGATATTg                               >  1:3799895/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTgg                                                   >  1:2786408/1‑51 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGa                                             >  1:2410548/1‑57 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGAcc                                           >  1:4175926/1‑59 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:4087393/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:3339655/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:3390392/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:3400539/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:3537808/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:3592835/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:435621/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:2493215/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:1361417/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:159296/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:1695766/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:579857/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:2004596/1‑71 (MQ=255)                                   gCAAAGTTATAATGAGCCGATGACCATGAGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTg                               >  1:1700444/1‑71 (MQ=255)                                     aaaGTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGtt                             <  1:4510183/71‑1 (MQ=255)                                        gTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTtctc                          <  1:3321493/71‑1 (MQ=255)                                        gTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTtctc                          <  1:4514202/71‑1 (MQ=255)                                        gTTATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTtctc                          <  1:3419759/71‑1 (MQ=255)                                         ttATAATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCg                         <  1:702872/71‑1 (MQ=255)                                          tataATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGa                        >  1:2072646/1‑71 (MQ=255)                                          tataATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGa                        >  1:4189512/1‑71 (MQ=255)                                          tataATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGATGGTGTCGACCCATCAGATATTGTTCTCGa                        >  1:1476426/1‑71 (MQ=255)                                           ataATGAGCCGATGACCATGCGCAAAGATTTAATCgc                                                         >  1:4265473/1‑37 (MQ=255)                                           ataATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGAt                       >  1:4501923/1‑71 (MQ=255)                                           ataATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGAt                       >  1:837726/1‑71 (MQ=255)                                           ataATGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGAt                       >  1:1110480/1‑71 (MQ=255)                                             aaTGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCATATATTGTTCTCGATTa                     >  1:742551/1‑71 (MQ=255)                                              aTGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTAc                    >  1:875190/1‑71 (MQ=255)                                               tGAGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACg                   >  1:1170310/1‑71 (MQ=255)                                                 aGCTGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACg                   >  1:3720746/1‑69 (MQ=255)                                                 aGCCGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACg                   >  1:2709788/1‑69 (MQ=255)                                                   ccGATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAgg               >  1:3136801/1‑71 (MQ=255)                                                     gATGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCt             <  1:3654450/71‑1 (MQ=255)                                                      aTGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGAt                       <  1:3637906/60‑1 (MQ=255)                                                      aTGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCtt            >  1:2627658/1‑71 (MQ=255)                                                       tGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGAt                       <  1:1927052/59‑1 (MQ=255)                                                       tGACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGc              <  1:2540652/68‑1 (MQ=255)                                                        gACCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTc          >  1:3324289/1‑71 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTTGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         >  1:3735915/1‑71 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTGGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         >  1:3615143/1‑71 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         >  1:3274025/1‑71 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         >  1:1592293/1‑71 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         <  1:3317395/71‑1 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         <  1:1667659/71‑1 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         <  1:4126887/71‑1 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         <  1:2268246/71‑1 (MQ=255)                                                         aCCATGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         <  1:2097366/71‑1 (MQ=255)                                                             tGCGCAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCg         <  1:3362434/67‑1 (MQ=255)                                                               cgcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCt   >  1:1146020/1‑71 (MQ=255)                                                               cgcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCt   >  1:1770578/1‑71 (MQ=255)                                                               cgcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCt   >  1:1014498/1‑71 (MQ=255)                                                               cgcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCt   >  1:577476/1‑71 (MQ=255)                                                               cgcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCt   >  1:2061925/1‑71 (MQ=255)                                                               cgcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCt   >  1:4261366/1‑71 (MQ=255)                                                               agcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGc    >  1:1767120/2‑70 (MQ=255)                                                               agcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGc    >  1:2959582/2‑70 (MQ=255)                                                                gcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTg  <  1:4252802/71‑1 (MQ=255)                                                                gcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTg  <  1:4181550/71‑1 (MQ=255)                                                                gcAAAAATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTg  <  1:4018800/71‑1 (MQ=255)                                                                     |                                                                  TAAATTATCTATTACTGAGCGGCGATAACGCATTGCAAAGTTATAATGAGCCGATGACCATGCGCAAAGATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTG  >  W3110S.gb/2236456‑2236589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑