Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,858,954	T→C	6.9%	T83T (ACA→ACG)	yibK ←	predicted rRNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,858,954	0	T	C	6.9%	151.7 / 2.6	58	T83T (ACA→ACG)	yibK	predicted rRNA methylase
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base C (2/2);  total (23/35)
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.38e-01

GCCGCGTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGATAGTCATGAT  >  W3110S.gb/3858888‑3859024                                                                   |                                                                       gCCGCGTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTa                                                                     >  1:2426487/1‑70 (MQ=255) gCCGCGTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTAc                                                                    <  1:3293927/71‑1 (MQ=255) gCCGCGTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTAc                                                                    <  1:3111550/71‑1 (MQ=255) gCCGCGTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTAc                                                                    <  1:3036381/71‑1 (MQ=255)  ccGCGTGTTTCC‑GGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCt                                                                  <  1:1511936/71‑1 (MQ=255)   cgcgTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGgtg                                                                       >  1:238451/1‑66 (MQ=255)   cgcgTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCt                                                                  <  1:474212/71‑1 (MQ=255)   cgcgTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCt                                                                  <  1:2807936/71‑1 (MQ=255)   cgcgTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCt                                                                  >  1:1371759/1‑71 (MQ=255)       tgtTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCttt                                                                <  1:3383381/69‑1 (MQ=255)       tgtTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCg                                                              <  1:902812/71‑1 (MQ=255)          ttCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGCGTACCTTTCgtgg                                                           >  1:3450290/1‑71 (MQ=255)          ttCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGCGTACCTTTCgtgg                                                           >  1:462179/1‑71 (MQ=255)            ccGGTCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          >  1:392153/1‑70 (MQ=255)            ccGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:1498279/70‑1 (MQ=255)            ccGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:1385289/70‑1 (MQ=255)            ccGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:1949444/70‑1 (MQ=255)             cGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggtg                                                         <  1:2724909/70‑1 (MQ=255)              gggCCGAACATCAGATAATCGCCATCCTTATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAg                                                       <  1:1254589/71‑1 (MQ=255)              gggCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAg                                                       <  1:1589928/71‑1 (MQ=255)               ggCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTa                                                                     >  1:1946674/1‑56 (MQ=255)               ggCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTa                                                                     >  1:3741105/1‑56 (MQ=255)                gCCGAACATCAGATAATCGCCATCCTGATAGATTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAggg                                                     >  1:3578098/1‑71 (MQ=255)                      cATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:1051494/60‑1 (MQ=255)                      cATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACa                                               <  1:131014/71‑1 (MQ=255)                        tCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGa                                                        >  1:1376941/1‑60 (MQ=255)                        tCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGaa                                                 <  1:1306908/67‑1 (MQ=255)                          aGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGTTGAGGGCGAACAGgcg                                           >  1:2462394/1‑71 (MQ=255)                          aGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGgcg                                           >  1:446160/1‑71 (MQ=255)                          aGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGgcg                                           >  1:2404984/1‑71 (MQ=255)                           gATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGgcgc                                          <  1:1433345/71‑1 (MQ=255)                            aTAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTAc                                                                    <  1:3767873/44‑1 (MQ=255)                            aTAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTAc                                                                    <  1:1321452/44‑1 (MQ=255)                                tCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTgggg                                     <  1:2016334/71‑1 (MQ=255)                                tCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTggg                                      >  1:1021989/1‑70 (MQ=255)                                     aTCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCg                                                              >  1:1755766/1‑41 (MQ=255)                                     aTCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCg                                                              >  1:1760223/1‑41 (MQ=255)                                     aTCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCg                                                              >  1:3470488/1‑41 (MQ=255)                                      tCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:2074489/44‑1 (MQ=255)                                      tCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:2863492/44‑1 (MQ=255)                                        cTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtggt                                                          <  1:259483/42‑1 (MQ=255)                                        cTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgtgg                                                           >  1:2494185/1‑41 (MQ=255)                                        cTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTg                             >  1:3400587/1‑71 (MQ=255)                                        cTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTg                             >  1:1784299/1‑71 (MQ=255)                                          gATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGAtt                                  >  1:3204354/1‑64 (MQ=255)                                          gATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGAtt                                  >  1:579317/1‑64 (MQ=255)                                           aTAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGc                            <  1:107690/69‑1 (MQ=255)                                            tAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgt                                                             <  1:1552219/35‑1 (MQ=255)                                            tAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCgt                                                             <  1:2050514/35‑1 (MQ=255)                                            tAGCTTACGGCGCTGTGAGCAGGCGTACCTTTCGTGGTGAGGGCAAACAGACGCTGGGGATTTTCTGCTTc                         <  1:1197503/71‑1 (MQ=255)                                            tAGCTTACGGCGCTGTGAGCAGGCGTACCTTTCGTGGTGAGGGCAAACAGACGCTGGGGATTTTCTGCTTc                         <  1:2464645/71‑1 (MQ=255)                                              gCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCg                        <  1:3481636/70‑1 (MQ=255)                                              gCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCg                        <  1:2540774/70‑1 (MQ=255)                                                 tACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAgg                     >  1:3510869/1‑70 (MQ=255)                                                  aCGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAgg                     <  1:1532968/69‑1 (MQ=255)                                                  aCGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGa                    <  1:2133370/70‑1 (MQ=255)                                                           tgAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGa            <  1:3488719/69‑1 (MQ=255)                                                           tgAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGa            >  1:1277725/1‑69 (MQ=255)                                                                 ggTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGATAGTCatg    <  1:2840669/71‑1 (MQ=255)                                                                  gtgtACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGATAGTCatga   >  1:283777/1‑71 (MQ=255)                                                                   tgtACCTTTCGTGGTGAGGGCTAACAGGCGCTGGGGAt                                   <  1:1658847/38‑1 (MQ=255)                                                                   tgtACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGATAGTCatgat  <  1:1678217/71‑1 (MQ=255)                                                                   |                                                                       GCCGCGTGTTTCCGGGCCGAACATCAGATAATCGCCATCCTGATAGCTTACGGCGCTGTGAGCAGGTGTACCTTTCGTGGTGAGGGCGAACAGGCGCTGGGGATTTTCTGCTTCGAGGAACGCGCGATAGTCATGAT  >  W3110S.gb/3858888‑3859024

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