Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,010 C→T 100% E53E (GAG→GAA U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,0100CT99.5% 110.8 / ‑6.9 38E53E (GAG→GAAU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/31);  minor base G (1/0);  total (7/31)
Fisher's exact test for biased strand distribution p-value = 1.62e-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.

CCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCC  >  CP006881/3724883‑3725137
                                                                                                                               |                                                                                                                               
ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATATAGCGCGGCGGCTTTAGCTTTCTGGTTTTATTCCTCATTCACAGGCATCCGCCATTCATAAAACg                                                                                                                       >  7:555/1‑138 (MQ=255)
   tGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTGTGCTTTTCTTCCTGATTCACCGCCATCCGCCATTCATAAAACg                                                                                                                       >  3:212341/1‑135 (MQ=255)
          tggCGGGCCCGCTCTCTCCTGGTTTGGCTATTCCTTGACGGATTGCAAACATATTCTTTAATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCAAAGCCATCCCCCATTCATAAAACGAAACGAGTcct                                                                                                            <  5:231983/137‑2 (MQ=255)
          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                            <  1:120892/139‑2 (MQ=255)
          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                            <  7:291237/139‑2 (MQ=255)
             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTAGCGTTGTGGTTTTATTACTCATTAACAGCCCTCCGCCAGTCATAAAACGAAACGCGTCCagcg                                                                                                         >  7:259916/1‑139 (MQ=255)
             cGGGCCCGATCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATGCGCCTTTCATAAAACGAAACGAGTCCtgcg                                                                                                         <  5:297120/139‑1 (MQ=255)
               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTAGCTTTGTGGGTTTGTTACTCATTCACAGCCAACCGCCCTTAAAAAAAAGCAACGAGTCCAGCCGc                                                                                                       >  2:43202/1‑139 (MQ=255)
                gCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCGTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                         >  3:140551/1‑136 (MQ=255)
                     ctctctCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAc                                                                                                 <  5:167564/139‑1 (MQ=255)
                          tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagc                                                                                              <  1:12368/137‑1 (MQ=255)
                          tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcag                                                                                               <  8:339325/136‑3 (MQ=255)
                            cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTGTTGCTCATTGACAGCCTTCCGCCATTCATTAAACGAAACGAGTCCTGCGGCCTGCAcagccg                                                                                            <  2:11673/137‑1 (MQ=255)
                                  tgggcaTTCCTTGATGGCTTGCAAAAATCTTATTTAATAATCCGCGGCGGCTTTGGCTTTATGCTTTTCTTCCTCATTCACAGCCATCCGCCATTAATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctt                                                                                        <  3:218520/130‑1 (MQ=255)
                                  tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                     <  1:43202/138‑1 (MQ=255)
                                   ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTc                                                                                      <  4:51240/136‑1 (MQ=255)
                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTc                                                                                      <  3:223437/135‑1 (MQ=255)
                                        ttCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTAGCTTCCACATTCACAGCCATCCCCCATTCATAAAACGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTAt                                                                              <  6:129153/139‑1 (MQ=255)
                                              gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                         <  5:264210/138‑1 (MQ=255)
                                               agacagtGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                       <  7:11585/133‑1 (MQ=255)
                                               aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                       <  1:1230/139‑1 (MQ=255)
                                               aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                       <  5:290103/139‑1 (MQ=255)
                                               aaGGCTGGCAAACATCTTATTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                       <  2:220950/139‑1 (MQ=255)
                                                  gCTTGCAAAGATCTTGTTTTATAAAGGGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAGCCGCCATTCAAAACACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                    <  1:135093/139‑1 (MQ=255)
                                                     tGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTATGCTTTTGTTGCTCATTCACAGCCAACCGCCAATAATAAAACGAAACGAGTCCAGCGGCCAGCCCAGCCGCGCTTATTTTGTTTGTTCGCGGAGGCg                                                                 >  6:211013/1‑139 (MQ=255)
                                                                       taACGCGCGGCGGCTTGCGGTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                               <  5:72796/139‑1 (MQ=255)
                                                                               ggcggcTTTCGCTTTCTGCTTTTCTTCCTCATTCCCCGCCCTCCGCCATTCATAAAACGAAACGCGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                       <  6:18213/139‑1 (MQ=255)
                                                                               ggcggcTTGCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                       <  2:291695/139‑1 (MQ=255)
                                                                                    cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCt                                     <  4:140551/136‑1 (MQ=255)
                                                                                         gCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCg                             <  5:211013/139‑1 (MQ=255)
                                                                                                         ccTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                     <  1:142582/67‑1 (MQ=255)
                                                                                                         ccTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCt                                                                                     >  2:142582/1‑67 (MQ=255)
                                                                                                             aTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaaa          <  8:259916/138‑1 (MQ=255)
                                                                                                             aTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGGATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaa           <  7:22270/137‑1 (MQ=255)
                                                                                                              ttCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaa           <  4:306015/136‑1 (MQ=255)
                                                                                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCAt    <  3:314848/136‑1 (MQ=255)
                                                                                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCAt    <  3:26036/136‑1 (MQ=255)
                                                                                                                       aTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATcc  <  3:353639/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                               
CCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCC  >  CP006881/3724883‑3725137

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: