Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,071,548	T→G	7.5%	T214P (ACA→CCA)	rarA ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,071,548	0	T	G	7.5%	198.6 / 9.2	80	T214P (ACA→CCA)	rarA	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base G (3/3);  total (39/41)
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.92e-01

TCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTA  >  W3110S.gb/1071480‑1071618                                                                     |                                                                       tCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGtt                                                                       <  1:4518414/70‑1 (MQ=255) tCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGtt                                                                       <  1:2524563/70‑1 (MQ=255) tCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGtt                                                                       <  1:321611/70‑1 (MQ=255)    gCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCa                                                                   <  1:2855773/71‑1 (MQ=255)      cATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGC‑‑CATGAAGTTCACTGGAACATGCTGTTGGCAc                                                                  >  1:2519515/1‑68 (MQ=255)        tAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACcag                                                               <  1:5807988/71‑1 (MQ=255)              aTCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatc                                                         <  1:3155551/71‑1 (MQ=255)              aTCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGat                                                          >  1:3556516/1‑70 (MQ=255)               tCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACcagcag                                                            <  1:3919389/67‑1 (MQ=255)               tCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACcagcag                                                            <  1:2335912/67‑1 (MQ=255)               tCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatca                                                        >  1:3843576/1‑71 (MQ=255)                 aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       <  1:1255816/70‑1 (MQ=255)                 aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       >  1:4662611/1‑70 (MQ=255)                 aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       <  1:4228117/70‑1 (MQ=255)                 aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       <  1:4669931/70‑1 (MQ=255)                 aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       <  1:3314597/70‑1 (MQ=255)                  ccATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGAtcatca                                                     <  1:3106383/71‑1 (MQ=255)                   cATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcatc                                                      <  1:6328326/69‑1 (MQ=255)                   cATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCAc                                                    <  1:811036/71‑1 (MQ=255)                    aTTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcatc                                                      <  1:40515/68‑1 (MQ=255)                    aTTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACt                                                   >  1:6032393/1‑71 (MQ=255)                     ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  <  1:2480433/71‑1 (MQ=255)                     ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  >  1:7428547/1‑71 (MQ=255)                     ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  >  1:7116929/1‑71 (MQ=255)                     ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  <  1:6472180/71‑1 (MQ=255)                     ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  >  1:3133073/1‑71 (MQ=255)                      tttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCAcc                                                                 <  1:7014341/55‑1 (MQ=255)                      tttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGATGTTGGCAcc                                                                 <  1:655100/55‑1 (MQ=255)                       ttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  >  1:5382995/1‑69 (MQ=255)                        tCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGt                                                                        <  1:7223018/46‑1 (MQ=255)                            gCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                  >  1:3558427/1‑64 (MQ=255)                             cTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgat                                          >  1:1703802/1‑71 (MQ=255)                             cTATCGGGCAGGGCGGCATGAAGTTCAATGGAACATGATGTTGGc                                                                    <  1:6922522/45‑1 (MQ=255)                              tATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgat                                          >  1:6233282/1‑70 (MQ=255)                               aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                     >  1:5506970/1‑42 (MQ=255)                               aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                     >  1:526625/1‑42 (MQ=255)                               aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                     >  1:955631/1‑42 (MQ=255)                               aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                     >  1:2627328/1‑42 (MQ=255)                                tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                       <  1:4561761/71‑1 (MQ=255)                                tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                       <  1:6476011/71‑1 (MQ=255)                                tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                       <  1:6346084/71‑1 (MQ=255)                                tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                       <  1:5775104/71‑1 (MQ=255)                                  gggCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcatc                                                      <  1:2629350/54‑1 (MQ=255)                                   ggCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGAtcatca                                                     >  1:799366/1‑54 (MQ=255)                                   ggCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGAtcatca                                                     >  1:7774022/1‑54 (MQ=255)                                    gCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTg                                    >  1:6679099/1‑70 (MQ=255)                                      aGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCAc                                 <  1:2753458/71‑1 (MQ=255)                                       gggCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       <  1:2517842/48‑1 (MQ=255)                                       gggCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                       <  1:3164912/48‑1 (MQ=255)                                          cggcATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCggg                             >  1:5877174/1‑71 (MQ=255)                                             cATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatca                                                        <  1:7233450/41‑1 (MQ=255)                                                gAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGc                            <  1:3929983/66‑1 (MQ=255)                                                gAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCg                        >  1:6032556/1‑70 (MQ=255)                                                 aaGTTCACTGGAACATGCTGTTGGCACCAGCAGatcatc                                                      >  1:795982/1‑39 (MQ=255)                                                 aaGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCgg                       <  1:7350356/70‑1 (MQ=255)                                                 aaGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                      <  1:5511826/71‑1 (MQ=255)                                                   gTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAttt                                     >  1:2999039/1‑54 (MQ=255)                                                    ttCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCg                        <  1:1977253/66‑1 (MQ=255)                                                    ttCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGc                   <  1:524046/71‑1 (MQ=255)                                                      cACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGAt                     >  1:586037/1‑67 (MQ=255)                                                      cACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCg                  <  1:6567320/70‑1 (MQ=255)                                                      cACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCg                  <  1:3614096/70‑1 (MQ=255)                                                      cACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGa                 >  1:4181197/1‑71 (MQ=255)                                                        cTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATc               <  1:5881263/71‑1 (MQ=255)                                                         tGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                      >  1:1178518/1‑63 (MQ=255)                                                         tGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                      >  1:7047162/1‑63 (MQ=255)                                                          ggAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCg             >  1:1333190/1‑71 (MQ=255)                                                             aCATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCAt          >  1:5014439/1‑71 (MQ=255)                                                              tctGCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCAc                                 <  1:3913678/45‑1 (MQ=25)                                                              tctGCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCAc                                 <  1:729546/45‑1 (MQ=25)                                                              cATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATg         >  1:2170777/1‑71 (MQ=255)                                                                tGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGt       <  1:1866169/71‑1 (MQ=255)                                                                 gCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                      <  1:7928317/55‑1 (MQ=255)                                                                 gCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCa           <  1:7588223/66‑1 (MQ=255)                                                                 gCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGt       >  1:538363/1‑70 (MQ=255)                                                                 gCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGt       >  1:793402/1‑70 (MQ=255)                                                                 gCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGt       >  1:4918991/1‑70 (MQ=255)                                                                  cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                      >  1:4473195/1‑38 (MQ=255)                                                                  cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                      >  1:6911902/1‑38 (MQ=255)                                                                  cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                      >  1:3897246/1‑38 (MQ=255)                                                                  cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                      >  1:7159678/1‑38 (MQ=255)                                                                  cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                      >  1:331178/1‑38 (MQ=255)                                                                   tGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCGGCGGATGCGAt                >  1:6395971/1‑59 (MQ=255)                                                                   tGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                      <  1:6609312/53‑1 (MQ=255)                                                                     ttGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTa  <  1:4538269/71‑1 (MQ=255)                                                                     |                                                                       TCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTA  >  W3110S.gb/1071480‑1071618

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