Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,247,5860TC71.9% 86.9 / 38.6 57L456L (TTA→TTGtreAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base C (0/41);  total (7/50)
Fisher's exact test for biased strand distribution p-value = 4.33e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCG  >  W3110S.gb/1247529‑1247635
                                                         |                                                 
gTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAAc                                                <  1:3058144/61‑1 (MQ=255)
gTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAAc                                                <  1:2364312/61‑1 (MQ=255)
gTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAAc                                                <  1:1515807/61‑1 (MQ=255)
gTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAAc                                                <  1:2481398/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcg                                          <  1:389370/61‑2 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2644477/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2162822/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:388472/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2298620/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2379346/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2427367/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:3079398/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2073319/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2649072/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2737893/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2744517/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2761196/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2791008/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2976053/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2914902/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1622642/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1089550/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:942446/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:884333/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:739635/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1388881/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1403726/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1519034/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:723146/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:602806/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2083266/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:575701/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1704512/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:177353/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:465131/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1821548/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1888062/61‑1 (MQ=255)
      cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:190962/61‑1 (MQ=255)
       gCCACCTCTTTTTGCCCGTAGTTTTTTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1183846/60‑1 (MQ=255)
       gCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1196980/60‑1 (MQ=255)
                ttttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:3068060/51‑1 (MQ=255)
                 tttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2267652/50‑1 (MQ=255)
                    gCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1580790/47‑1 (MQ=255)
                        gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATccc                       >  1:2793143/1‑62 (MQ=255)
                        gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATccc                       >  1:132532/1‑62 (MQ=255)
                        gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATccc                       >  1:2451798/1‑62 (MQ=255)
                        gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATccc                       >  1:180221/1‑62 (MQ=255)
                        gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATccc                       >  1:1651875/1‑62 (MQ=255)
                        gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATccc                       >  1:153609/1‑62 (MQ=255)
                         tAGTTTTGTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:1270252/42‑1 (MQ=255)
                                gTAATCCTTCTGTCGCGACCCACTGCAACGGTGcc                                          <  1:2964326/35‑1 (MQ=255)
                                             cgcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCg  >  1:369710/1‑62 (MQ=255)
                                              gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc    <  1:2135644/59‑1 (MQ=255)
                                              gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc    <  1:635576/59‑1 (MQ=255)
                                              gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc    <  1:3063395/59‑1 (MQ=255)
                                              gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc    <  1:3018473/59‑1 (MQ=255)
                                              gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc    <  1:2830888/59‑1 (MQ=255)
                                                         |                                                 
GTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCG  >  W3110S.gb/1247529‑1247635

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: