Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,094,419 A→G 12.1% H392H (CAT→CAC)  tyrR_2 ← Transcriptional regulatory protein TyrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,094,4190AG12.1% 86.8 / 4.7 33H392H (CAT→CAC) tyrR_2Transcriptional regulatory protein TyrR
Reads supporting (aligned to +/- strand):  ref base A (12/17);  new base G (2/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGCACACAACTCGGACAAATCGACCT  >  ACB122/5094292‑5094549
                                                                                                                               |                                                                                                                                  
ccTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCAc                                                                                                                            <  8:46328/136‑1 (MQ=255)
ccTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCAc                                                                                                                            <  2:149539/136‑1 (MQ=255)
ccTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCAc                                                                                                                            <  2:3298/136‑1 (MQ=255)
ccTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCAc                                                                                                                            <  3:42078/136‑1 (MQ=255)
 cTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCAcc                                                                                                                           <  6:70282/136‑1 (MQ=255)
 cTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCAcc                                                                                                                           <  5:174795/136‑1 (MQ=255)
        aaaaCGTATTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCt                                                                                                                    <  4:147299/136‑1 (MQ=255)
          aaCGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTcc                                                                                                                  <  6:89241/136‑1 (MQ=255)
             gTTTTCAAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCAATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCACCCCCTCCAgc                                                                                                               >  5:190767/1‑135 (MQ=255)
                  ccAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCTGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGAccatc                                                                                                          >  7:135035/1‑136 (MQ=255)
                  ccAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGTTGCACCAACCCCTCCAGAccatc                                                                                                          >  6:152804/1‑136 (MQ=255)
                            tACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGAATTcg                                                                                                >  8:133920/1‑136 (MQ=255)
                                      ggCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCACCCCCTCCAGACCATCCATGCATTCGCGCAAAGGCg                                                                                      >  6:121019/1‑136 (MQ=255)
                                      ggCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCTAAGGCg                                                                                      >  7:171234/1‑136 (MQ=255)
                                      ggCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCg                                                                                      >  1:57203/1‑136 (MQ=255)
                                      ggCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCg                                                                                      >  1:158882/1‑136 (MQ=255)
                                      ggCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCAGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCg                                                                                      <  6:118366/136‑1 (MQ=255)
                                           gTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAGTGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAAtg                                                                                 >  7:6150/1‑136 (MQ=255)
                                           gTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAGTGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAAtg                                                                                 >  7:6149/1‑136 (MQ=255)
                                           gTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCCCCACCCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAAtg                                                                                 >  5:17180/1‑136 (MQ=255)
                                           gTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAAtg                                                                                 <  8:169563/136‑1 (MQ=255)
                                             gATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAAtgtg                                                                               <  4:186477/136‑1 (MQ=255)
                                              aTACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGGTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAAtgtg                                                                               <  5:152804/135‑1 (MQ=255)
                                               tACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCa                                                                             <  2:158882/136‑1 (MQ=255)
                                                    cccAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCg                                                                          <  2:57203/134‑1 (MQ=255)
                                                       aaGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAg                                                                      >  6:138193/1‑135 (MQ=255)
                                                            tGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACAtt                                                                >  6:148936/1‑136 (MQ=255)
                                                                                            gcgcAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCg                                 <  5:148936/135‑1 (MQ=255)
                                                                                                   cAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGc                         <  5:138193/136‑1 (MQ=255)
                                                                                                               tGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGCACACAACTCGGa             >  8:125150/1‑136 (MQ=255)
                                                                                                                 gCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCTAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGCACACAACTCGGAc            <  8:171234/135‑1 (MQ=255)
                                                                                                                          agaAAGTGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGCACACAACTCGGACAAATCGACCt  <  8:6149/136‑1 (MQ=255)
                                                                                                                          agaAAGTGCGTCACCAACCCATCCAGAGAATCCATGCATTCGCGCAACGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGCACACAACTCGGACAAATCGACCt  <  8:6150/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                                  
CCTGGAACAAAACGTTTTCCAACTGCCTTACATTACCTGGCCAGTGATACTGCCCAAGTTTGTCCATTGCCGCCGGTGCCAGGCGCGGCATGGCGCAGCCAATCTGCCGACTGGCCTGATCGAGAAAATGCTGCACCAACCCCTCCAGACCATCCATGCATTCGCGCAAAGGCGGAATGTGCAGCGACAGTACATTGAGCCGGTGATACAGGTCCTGGCGAAACTCGCCACGGGCACACAACTCGGACAAATCGACCT  >  ACB122/5094292‑5094549

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