Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	196,437	A→G	42.3%	R388R (CGT→CGC)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	196,437	0	A	G	42.3%	28.9 / 74.4	97	R388R (CGT→CGC)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (21/35);  new base G (22/19);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTT  >  minE/196371‑196498                                                                   |                                                              ttCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACg                                                             <  1:16929/69‑1 (MQ=255)               tGACCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACtt                                               >  1:457008/1‑69 (MQ=255)                       agcCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGtc                                       <  1:3258542/69‑2 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:995653/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:1962766/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:2545680/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:2970881/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:1578947/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:1427781/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:1358060/69‑1 (MQ=255)                         ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                     <  1:723463/69‑1 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:1497084/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:267084/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:2996136/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    <  1:1610173/50‑1 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    <  1:521050/50‑1 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:3403040/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:3059702/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    <  1:1124442/50‑1 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:492974/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:2386531/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:1333194/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:394532/1‑50 (MQ=255)                             cccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGc                                                    >  1:1838617/1‑50 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:3029387/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:840895/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:2821701/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:704513/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:323423/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:3284608/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:3332740/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:1634278/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:1219812/68‑1 (MQ=255)                                    cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:1357938/68‑1 (MQ=255)                                     aTCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                           <  1:720098/67‑1 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCggag                                         >  1:1519186/1‑48 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATccc                       >  1:2809779/1‑68 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATccc                       >  1:2599746/1‑68 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:1284012/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:1272162/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:3228723/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:679728/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:3435222/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:1214967/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:152378/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:3408442/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:965421/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:992969/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:2182912/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:2334160/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:814584/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:2658913/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:1818568/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:173710/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:2536775/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:2492017/1‑69 (MQ=255)                                        cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATAAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                      >  1:2655389/1‑69 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:1124207/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:859910/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:499460/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:1179920/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:1162671/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:610998/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:397888/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:1170213/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:632957/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:3434536/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:3376567/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:713857/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:3131164/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:964855/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:2921566/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:732566/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:2710709/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:1693512/69‑1 (MQ=255)                                           gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                   <  1:215814/69‑1 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTcc                                             >  1:450035/1‑39 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   >  1:1095272/1‑65 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   >  1:2000880/1‑65 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   >  1:959093/1‑65 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   >  1:417483/1‑65 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   >  1:347297/1‑65 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   >  1:2773058/1‑65 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATg               >  1:522291/1‑69 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATg               >  1:3121327/1‑69 (MQ=255)                                               ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATg               >  1:1039401/1‑69 (MQ=255)                                                            ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt  <  1:481848/69‑1 (MQ=255)                                                            ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt  <  1:386499/69‑1 (MQ=255)                                                            ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt  <  1:3300884/69‑1 (MQ=255)                                                            ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt  <  1:2364150/69‑1 (MQ=255)                                                              acGTTACGTATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   <  1:1016259/50‑1 (MQ=255)                                                              acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCGCAGtt                   <  1:2020108/50‑1 (MQ=255)                                                              acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   <  1:1027269/50‑1 (MQ=255)                                                              acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   <  1:1692720/50‑1 (MQ=255)                                                              acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   <  1:2967524/50‑1 (MQ=255)                                                               cGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                   <  1:532426/49‑1 (MQ=255)                                                                 ttACGGATATCTGCTACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt  <  1:1015335/64‑1 (MQ=255)                                                                   |                                                              TTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTT  >  minE/196371‑196498

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: ‑