Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,835,8490TG52.6% 7.3 / 19.0 19R667R (CGA→CGC) barASignal transduction histidine‑protein kinaseBarA
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base G (0/10);  total (4/15)
Fisher's exact test for biased strand distribution p-value = 3.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGG  >  ACB122/1835713‑1835978
                                                                                                                                        |                                                                                                                                 
cGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgt                                                                                                                                   >  5:342796/1‑137 (MQ=255)
cGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgt                                                                                                                                   >  8:36596/1‑137 (MQ=255)
cGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgt                                                                                                                                   >  2:435558/1‑137 (MQ=255)
    cAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               <  6:319251/137‑1 (MQ=255)
    cAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               <  3:210335/137‑1 (MQ=255)
    cAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGGTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGTGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               <  4:127371/137‑1 (MQ=255)
     aaTTCTTCCTGCACGGCCTGTACCGCCGCATACCCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGGGGCATTGTCGTCGACGCAGAGTATTTTCGGgcgtcg                                                                                                                              <  7:175623/137‑1 (MQ=255)
     aaTTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGtcgtcg                                                                                                                              >  6:320630/1‑137 (MQ=255)
     aaTTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACAAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtgg                                                                                                                              >  7:375130/1‑135 (MQ=255)
     aaTTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCAACCAGGTCTTAGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGtcgtcg                                                                                                                              >  2:294752/1‑137 (MQ=255)
            ccTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCGTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGCCGCAGAGTATTTTCGGTCGTGGCGGGcgg                                                                                                                       >  8:95100/1‑135 (MQ=255)
                  cGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTgggg                                                                                                                 <  5:230589/137‑1 (MQ=255)
                   ggCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGt                                                                                                                <  4:35024/137‑1 (MQ=255)
                      cTGTACCGCCGCATAACCATTGTAGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGc                                                                                                             <  5:11885/137‑1 (MQ=255)
                       tGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTGGCGGGGCGGTGGGGGTggt                                                                                                            >  2:6407/1‑135 (MQ=255)
                        gTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTGGCGGGGCGGTggggggtgtg                                                                                                           >  8:368699/1‑131 (MQ=255)
                        gTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTc                                                                                                           <  2:88677/137‑1 (MQ=255)
                                gcATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               >  6:422066/1‑109 (MQ=255)
                                gcATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               <  5:422066/109‑1 (MQ=255)
                                  aTAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTGGCGGGGCGGGGGGGGTGGGCCCCgccgccc                                                                                                 >  8:229049/1‑136 (MQ=255)
                                                    cAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               <  2:40666/89‑1 (MQ=255)
                                                    cAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtcgtc                                                                                                                               >  1:40666/1‑89 (MQ=255)
                                                                          tCGACCCGGGTTTTAACCAGCCGCAGGTGGGGGGCATGGTCGGCGACGCAGAGTTTTTTCGGGCGGCGCTGGCCGTGGGGGTGGCCCCCCGCGGCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGcgccg                                                         <  4:366051/137‑1 (MQ=255)
                                                                               ccgggTTTTAACCAGCCCCAGGTGGGCGGCATTGCGGGCGACGCAGAGATTTTTCGGGCGGCGGTGGCCGTTGGGGTGGCCCCCCTCGGCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCt                                                    <  6:323066/135‑1 (MQ=255)
                                                                                                        ggggcATGGTCGGCCCCGCAGCGTTTTTGGGGGCGGCGCGGGCCGTTGGGGTTGCCCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAg                           <  6:341277/135‑1 (MQ=255)
                                                                                                          ggaaTTGTCGTCGACGCTCATTTTTTTCGGGCGGCGGTGGCCGTTGGGTTGGCGCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCa                         <  6:440900/134‑1 (MQ=255)
                                                                                                               tggcgtcgCCGCCGTTTATTTTGGGGCGGCGCGGGCCGTTGGGGTTGCCCCCCGCGCCCAGGGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTgg                    <  7:229049/134‑1 (MQ=255)
                                                                                                                gtcgtcGACGCAGAGTATTTTCGGTCGTCGCGGGCGGTTGGGGTGGCGCCCCTCGGCCAGCGGTCGGCGAGGTTGCCCCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGc                   <  7:368699/137‑1 (MQ=255)
                                                                                                                          ctgtttATTTTCGGGCGGCGGGGGCGGGGGGGGTGGCCCCCCTCGGCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCGTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGtt         <  7:315219/132‑1 (MQ=255)
                                                                                                                          aagtttATTTTCGGGGGGCGCGGGCGGTTGGGGTGGCCCCCCGGGGCCGGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGtt         <  1:6407/132‑1 (MQ=255)
                                                                                                                               ttttttGGGGCGGCGCGGGCGGGTGGGGGGGCTCCCCTCGCCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCa    <  7:347854/135‑1 (MQ=255)
                                                                                                                               tATTTGCGGGGGTCGGTGGCCGTGGGGGTGGCCCCCCGCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCa    <  5:320630/137‑1 (MQ=255)
                                                                                                                                 tttgCGGGCGTCGGGGGCGGTGGGGGTTGCCCCCCTCGCCCAGGGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAgg  <  1:291377/137‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
CGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGG  >  ACB122/1835713‑1835978

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