Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,123,2710CG34.2% 41.7 / 10.6 38intergenic (‑273/+31)gcvH_1/CJ019_01062Glycine cleavage system H protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/7);  new base G (13/0);  total (31/7)
Fisher's exact test for biased strand distribution p-value = 7.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACGCATCG  >  ACB143/1123145‑1123404
                                                                                                                              |                                                                                                                                     
gtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTAAGACGGAGCGAAAACGATACATTGGGTGGCTGGACTgg                                                                                                                              >  2:100687/1‑136 (MQ=255)
gtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGCCGTAGCGCAAACGATACATGTGGTGCCGGGTCTgg                                                                                                                              >  8:182328/1‑136 (MQ=255)
gtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGCTACAAGTGGTGCCGGGGCTgg                                                                                                                              >  6:357612/1‑136 (MQ=255)
  ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACAAGTGGTGCCTGGTCTGGCg                                                                                                                            >  1:240067/1‑136 (MQ=255)
    aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTTGCGCAAACGATTCATGTGGTGGCGGGCCTGGCGAg                                                                                                                          >  5:306472/1‑136 (MQ=255)
       gATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTAAGACCTAGCGAAAAAGCTACATGTGGTGCCGTGCCTGGCGAGGTg                                                                                                                       >  8:210555/1‑136 (MQ=255)
              ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACCACGTAGCGAAAACCAAACATTTGGTGCCTGGACTGGCGAGATGCCAATgg                                                                                                                >  2:313442/1‑136 (MQ=255)
                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGGACGACGGAGCGAAAACGATACATGTGGTGGCGGGACTGGCGAGGTTGCAATGGccg                                                                                                             >  5:195098/1‑135 (MQ=255)
                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGCCGGTGCGTAAACGATACATTTGGTTCCTGGGCTGGCGGGGTGCCAATGGcc                                                                                                              >  3:45939/1‑135 (MQ=255)
                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGCCGGTGCGAAAAAGCAAAAAGTGGTTGCGGGGCGGGGGCGGTGGCAATGGccg                                                                                                             >  7:322734/1‑135 (MQ=255)
                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTAAGACGTAGCGCAAACGATACATGTGGTGCCTGGTCTGGCGAGGTGCCCATTGccg                                                                                                             >  6:350434/1‑135 (MQ=255)
                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTAAGACGTAGCGAAAACGATACATGTGGTGCCTGGCCTGGCGAGGTGCCAATGGCCt                                                                                                             >  7:271820/1‑136 (MQ=255)
                   ggTGCCTCTGAGAGGAAACGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACATGTGGTGGCTGGACTTGCGAGGTGCCAATTGCCGTc                                                                                                           >  1:332718/1‑136 (MQ=255)
                          ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGCCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc                                                                                                    >  8:360199/1‑136 (MQ=255)
                          ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGGTTTAGCGTCGCCGGAGCGCAAACGATAAAAGTGGTGGCGGGGCGGGGGGGGTGGCAATTGCCGTGGGtga                                                                                                      >  4:162509/1‑132 (MQ=255)
                            gagaGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTTAGACGGAGCGCAAACGATACATTTGGTGCCTGGACTGGCGCGGTGGCAATGGCCTTCgggaatgt                                                                                                   >  5:371688/1‑129 (MQ=255)
                                    aaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGTGACATGTGGTGGCGGGGCTGGCGCGGTGCCAATGGCCTTCGGTGATGTCCCGTTcag                                                                                          >  2:95059/1‑136 (MQ=255)
                                     aGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGTCGACGGTGGGCAAACGATACATGTGGTGGCGGGGCTTGCGGGGTGCCAATTGCCGTCGGTCCTCTCACTTTgag                                                                                          >  8:211838/1‑135 (MQ=255)
                                     aGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTACCAACCCCGGAACGAAAACGCTACAGGGGGGGGCAGGGCAGGCGAGGTGCCAAAGGCCTTAGGTACTCACACTTTcac                                                                                          >  4:304424/1‑134 (MQ=255)
                                      gggTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTTGCGACGACGTTGCGCAAACGAAACATGTGGTGGCTGGGCTGGCGAGGGGGCAATGGGCCTAGGTGCTCTCACGTTcagc                                                                                         >  4:374535/1‑135 (MQ=255)
                                              gCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATTCATGTGGTGGCTGGACTGGGGAGGTGCCAATGGGCTTCGGTGATCTCACTTTGAGCAGATGATc                                                                                 >  8:251988/1‑135 (MQ=255)
                                               cATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGGAAACGATACATGTGGTGCCTGGACTGGCGAGGTTCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCa                                                                               >  1:296086/1‑136 (MQ=255)
                                                             gCACCGACGGCGAGGAGGCTGTGTATCGTTTTTTGTAAGACGTTGCGAAAAAGAGAAATGTTGTGGCAGGACCTGCGAGGGGGCCATGGCCTTAGGTCATATCACATTCAGCAGATTATCCATCGAAAACGccccc                                                                 >  2:233228/1‑135 (MQ=255)
                                                              cACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGCCGTATCGAAAAACATACAAGTGGTGCCTGGACTGGCGAGGTGCCAAAGGCCTTCGGTCATCTCACTTTAAGCAGATGATCAACCGAAAAACCGCCgg                                                                >  8:360899/1‑136 (MQ=255)
                                                                    cGGCGAGGAGGCTGTGTATCGTTTTAGCGAACACGTAGCGCAAACCAAAAATGTGGTTGCTGGGCTGGGGCGGTGCCAATTGCCGTAGGTGATCTCACATTAAGCAGATGATCCAGCGAAAACGCCCCGGCGCCTc                                                          >  2:9605/1‑136 (MQ=255)
                                                                                          tttCGATACGACTTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAaga                                     <  7:251988/135‑1 (MQ=255)
                                                                                                           gAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcaggt                   >  1:15071/1‑136 (MQ=255)
                                                                                                           gAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcaggt                   >  1:262959/1‑136 (MQ=255)
                                                                                                                agATACATGTGTGGCCTGGACTGGCGTGGGGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGTTTATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtgggtg              <  2:38291/135‑1 (MQ=255)
                                                                                                                   tACATGGGGGGCCTGGACTGGCGAGGTGCCATTGGCCTTCGGTCATCTCACTTTCAGCAGATGCTCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGtgggtggg            <  7:301728/135‑1 (MQ=255)
                                                                                                                    acttgtgGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGTTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGtgggtggg            <  3:248721/131‑1 (MQ=255)
                                                                                                                      aTGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCa         >  3:238125/1‑135 (MQ=255)
                                                                                                                      aTGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCAc        >  7:915/1‑136 (MQ=255)
                                                                                                                       tgtgGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACg       <  8:96157/136‑1 (MQ=255)
                                                                                                                        gtggtgCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACg       <  7:182328/135‑1 (MQ=255)
                                                                                                                        gtggtgCCTGGACTGGCGAGGTGCAAATGGCCTGCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACg       <  2:332718/135‑1 (MQ=255)
                                                                                                                           gtgCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACGCATc   >  7:344672/1‑136 (MQ=255)
                                                                                                                            tgCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACGCATCg  >  1:267834/1‑136 (MQ=255)
                                                                                                                              |                                                                                                                                     
GTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGGGCCACGCATCG  >  ACB143/1123145‑1123404

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