Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 1,835,872 T→C S660G (AGC→GGC)  barA ← Signal transduction histidine‑protein kinaseBarA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,835,8720TC75.0% 10.3 / 5.3 18S660G (AGC→GGC) barASignal transduction histidine‑protein kinaseBarA
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (0/12);  total (3/15)
Fisher's exact test for biased strand distribution p-value = 7.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-02

CATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTACGGGTGGTACAGCGCCTGCTCGGTGGT  >  ACB122/1835746‑1836007
                                                                                                                               |                                                                                                                                       
cATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTGGCGGTGGCGTTTGGGTGGGTcccccgggg                                                                                                                                >  2:196239/1‑132 (MQ=255)
 aTAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGGCGGTGGGGTTGGTCCCCACCGccc                                                                                                                              >  5:1059/1‑136 (MQ=255)
                  ccAGTACCTCGGCACCCAGGTCTTCGAGCAGGTTTTTTACCAGCAGCAGGTTGGCGGCTTTGTCGTCGACGCAGAGTATTTTCGGTCGGCGCTGGCCGTTGGGGTGGCCCCCCTCGGCCAGCGGACGGCGAGGTTGc                                                                                                             <  7:266554/137‑1 (MQ=255)
                                        tgcGAGCCGGGTTTTACCCAGCAGCAGGTGGGCGCCATGGGCGTCCCCGCAGTGTATTTGCGGGCGGCGCGGGCCGGTGGGGGGCCCCCCCGCGCCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGcgcc                                                                                       <  2:204359/135‑1 (MQ=255)
                                                       aaccaccaccTGGTGGGCGGATTTGGCGGCGACGCAGTGTTTTTTCGGGCGGCGGTGGCGTTGGGGGTGGCGCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGc                                                                        <  8:19427/136‑1 (MQ=255)
                                                            cccgcagGTGGGCGGCTTTGGCGTCGCCGCTGTTTTTTTTGGGGCGGGGCGGGCGTTGGGGGTGGCCCCCCCCGGCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTc                                                                   <  4:110537/134‑1 (MQ=255)
                                                             cagccgGTGGGCGGCATTGGCGTCGACGCAGCTTATTTTCGGCCGGCGGTGGCCGTTGGGGTGCCCCCCCGCGGCCAGGGGACGGCGAGGGTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCg                                                                  <  6:336253/137‑1 (MQ=255)
                                                                               gtggtcgACGCAGTTTATTTGGGGGGGTCGGTGGCGGTTGGGGTGGCCCCCCGCGGCCTGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGc                                                <  3:252678/134‑1 (MQ=255)
                                                                               gtcgtcGACGCAGAGTATTTTCGGGCGGCGCGGGCCGTTGGGGTGGCGCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGc                                                <  2:264112/137‑1 (MQ=255)
                                                                                    cgccGCTGAGTTTTTTCGGGCGGGGGTGGCGGTTGGGGTGGCCCCCCGCGGCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATgg                                           <  4:102901/134‑1 (MQ=255)
                                                                                         cAGAGTTTTTTCGGGCGTCGCTTTCCGTTGGGGGTGCCCCCCGCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGtt                                      <  6:1059/137‑1 (MQ=255)
                                                                                              aaTTTTCGGGCGGCGCTGGCCGTGGGGGTTGACCCCCTCGGCCAGCGGACGGCGAGGTTGCACCATCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCa                                 <  8:100705/136‑1 (MQ=255)
                                                                                                   tCGGTCGTCGCTGACCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTAc                            <  8:120756/137‑1 (MQ=255)
                                                                                                   ggggTCGGCGCGGGCGGTGGGGGTGCCCCCCCGGGCCCCGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTAc                            <  3:179015/135‑1 (MQ=255)
                                                                                                   gcGGTGGGCGCGGGCCGTTGGGGTGGCCCCCCGCGGCCGGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGCGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTat                            <  8:115893/136‑2 (MQ=255)
                                                                                                   gcGGGCGGCGCGGGCCGTTGGGGTGGCCCCCCTCCGCCAGCGCACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTAc                            <  6:279082/136‑1 (MQ=255)
                                                                                                                       ggggGTGCCCCCGGCGCCCCGCGGACGGCGAGGGTGCACCAGCTCCAGCAAGAGGCGCCGCAGCGTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTACGGGTGGTACAGCGCCTGCt         <  8:179407/137‑1 (MQ=255)
                                                                                                                              cTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTACGGGTGGTACAGCGCCTGCTCggtggt  >  5:273261/1‑137 (MQ=255)
                                                                                                                               |                                                                                                                                       
CATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGGCCATGGCCGTTGGGCAGGTACGGGTGGTACAGCGCCTGCTCGGTGGT  >  ACB122/1835746‑1836007

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: