Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,071,551	G→C	7.6%	P213A (CCA→GCA)	rarA ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,071,551	0	G	C	7.6%	229.0 / 4.6	77	P213A (CCA→GCA)	rarA	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base G (32/39);  new base C (3/3);  total (35/42)
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

GCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTAAA  >  W3110S.gb/1071483‑1071620                                                                     |                                                                      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:4228117/70‑1 (MQ=255)              aCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                                         >  1:4662611/1‑70 (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:3133073/1‑71 (MQ=255)                  ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                    <  1:2480433/71‑1 (MQ=255)                  ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                    >  1:7116929/1‑71 (MQ=255)                  ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                    <  1:6472180/71‑1 (MQ=255)                  ttttCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTg                                                    >  1:7428547/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)                         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:526625/1‑42 (MQ=255)                            aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                       >  1:5506970/1‑42 (MQ=255)                            aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                       >  1:955631/1‑42 (MQ=255)                            aTCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTgg                                                                       >  1:2627328/1‑42 (MQ=255)                             tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                         <  1:5775104/71‑1 (MQ=255)                             tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                         <  1:4561761/71‑1 (MQ=255)                             tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                         <  1:6346084/71‑1 (MQ=255)                             tCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                                         <  1:6476011/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:7159678/1‑38 (MQ=255)                                                               cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                        >  1:6911902/1‑38 (MQ=255)                                                               cTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAtt                                        >  1:3897246/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)                                                                    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)                                                                     |                                                                      GCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGCGGATGCGATCCGCATGGTGACTAAA  >  W3110S.gb/1071483‑1071620

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