Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,977 T→C 100% 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,9770TC98.0% 85.2 / ‑6.5 31Q64Q (CAA→CAGU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (6/23);  minor base A (1/1);  total (7/24)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724851‑3725105
                                                                                                                              |                                                                                                                                
cTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTcc                                                                                                                      <  7:95271/139‑1 (MQ=255)
                tACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCGTTCTGGTTTTCTTCCGCATTAACAGCCCTcc                                                                                                      >  8:67313/1‑139 (MQ=255)
                   ctcctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGcc                                                                                                   <  4:236843/139‑1 (MQ=255)
                     cctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTCTGCTTTTATTACACATTAACAGCCATCCGCCa                                                                                                  >  6:228921/1‑138 (MQ=255)
                     cctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGCGTTCTGCTTTTATTCCTCATTCACAGCCATCCGCCat                                                                                                 >  2:184060/1‑138 (MQ=255)
                             ccGCCATGTAAAGAGGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATaaaa                                                                                         <  5:191012/139‑1 (MQ=255)
                                ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTAGCTTTATGGTTTTAATCCTGATTAACAGCCAACAGCCATTCaaaaaaag                                                                                       >  1:2628/1‑136 (MQ=255)
                                   tGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                   <  8:59161/139‑1 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  1:184060/139‑2 (MQ=255)
                                               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                       >  8:184179/1‑139 (MQ=255)
                                               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                       <  7:67313/139‑1 (MQ=255)
                                                        tctcCTGCTTTGGCTTTTCATTGAAGGCTTGCAAACATCTTCTTTTATAACGCGCGGGGGCTTGCGCTTTCTGCTTTTCTTCCTAATTCACAGCCATCCGCAATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagc                                                              <  6:36810/139‑1 (MQ=255)
                                                            gtGCGTTGGCGATTCATTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTATGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagccg                                                            <  6:160812/136‑1 (MQ=255)
                                                               cTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGGGTTCTTCGTTTGGTACTCATTAAAAACCATACCCCATTAATAAAAAGAAAAGAGTCCTGGCGCCGGCAAcaccgccgtg                                                       >  1:208989/1‑135 (MQ=255)
                                                                  tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCGTTTCTTCCTCATTCACAGCCCTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCtt                                                    >  5:120890/1‑139 (MQ=255)
                                                                   ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  3:212744/139‑1 (MQ=255)
                                                                   ggCGATTACTTTAAGGCTTGCAAACATCTTCTTTTATAACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACCGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  6:183007/139‑1 (MQ=255)
                                                                    gCGTTTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTg                                                  <  6:151360/139‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  7:184179/138‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  5:276360/138‑1 (MQ=255)
                                                                              gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                         <  5:228921/138‑1 (MQ=255)
                                                                               aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  3:40567/139‑1 (MQ=255)
                                                                               aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  7:241712/139‑1 (MQ=255)
                                                                                aGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCTGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  2:2628/138‑1 (MQ=255)
                                                                                  gCTTGCAAACATCTTCTTTTATTAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCTTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  7:136095/139‑1 (MQ=255)
                                                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCTTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  8:275274/139‑1 (MQ=255)
                                                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  2:253793/139‑1 (MQ=255)
                                                                                                  ttttATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGc                    <  6:73877/139‑1 (MQ=255)
                                                                                                          aGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAg            <  7:54206/139‑1 (MQ=255)
                                                                                                          aGCGCGGCGGCGTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAGAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAg            <  8:203005/139‑1 (MQ=255)
                                                                                                                    cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCAGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc  <  2:208989/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                
CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724851‑3725105

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: