Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,569,965	T→C	46.8%	Q191Q (CAA→CAG)	yddW ←	predicted liprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,569,965	0	T	C	46.8%	4.9 / 84.1	62	Q191Q (CAA→CAG)	yddW	predicted liprotein
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base C (12/17);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CCTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACC  >  W3110S.gb/1569898‑1570024                                                                    |                                                            ccTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTg                                                           >  1:3310291/1‑70 (MQ=255) ccTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTg                                                           >  1:3188300/1‑70 (MQ=255) ccTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTg                                                           >  1:2698790/1‑70 (MQ=255) ccTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTg                                                           >  1:1293144/1‑70 (MQ=255) ccTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTg                                                           >  1:2149644/1‑70 (MQ=255) ccTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTg                                                           >  1:205826/1‑70 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTgga                                                         >  1:1624740/1‑69 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:897892/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:896690/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:3250238/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:221002/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:2124337/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:2075965/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:1888820/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                         >  1:1157880/1‑71 (MQ=255)  cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca                                                         <  1:474253/71‑1 (MQ=255)                       gTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt                                    <  1:764955/71‑1 (MQ=255)                       gTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt                                    <  1:2208786/71‑1 (MQ=255)                       gTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt                                    <  1:1513423/71‑1 (MQ=255)                       gTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt                                    <  1:397890/71‑1 (MQ=255)                        tCGAGGACAAAGCGATCGCCAGACCTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt                                    <  1:827707/70‑1 (MQ=255)                                aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGCCGCTCGCCGGTTGTTGAGAc                               >  1:3342177/1‑67 (MQ=255)                                aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACa                              >  1:1814693/1‑68 (MQ=255)                                aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACa                              >  1:1495978/1‑68 (MQ=255)                                aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACa                              >  1:498367/1‑68 (MQ=255)                                aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACa                              >  1:742146/1‑68 (MQ=255)                                aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACa                              >  1:1052275/1‑68 (MQ=255)                                 aaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACa                              <  1:2951650/67‑1 (MQ=255)                                  aGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTg                         <  1:355861/71‑1 (MQ=255)                                  aGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTg                         <  1:351901/71‑1 (MQ=255)                                  aGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTg                         <  1:1072417/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTg                         <  1:1007111/70‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:1133835/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:1997841/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:790233/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:1535538/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:679515/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:66058/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:519378/71‑1 (MQ=255)                                   gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:3313133/71‑1 (MQ=255)                                                   cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACTGAGTGCTATTCAGTTCCCTGa         <  1:603397/70‑1 (MQ=255)                                                   cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa         <  1:3326696/70‑1 (MQ=255)                                                   cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa         <  1:2466905/70‑1 (MQ=255)                                                   cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa         <  1:253629/70‑1 (MQ=255)                                                   cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa         <  1:2066282/70‑1 (MQ=255)                                                   cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa         <  1:779124/70‑1 (MQ=255)                                                    tGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa         <  1:2749562/69‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTTTTCAGTTCCCTGATAGt     <  1:1120845/71‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt     <  1:2406777/71‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt     <  1:1347273/71‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt     <  1:2098089/71‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt     <  1:2080557/71‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt     <  1:522752/71‑1 (MQ=255)                                                      aTCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt     <  1:2972810/71‑1 (MQ=255)                                                         cAGTCGCGGTGTTTCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTAc   >  1:108359/1‑70 (MQ=255)                                                         cAGTCGCGGTGTTGCACATAGACGCTCGCCGGttgttg                                   >  1:844435/1‑38 (MQ=255)                                                         cAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTAcc  >  1:3164305/1‑71 (MQ=255)                                                         cAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTAcc  >  1:746664/1‑71 (MQ=255)                                                         cAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTAcc  >  1:67040/1‑71 (MQ=255)                                                         cAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTAcc  >  1:1924878/1‑71 (MQ=255)                                                         cAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTAc   >  1:1586783/1‑70 (MQ=255)                                                           ttCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc                        <  1:2367360/46‑1 (MQ=255)                                                                    |                                                            CCTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACC  >  W3110S.gb/1569898‑1570024

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