Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,375,132 C→G 9.3% A277G (GCC→GGC)  dinB → DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,375,1320CG 9.3% 130.1 / 7.0 43A277G (GCC→GGC) dinBDNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base G (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

GTGCGGTGCATAACGACAGCCGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAATTGCTGGG  >  ACB122/1375003‑1375258
                                                                                                                                 |                                                                                                                              
gTGCGGTGCATAACGACAGCCGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAAt                                                                                                                          <  7:172061/136‑1 (MQ=255)
 ggCGGTGCATAACGACAGCCGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATg                                                                                                                         <  6:90718/135‑1 (MQ=255)
                   ccGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGc                                                                                                         >  4:153099/1‑134 (MQ=255)
                   ccGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGc                                                                                                         >  2:144870/1‑134 (MQ=255)
                   ccGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGc                                                                                                         >  8:168717/1‑134 (MQ=255)
                   ccGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGc                                                                                                         >  5:182549/1‑134 (MQ=255)
                   ccGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGATGCTGGATAGCCTGAACGAACGCCTCGCCCGAATGGACAGCAGGTACCGGc                                                                                                         >  5:166824/1‑134 (MQ=255)
                    cGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAg                                                                                                      <  1:45210/136‑1 (MQ=255)
                     gACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGCCCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGa                                                                                                     >  1:105557/1‑136 (MQ=255)
                     gACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGGCCGAATGGACAGCAGTTACCGGCCAGa                                                                                                     >  4:38851/1‑136 (MQ=255)
                     gACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGGCCGAATGGACAGCAGTTACCGGCCAGa                                                                                                     >  4:36005/1‑136 (MQ=255)
                     gACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGa                                                                                                     >  6:79293/1‑136 (MQ=255)
                                gTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCg                                                                                          <  1:147440/136‑1 (MQ=255)
                                       tgGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGGCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCgtcaaggt                                                                                   <  3:38851/136‑1 (MQ=255)
                                       tgGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGGCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCgtcaaggt                                                                                   <  3:36005/136‑1 (MQ=255)
                                          aaaaCACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCgtcaaggtc                                                                                  >  2:193953/1‑134 (MQ=255)
                                                  tACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCAt                                                                         >  8:39795/1‑135 (MQ=255)
                                                       cacCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACtt                                                                    >  2:15211/1‑135 (MQ=255)
                                                            aCCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCa                                                              >  4:151291/1‑136 (MQ=255)
                                                            aCCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCa                                                              <  7:170492/136‑1 (MQ=255)
                                                            aCCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCa                                                              <  5:35844/136‑1 (MQ=255)
                                                            aCCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCa                                                              <  1:91902/136‑1 (MQ=255)
                                                                 ccAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGacca                                                         <  1:15211/136‑1 (MQ=255)
                                                                 ccAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGacca                                                         <  7:7548/136‑1 (MQ=255)
                                                                 ccAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGacca                                                         <  7:105636/136‑1 (MQ=255)
                                                                      ccTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGCACGCACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGGTCGTCAAGGTCAAGTTCCATGACTTCAGCCCGACCACCATg                                                    <  1:123074/136‑1 (MQ=255)
                                                                            cAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAg                                              >  8:148216/1‑136 (MQ=255)
                                                                                 gCCTGGCGCGGCTGCCAGAGTTGCGGGAGAGCCGGAACGAACGCATCGCCCGAAGGGACAGCTGTAACCGGCCAGACAAACCGCTCGTCAATGTCACGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCg                                          <  1:181928/135‑1 (MQ=255)
                                                                                 gCCGGGCGCGGCTGCCAGAGTGGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCg                                          <  5:36250/135‑1 (MQ=255)
                                                                                       cgcgGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGgggcggg                                   >  3:167814/1‑136 (MQ=255)
                                                                                        gcgGCTGCCAGAGTTGCGGGATAGCCTGCACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGgggcgggc                                  <  1:95123/136‑1 (MQ=255)
                                                                                          ggCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTa                                                                                                              >  1:20115/1‑58 (MQ=255)
                                                                                          ggCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTa                                                                                                              >  5:118964/1‑58 (MQ=255)
                                                                                          ggCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTa                                                                                                              <  2:20115/58‑1 (MQ=255)
                                                                                          ggCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTa                                                                                                              <  6:118964/58‑1 (MQ=255)
                                                                                                         tGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTa                 >  8:34392/1‑136 (MQ=255)
                                                                                                                 cTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAAt         <  5:79293/136‑1 (MQ=255)
                                                                                                                 cTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGGCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAAt         <  3:155579/136‑1 (MQ=255)
                                                                                                                   gaacgaacGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGACAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAATTg       >  1:60691/1‑136 (MQ=255)
                                                                                                                        aacgCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAATTGCTggg  >  4:163959/1‑136 (MQ=255)
                                                                                                                        aacgCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAATTGCTggg  >  8:150973/1‑136 (MQ=255)
                                                                                                                             aTCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAg                                                                 <  7:114801/68‑1 (MQ=255)
                                                                                                                             aTCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAg                                                                 >  8:114801/1‑68 (MQ=255)
                                                                                                                                 |                                                                                                                              
GTGCGGTGCATAACGACAGCCGACGCCAGTCGGTGAGTGTGGAAAACACCTACGACACCGACCTGCCAGACCTGGCCAGTTGCCTGGCGCGGCTGCCAGAGTTGCTGGATAGCCTGAACGAACGCATCGCCCGAATGGACAGCAGTTACCGGCCAGACAAACCGTTCGTCAAGGTCAAGTTCCATGACTTCAGCCAGACCACCATGGAGCAGGCCGGGGCGGGCCGGGACCTGGAGAGTTATCGGCAATTGCTGGG  >  ACB122/1375003‑1375258

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: