Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,724,977 T→C Q64Q (CAA→CAG U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,9770TC85.7% 62.7 / ‑4.3 21Q64Q (CAA→CAGU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (2/16);  minor base G (3/0);  total (5/16)
Fisher's exact test for biased strand distribution p-value = 7.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTT  >  CP006881/3724851‑3725100
                                                                                                                              |                                                                                                                           
cTTTTGACTGGGATTTTACCTCCTTTTGCCCGCCATGTAAAGATGCGGGCCCGCGCTCTCCTGCTTTGGCGATTCCTTGACGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTttcttc                                                                                                                  <  7:3540/138‑1 (MQ=255)
          ggATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTATCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGGTTTGGGGTGTTGGTGCTCATTcaccg                                                                                                       >  8:44380/1‑137 (MQ=255)
                     cctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGCTTTGTGCTTTTGTTCCTCATTCACAGCCATGCGCCat                                                                                            >  6:111284/1‑138 (MQ=255)
                             ccGCCATGAAAAGTTGGTGGCCCGCTCGCTCGTGCTTTGGCGATTCCTTGAAGGATTGCAAACATCTTCTTTTTTAAAGCGCGGCGGATTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAGGCGCCATTCATaaaa                                                                                    <  4:112707/139‑1 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                       <  2:162391/139‑2 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                       <  8:34509/139‑2 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCACCATTCATAAAACGAAACGAGTcct                                                                       <  5:80437/139‑2 (MQ=255)
                                               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTGGCGTTGTGGTTTTGTTGCTGATTAACAGCCCTCCGCCCATCATAAAACGCAACGCGTACTGCGGc                                                                  >  8:54044/1‑139 (MQ=255)
                                                  ccgCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTg                                                               <  3:32751/137‑1 (MQ=255)
                                                       ctctcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcag                                                          <  7:54044/139‑3 (MQ=255)
                                                          tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCCAAGCCATCCGCCATTCATAAAACGAAACCAGTCCAGCGGCCTGCCccgccc                                                       >  3:217015/1‑138 (MQ=255)
                                                            cTGCTTTGGGGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATATAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGACATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                     <  1:229576/139‑2 (MQ=255)
                                                                   ggCGATGCCTTGATGGCTTGCAAACATCTTCTTTTTTAACGCGCGGCGGATTTCGATTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                              <  4:1375/139‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                              <  4:63746/138‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                              <  6:132074/138‑1 (MQ=255)
                                                                     cGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTgg                                            <  7:44380/139‑1 (MQ=255)
                                                                                      gCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                  >  7:97554/1‑132 (MQ=255)
                                                                                      gCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                  <  8:97554/132‑1 (MQ=255)
                                                                                                      ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc            <  5:274366/138‑1 (MQ=255)
                                                                                                       taAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg          <  2:39611/139‑1 (MQ=255)
                                                                                                               ggcggcTTTCGATTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt  <  4:102576/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                           
CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTT  >  CP006881/3724851‑3725100

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: