Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,025 A→G 100% L48L (CTT→CTC U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,0250AG95.6% 100.1 / ‑5.7 30L48L (CTT→CTCU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (4/24);  minor base C (2/0);  total (6/24)
Fisher's exact test for biased strand distribution p-value = 3.45e-02
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.

ATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGA  >  CP006881/3724893‑3725150
                                                                                                                                    |                                                                                                                             
aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                                         <  1:194577/139‑2 (MQ=255)
   cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                                      <  6:172027/139‑1 (MQ=255)
   cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTAGCTTTCTGGTTTTATTACACATTACCAGCCCAACGCCAATGATAAAACCAAACGAGTCCagcc                                                                                                                      >  6:64368/1‑138 (MQ=255)
     ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTATGCTTTTATTACTCATTCACAGCCATCCGCCAATCATAAAACGAAACCAGTCCTGCGGc                                                                                                                    >  8:52751/1‑139 (MQ=255)
           ctctctCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTATGGTTTTATTACTCATTAACACCCAACCGCCCTTCATAAAAAGAAACCATTCCTGCGGCCAGccc                                                                                                              >  8:162333/1‑137 (MQ=255)
                tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagc                                                                                                           <  3:115503/137‑1 (MQ=255)
                  cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagccg                                                                                                         <  3:3581/137‑1 (MQ=255)
                        tGGCGATTCCTTGATGTCTTGCAAACATCTTATTTTATAAAGCGCGGCGGCTTTCTCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCTTAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                  <  5:11982/138‑1 (MQ=255)
                        tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                  <  7:52751/138‑1 (MQ=255)
                        tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                  <  5:287798/138‑1 (MQ=255)
                          gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAACCGACACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                <  1:274333/138‑1 (MQ=255)
                          gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                <  5:203995/138‑1 (MQ=255)
                          gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                <  5:208236/138‑1 (MQ=255)
                               tCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATg                                                                                          >  7:81861/1‑139 (MQ=255)
                                ccTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATg                                                                                          <  2:77232/138‑1 (MQ=255)
                                    gatggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                      <  5:128284/135‑1 (MQ=255)
                                    gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                      <  3:150048/138‑1 (MQ=255)
                                     aaGGCTGGCAAACATCTTCTTTTATAAAGCGCGCCGGCTTTCGCTTTCTGCTTTTCTTCCTCAATCCCAGCCATCCGCCATTCATAAAACGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                    <  2:77764/139‑1 (MQ=255)
                                          ttGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTATTCCTGATTCACAGCCATCCGGCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGGCGCCCTTCGTTTGTTTATACGCtatacc                                                                               >  3:167996/1‑134 (MQ=255)
                                                          ttATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCgg                                                                  <  7:162333/136‑1 (MQ=255)
                                                             taACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                            <  3:202022/139‑1 (MQ=255)
                                                                          cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                               <  8:111415/139‑1 (MQ=255)
                                                                               gCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCg                                          <  8:65004/139‑1 (MQ=255)
                                                                                          ttcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                                      >  7:229501/1‑52 (MQ=255)
                                                                                          ttcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                                      <  8:229501/52‑1 (MQ=255)
                                                                                                   aTTCACAGCCAACCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaa                        <  4:300132/137‑1 (MQ=255)
                                                                                                    ttCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                     <  5:64368/139‑1 (MQ=255)
                                                                                                    ttCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                     <  3:231437/139‑1 (MQ=255)
                                                                                                    ttCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                     <  8:81861/139‑1 (MQ=255)
                                                                                                                        aTAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGa  <  5:288088/138‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
ATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGA  >  CP006881/3724893‑3725150

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: