Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,071,554	C→A	7.6%	V212L (GTG→TTG)	rarA ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,071,554	0	C	A	7.6%	255.8 / 4.3	78	V212L (GTG→TTG)	rarA	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base A (3/3);  total (37/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 = 1.00e+00

CATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTAAAGT  >  W3110S.gb/1071485‑1071622                                                                      |                                                                     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:4669931/70‑1 (MQ=255)            aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                           >  1:4662611/1‑70 (MQ=255)            aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                           <  1:4228117/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:3133073/1‑71 (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:2480433/71‑1 (MQ=255)                 tttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCAcc                                                                     <  1:7014341/55‑1 (MQ=255)                 tttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGATGTTGGCAcc                                                                     <  1:655100/55‑1 (MQ=255)                  ttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                      >  1:5382995/1‑69 (MQ=255)                       gCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                      >  1:3558427/1‑64 (MQ=255)                        cTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgat                                              >  1:1703802/1‑71 (MQ=255)                         tATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgat                                              >  1:6233282/1‑70 (MQ=255)                           tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                           <  1:6476011/71‑1 (MQ=255)                           tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                           <  1:4561761/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:3164912/48‑1 (MQ=255)                                  gggCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                           <  1:2517842/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:7047162/1‑63 (MQ=255)                                                    tGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                          >  1:1178518/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:4918991/1‑70 (MQ=255)                                                            gCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGt           >  1:793402/1‑70 (MQ=255)                                                            gCTGGTGCCAACAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGt           >  1:538363/1‑70 (MQ=255)                                                             cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                          >  1:7159678/1‑38 (MQ=255)                                                             cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                          >  1:3897246/1‑38 (MQ=255)                                                             cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                          >  1:331178/1‑38 (MQ=255)                                                             cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                          >  1:6911902/1‑38 (MQ=255)                                                             cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                          >  1:4473195/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)                                                                  ggCACCAGCAGATCATCACTGGCGCAGATGATTTGCa                                      >  1:4945120/1‑37 (MQ=255)                                                                  ggCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATg             >  1:2129786/1‑62 (MQ=255)                                                                  ggCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTa      >  1:4678980/1‑69 (MQ=255)                                                                   gCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTaaa    <  1:4554114/70‑1 (MQ=255)                                                                    cACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGa                          <  1:3963339/47‑1 (MQ=255)                                                                    cACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTAAAGt  <  1:7638466/71‑1 (MQ=255)                                                                      |                                                                     CATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTAAAGT  >  W3110S.gb/1071485‑1071622

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