Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,759,336	C→T	42.7%	T45I (ACA→ATA)	yfjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,759,336	0	C	T	42.7%	37.2 / 91.1	75	T45I (ACA→ATA)	yfjJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/6);  new base T (18/14);  total (55/20)
Fisher's exact test for biased strand distribution p-value = 7.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAG  >  W3110S.gb/2759272‑2759402                                                                 |                                                                   gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAg                                                                  >  1:1306292/1‑67 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:1525920/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:518262/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:1704968/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:451976/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:739799/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:338685/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:3323773/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:2741089/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                              >  1:2234735/1‑71 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaa                                                               >  1:1438410/1‑70 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaa                                                               >  1:670278/1‑70 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaa                                                               >  1:1631425/1‑70 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                                >  1:204036/1‑69 (MQ=255) gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATAGATATAGaaa                                                               >  1:606121/1‑70 (MQ=255)  aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGCTACAGaaaa                                                              <  1:340176/70‑1 (MQ=255)  aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaaa                                                              <  1:1127443/70‑1 (MQ=255)                      aaCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTc                                          <  1:1600772/70‑1 (MQ=255)                      aaCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTc                                          <  1:2850647/70‑1 (MQ=255)                       aCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTc                                          <  1:216950/69‑1 (MQ=255)                          cAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                           >  1:504541/1‑49 (MQ=255)                          cAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                           >  1:1532636/1‑49 (MQ=255)                          cAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                           >  1:764484/1‑49 (MQ=255)                          cAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                           >  1:1321093/1‑49 (MQ=255)                          cAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                           >  1:1252759/1‑49 (MQ=255)                          cAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                           >  1:2200984/1‑49 (MQ=255)                           aaTGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAcc                                      <  1:96288/69‑1 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTct                                         <  1:2423777/64‑2 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:1680633/64‑1 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:1733339/64‑1 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:2256324/64‑1 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:1891764/64‑1 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:1941414/64‑1 (MQ=255)                             tGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:2521494/64‑1 (MQ=255)                                ttCTACGTGTTGATCTGCATGACCCAATTGATATAGa                                                                 <  1:3031819/37‑1 (MQ=255)                                ttCTACGTGTTGATCTGCATGACCCAATTGATATAGa                                                                 <  1:998564/37‑1 (MQ=255)                                ttCTACGTGTTGATCTGCATGACCCAATTGATATAGa                                                                 <  1:728182/37‑1 (MQ=255)                                ttCTACGTGTTGATCTGCATGACCCAATTGATATAGa                                                                 <  1:1733194/37‑1 (MQ=255)                                ttCTACGTGTTGATCTGCATGACCCAATTGATATAGa                                                                 <  1:453340/37‑1 (MQ=255)                                ttCTACGTGTTGATCTGCATGACCCAATTGATATAGa                                                                 <  1:2373676/37‑1 (MQ=255)                                  cTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGa                             >  1:2504686/1‑71 (MQ=255)                                  cTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGa                             >  1:2603289/1‑71 (MQ=255)                                  cTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGa                             >  1:356979/1‑71 (MQ=255)                                  cTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGa                             >  1:2766162/1‑71 (MQ=255)                                  cTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGa                             >  1:283477/1‑71 (MQ=255)                                        gtTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTgg                       >  1:2078388/1‑71 (MQ=255)                                        gtTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTgg                       >  1:1877586/1‑71 (MQ=255)                                        gtTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTgg                       >  1:2789237/1‑71 (MQ=255)                                        gtTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTgg                       >  1:895543/1‑71 (MQ=255)                                           gATCTGCATGTCCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                         <  1:839518/50‑1 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:1957104/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:2154270/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:2825053/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:2044445/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:2675975/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:368886/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:2189354/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:2228758/1‑46 (MQ=255)                                                 cATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                                       >  1:1595940/1‑46 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAgg                                  >  1:1369229/1‑40 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:2088786/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:1573236/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:1191266/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:2193312/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:523400/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:2088276/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:493607/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:42913/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:1094809/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:2573234/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCta                  >  1:246781/1‑56 (MQ=255)                                                            tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGGCTCTGGTGCt                   >  1:3264079/1‑55 (MQ=255)                                                             gATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAg  >  1:1726328/1‑71 (MQ=255)                                                             gATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAg  >  1:3033476/1‑71 (MQ=255)                                                             gATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAg  >  1:2891457/1‑71 (MQ=255)                                                                 |                                                                   GACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAG  >  W3110S.gb/2759272‑2759402

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