Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,480,143	T→C	5.5%	Y9H (TAT→CAT)	yjgL →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,480,143	0	T	C	5.5%	260.1 / 5.8	90	Y9H (TAT→CAT)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/42);  new base C (2/3);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

AAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGT  >  W3110S.gb/4480074‑4480208                                                                      |                                                                  aaaTTTGTAATGCTTAGACGTTCTTATTCCCTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAatta                                                                  >  1:4386536/1‑71 (MQ=255) aaaTTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAatta                                                                  >  1:249495/1‑71 (MQ=255) aaaTTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAatt                                                                   <  1:1063441/70‑1 (MQ=255)    tttGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtc                                                               >  1:4014620/1‑71 (MQ=255)    tttGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtc                                                               >  1:4200190/1‑71 (MQ=255)      tGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATCATtct                                                              >  1:3553645/1‑70 (MQ=255)               tAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAAcaca                                                       >  1:1010884/1‑68 (MQ=255)               tAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAAcaca                                                       >  1:106397/1‑68 (MQ=255)               tAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacatta                                                    >  1:1030830/1‑71 (MQ=255)               tAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacatta                                                    >  1:2071795/1‑71 (MQ=255)                 gACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  <  1:4537510/71‑1 (MQ=255)                 gACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  <  1:3471175/71‑1 (MQ=255)                 gACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  <  1:3816562/71‑1 (MQ=255)                 gACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  <  1:850348/71‑1 (MQ=255)                 gACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattac                                                   <  1:729413/70‑1 (MQ=255)                 gACGTTCTTATTCACTCAAGGAGATCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  <  1:3346402/71‑1 (MQ=255)                  aCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:4590191/59‑1 (MQ=255)                  aCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATCATTCTCAACacatta                                                    <  1:3164310/68‑1 (MQ=255)                  aCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATCATTCTCAACacatta                                                    <  1:1582258/68‑1 (MQ=255)                    gTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATCATTCTCAACacattaca                                                  >  1:652798/1‑68 (MQ=255)                     ttCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:2291345/56‑1 (MQ=255)                     ttCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:3784631/56‑1 (MQ=255)                      tCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacatta                                                    <  1:1038499/64‑1 (MQ=255)                      tCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGc                                             >  1:2522905/1‑71 (MQ=255)                       cTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattacatt                                                >  1:1895994/1‑67 (MQ=255)                       cTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattacat                                                 >  1:2973938/1‑66 (MQ=255)                       cTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGcc                                            >  1:2776742/1‑71 (MQ=255)                          aTTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:4371419/51‑1 (MQ=255)                           ttCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:941117/50‑1 (MQ=255)                           ttCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:1660487/50‑1 (MQ=255)                           ttCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc                                                             <  1:3831240/50‑1 (MQ=255)                           ttCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacat                                                      <  1:4925339/57‑1 (MQ=255)                             cACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAacac                                                        >  1:3133782/1‑53 (MQ=255)                              aCTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacat                                                      <  1:3009414/54‑1 (MQ=255)                               cTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAg                                              <  1:3599738/61‑1 (MQ=255)                               cTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAg                                              <  1:3605793/61‑1 (MQ=255)                                tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAAttt                                    >  1:3808660/1‑70 (MQ=255)                                tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATCATTCTCAac                                                          <  1:3288244/48‑1 (MQ=255)                                  aaGGAGTTCGTCATGAGCAATATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTa                                  >  1:1934988/1‑70 (MQ=255)                                  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCaa                                                           >  1:434982/1‑45 (MQ=255)                                  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCaa                                                           >  1:3379846/1‑45 (MQ=255)                                  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTaa                                 <  1:4591880/71‑1 (MQ=255)                                  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTa                                  >  1:952606/1‑70 (MQ=255)                                  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTa                                  >  1:3911511/1‑70 (MQ=255)                                  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTa                                  >  1:330361/1‑70 (MQ=255)                                    ggAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTaaa                                >  1:161635/1‑70 (MQ=255)                                     gAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACa                              <  1:3828200/71‑1 (MQ=255)                                      aGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattac                                                   <  1:553718/49‑1 (MQ=255)                                      aGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattac                                                   <  1:4843477/49‑1 (MQ=255)                                      aGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattac                                                   <  1:400665/49‑1 (MQ=255)                                        ttCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACa                                        <  1:928068/58‑1 (MQ=255)                                         tCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACaaaa                           >  1:5058311/1‑70 (MQ=255)                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCa                                                            >  1:4017426/1‑36 (MQ=255)                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCg                                           >  1:731375/1‑53 (MQ=255)                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTaaa                                >  1:3647505/1‑64 (MQ=255)                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTaa                                 >  1:919683/1‑63 (MQ=255)                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACaaaa                           >  1:962967/1‑69 (MQ=255)                                          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACaaa                            >  1:4633493/1‑68 (MQ=255)                                            tCATGAGCAAAATATCAGATTTAAATTATTCTCAaca                                                         <  1:400326/37‑1 (MQ=255)                                              aTGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTg                      <  1:2190762/70‑1 (MQ=255)                                               tGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  >  1:4645520/1‑41 (MQ=255)                                               tGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  <  1:2798691/41‑1 (MQ=255)                                               tGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  >  1:2467049/1‑41 (MQ=255)                                               tGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca                                                  >  1:1328086/1‑41 (MQ=255)                                                 aGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGt                  >  1:2564454/1‑71 (MQ=255)                                                 aGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGt                  >  1:583086/1‑71 (MQ=255)                                                   cAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGttt                <  1:3196497/71‑1 (MQ=255)                                                    aaaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGttt                <  1:1975669/70‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGtttt               >  1:4355277/1‑70 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:117640/71‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:3499161/71‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:3781527/71‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:123647/71‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:2698667/71‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:3525472/71‑1 (MQ=255)                                                     aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa              <  1:4919822/71‑1 (MQ=255)                                                      aaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTaa             <  1:888388/71‑1 (MQ=255)                                                      aaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTaa             <  1:1300224/71‑1 (MQ=255)                                                        tatCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTaaa            >  1:4228174/1‑70 (MQ=255)                                                        tatCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAAt           <  1:1301328/71‑1 (MQ=255)                                                        tatCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAAt           <  1:4584229/71‑1 (MQ=255)                                                          tCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATAc         >  1:4956872/1‑71 (MQ=255)                                                            aGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCt       >  1:548346/1‑71 (MQ=255)                                                            aGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCt       >  1:3008388/1‑71 (MQ=255)                                                              aTTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTgg     >  1:4218057/1‑71 (MQ=255)                                                               tttAAATTATTCTCAACACATTACATTAGCCGACAATTTTaaa                                >  1:2516768/1‑43 (MQ=255)                                                               tttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGc    >  1:1931819/1‑71 (MQ=255)                                                                ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACaaaaaa                         >  1:1207445/1‑49 (MQ=255)                                                                ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGc    >  1:4419407/1‑70 (MQ=255)                                                                ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGc    >  1:3198430/1‑70 (MQ=255)                                                                  aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCgt  <  1:3787243/70‑1 (MQ=255)                                                                  aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCgt  <  1:3784817/70‑1 (MQ=255)                                                                   aaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATAc         >  1:2756059/1‑62 (MQ=255)                                                                      tattCTCAACACATTACATTAGCCGACAATTTTAAACa                              <  1:4486129/38‑1 (MQ=255)                                                                      tattCTCAACACATTACATTAGCCGACAATTTTAAACa                              <  1:3116199/38‑1 (MQ=255)                                                                      |                                                                  AAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGT  >  W3110S.gb/4480074‑4480208

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