Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,356,5120GC25.6% 84.8 / 11.5 43L408V (CTG→GTG) glmMPhosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (27/5);  new base C (1/10);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 1.49e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGCAAGCGCCTTCAGACACAAACTTCGCCGACCAGTTTGGCCAGCGCTTCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACTTGCGCAGGG  >  ACB122/5356379‑5356643
                                                                                                                                     |                                                                                                                                   
gcTGGCAAGCGCCTTCAGACACAAACTTCGCCGACCAGTTTGGCCAGCGCTTCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGaac                                                                                                                                  >  3:14544/1‑137 (MQ=255)
  tGGCAAGCGCCTTCAGACACAAACTTCGCCGACCAGTTTGGCCAGCGCTTCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCacac                                                                                                                                >  7:225450/1‑137 (MQ=255)
                 gACACAAACTTCGCCGACCAGAATGGCCAGCGCTTCAGCAAGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGAACCCCACCAACGGCACGGTACCGCACTTGCGCAACACCACACCACCCCCCCCCGCCa                                                                                                                 >  5:378140/1‑137 (MQ=255)
                      aaaCTTCGCCGACCAGTTTGGCCAGCGCTTCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGTCCCCGCCCCGCCAGCGCt                                                                                                            >  4:276568/1‑137 (MQ=255)
                                         ggCCAGCGCTTCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTgg                                                                                         >  8:64509/1‑137 (MQ=255)
                                                  ttCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGa                                                                                >  1:63037/1‑137 (MQ=255)
                                                  ttCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGa                                                                                >  5:11737/1‑137 (MQ=255)
                                                  ttCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGa                                                                                >  3:341590/1‑137 (MQ=255)
                                                   tCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTGCGGTAACCTTGTCACTGGCTGCCTTGaa                                                                               >  3:269810/1‑136 (MQ=255)
                                                    cagcagGGCCCCGCCCCGGGCTTTCCGCCCCGCCCACACCAAGCCCCCCCCCCAACGGCTCGTTCCCGGACTTCCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGAcg                                                                              <  6:305820/137‑1 (MQ=255)
                                                    cAGCATGGGCCCGCACCTGTCTTGCTTCCTCGCACTCAACCATGCCCCCCACCACCGGCTCGGTACCGGAATTGCGCAACACCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGAcg                                                                              <  7:244509/137‑1 (MQ=255)
                                                    cAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATCACCCCCACCACCGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGAcg                                                                              <  8:99188/137‑1 (MQ=255)
                                                      gCATGGCCCCGCACCTGGCTTTCGTCCTCGCCCGCACCCATCACCCGCACCACCGGCTCGGTACCGGACTTGCGCAACAGCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGAcggc                                                                            <  8:225450/137‑1 (MQ=255)
                                                       cATGGCCCCGCACCTGGCTTGCGTCCTCCCCCCCAACCATCACCCGCACCACCGGCGCGGTACCGGACTTGCGCAACCCCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGAcggcg                                                                           <  5:90331/137‑1 (MQ=255)
                                                         ggggCCCGCCCGGGGATTGCGCCCTCCCCCTCAACCATGCCCCGCCCCACCGGCGCTGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCggg                                                                         <  7:189355/136‑1 (MQ=255)
                                                            gCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCACCGGCTCGGTACCGGACTTGCGCACCAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTg                                                                      >  2:213902/1‑137 (MQ=255)
                                                               cGCCCCTGGCTTTCGTCCTCGCCCTCACCCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCTCCAGCACACGACCCCGCCCCGCCAGCGCGTCGGTAACCTTGGCCCTGGCTTCCGTGACGGCGGGGTGGTc                                                                   >  6:306816/1‑137 (MQ=255)
                                                                         tttCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCg                                                         >  1:368459/1‑137 (MQ=255)
                                                                         ttgcGTCCCCGCCCCCACCCACCCCCCCCCCCACCGGCTCGTTCCCGGACTGCCGCAACACCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCg                                                         <  8:327514/134‑1 (MQ=255)
                                                                                              tgccccGAACCAACGGCTCGGTACCGGACTTGCGCAACCGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAgcgc                                    <  7:50543/134‑1 (MQ=255)
                                                                                              accccccCACCAACGGCCCGGTACCGGCCTCGCGAAACACCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGTGTCGACCTTGCTGGCGCCAAAgcgc                                    <  8:201130/134‑1 (MQ=255)
                                                                                                  ccGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACAt                                >  3:267918/1‑137 (MQ=255)
                                                                                                  ccGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACAt                                >  6:30046/1‑137 (MQ=255)
                                                                                                   ccccccAACGGCGCGGTACCGCACTTGCCCAACCCCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACAtt                               <  3:345229/133‑1 (MQ=255)
                                                                                                    gCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTg                              >  6:146824/1‑137 (MQ=255)
                                                                                                      aCCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGAt                            >  3:237162/1‑137 (MQ=255)
                                                                                                      aCCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGAt                            >  7:299274/1‑137 (MQ=255)
                                                                                                      aCCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGAt                            >  7:266054/1‑137 (MQ=255)
                                                                                                       ccAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATc                           >  7:320197/1‑137 (MQ=255)
                                                                                                       ccAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATc                           >  3:211265/1‑137 (MQ=255)
                                                                                                        cAACGGCTCGGTACCGGGCTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCa                          >  1:21477/1‑137 (MQ=255)
                                                                                                        cAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCa                          >  6:186012/1‑137 (MQ=255)
                                                                                                        cAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCa                          >  1:234349/1‑137 (MQ=255)
                                                                                                                   ccccGCACTCGCCCAACACCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCggg               <  1:162266/135‑1 (MQ=255)
                                                                                                                    cccGGACTTGCCCAACACCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGc              <  3:181412/136‑1 (MQ=255)
                                                                                                                        gactgccgcAACACCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACtt          <  2:161283/131‑1 (MQ=255)
                                                                                                                        gACTTGCCCAACCCCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACtt          <  8:27891/137‑1 (MQ=255)
                                                                                                                        caattGCCCACCACCCCACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACtt          <  1:262315/134‑1 (MQ=255)
                                                                                                                                cAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCGTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTTCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACTTGCGCAggg  >  3:137240/1‑137 (MQ=255)
                                                                                                                                cAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACTTGCGCAggg  >  7:272554/1‑137 (MQ=255)
                                                                                                                                cAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACTTGCGCAggg  >  8:112712/1‑137 (MQ=255)
                                                                                                                                cAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACTTGCGCAggg  >  8:69749/1‑137 (MQ=255)
                                                                                                                                cAACAGCACACGACCCCGCCCCGCCAGCGCAACGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGCCCAAGCTGGCGCCAAAGAGCACAATGATCCGCACCTCCGGGCACTTGCGCAggg  >  1:150096/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                   
GCTGGCAAGCGCCTTCAGACACAAACTTCGCCGACCAGTTTGGCCAGCGCTTCAGCATGGGCCCGCACCTGGCTTTCGTCCTCGCCCTCAACCATGACCCGCACCAACGGCTCGGTACCGGACTTGCGCAACAGCACACGACCCCGCCCCGCCAGCGCTTCGGTAACCTTGGCACTGGCTTCCTTGACGGCGGGGTGCTCCAGCGGGTCGACCTTGCTGGCGCCAAAGCGCACATTGATCAGCACCTGCGGGCACTTGCGCAGGG  >  ACB122/5356379‑5356643

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