Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,094,168 A→C 28.7% F476C (TTT→TGT)  tyrR_2 ← Transcriptional regulatory protein TyrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,094,1680AC28.7% 12.2 / 7.5 14F476C (TTT→TGT) tyrR_2Transcriptional regulatory protein TyrR
Reads supporting (aligned to +/- strand):  ref base A (7/3);  new base C (2/2);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

TCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCTTTTAACCACGCCGCCG  >  ACB122/5094063‑5094275
                                                                                                         |                                                                                                           
tCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCtt                                                                               >  5:70282/1‑136 (MQ=255)
              ggCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAACACTCg                                                                 >  8:93581/1‑136 (MQ=255)
                            gCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGCCTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCg                                                   >  3:186477/1‑136 (MQ=255)
                                        gTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATc                                       >  3:116750/1‑136 (MQ=255)
                                                      gCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAatat                         <  3:45832/136‑1 (MQ=255)
                                                        ccGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCACAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGcc                                                    <  1:86304/107‑1 (MQ=255)
                                                        ccGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCACAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGcc                                                    <  1:88970/107‑1 (MQ=255)
                                                        ccGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCACAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGcc                                                    >  2:86304/1‑107 (MQ=255)
                                                        ccGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCACAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGcc                                                    >  2:88970/1‑107 (MQ=255)
                                                                   gAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCGTTTAAc            >  6:174795/1‑136 (MQ=255)
                                                                    aaaCTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCTTTTAAcc           <  7:93581/136‑1 (MQ=255)
                                                                     aaCTCCCCCATCAAACTTTCCAGAACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCTTTTAACCa          <  1:145579/136‑1 (MQ=255)
                                                                             cATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCTTTTAACCACGccgccg  >  3:170847/1‑136 (MQ=255)
                                                                             cATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCTTTTAACCACGccgccg  >  7:172618/1‑136 (MQ=255)
                                                                                                         |                                                                                                           
TCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCACCGCCTTTTCAAAGCGCCCGACTATCTGCGAAAGATCCCCTTCCAGCGAAAACTCGCCCAACGGTTGCCGCACGCCATAATCCGGCAAACGAATATGCTCGCTTTTAACCACGCCGCCG  >  ACB122/5094063‑5094275

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