Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,010 C→T 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,0100CT97.7% 131.1 / NA 44E53E (GAG→GAAU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/30);  minor base A (1/0);  total (14/30)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACA  >  CP006881/3724886‑3725140
                                                                                                                            |                                                                                                                                  
tGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                      <  6:137583/139‑1 (MQ=255)
 ctAAAGAGCCGGGCCCGCGCTCGCCTGCTTTGGCGATTCCTTGAAGGGGTCCAAACATCTTATTTTATAAAGCGCGGCGGCTTTCTCTTTCTGCTTTTCTTCATCATTCACCGCAATCCGCCATTAATAAAACGAAAc                                                                                                                      <  6:45254/137‑1 (MQ=255)
       aTGCGGGCCCGCTCTCTCCTGCTTTGGGGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                               <  3:372145/139‑2 (MQ=255)
       aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                               <  5:210739/139‑2 (MQ=255)
       aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTACACATTAACAGCCATCCGCCATTCATAAAACCAAACGAGTcct                                                                                                               >  6:210739/1‑138 (MQ=255)
       aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTGAGGTTTTGTTACGAATTAACACACCTCCGCCATTCATAAAAAGAAACGAGGcct                                                                                                               >  4:372145/1‑138 (MQ=255)
          cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTGGCGTTCTGGTTTTCTTGCTAATTACCAGCCATCCGCCATTCATAAAAAGAAACGAGTCCagcg                                                                                                            >  6:107516/1‑139 (MQ=255)
          cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTGGCGTTCAGGTTTTGTTGCTCAATCACAGCCAGCCGCCCATCATAAAACGAAACGAGTccttcg                                                                                                            >  1:22522/1‑135 (MQ=255)
          cGGGCCCGCGCTCTCCTGCTTGGGCGATTCCTTGACGGCTTGAAAACATCTTCTTTTATAAAGCGCGGCGGCTTGCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                            <  6:283662/139‑1 (MQ=255)
            ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCGTTCTGCGTTTATTGCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                                                          >  1:42360/1‑139 (MQ=255)
                  ctctctCCTGCTTTGGAGATTCCTGGAATGCTTGCAAACATCTTCTTTTATAATGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACATCAATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAc                                                                                                    <  8:315908/139‑1 (MQ=255)
                       tcCTGCTTTGGCGCTTCCTTGAAGGCTTGCAAACTTCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagc                                                                                                 <  3:22270/137‑1 (MQ=255)
                         cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                                                             <  3:357353/139‑2 (MQ=255)
                            cTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCGTTTCTTCCTCCTTAACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCCccgcccccggt                                                                                          >  7:453559/1‑133 (MQ=255)
                                 gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                      >  6:339276/1‑106 (MQ=255)
                                 gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                      <  5:339276/106‑1 (MQ=255)
                                 gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                      <  6:214836/138‑1 (MQ=255)
                                     ttCCTTGAAGGCTTGCAAACATCCTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTTGTAt                                                                                 <  5:194938/139‑1 (MQ=255)
                                           gAAGGCTTGCAAACATCTTCTTTTATAATGCGGGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                            <  5:424822/138‑1 (MQ=255)
                                           catggCTTGAAAACATCTTTTTTTATAAAGCGCGCCGGATTTCGCTTTCTGCTTTTATTCCTCATTCACACCCATCCGCCATTCAAAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                            <  6:66428/135‑1 (MQ=255)
                                            aaGGCTTGCAAACATCTTCTTTTATAATGCGCGGCGGCTTGCGCTTTCTGGTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                          <  2:162974/139‑1 (MQ=255)
                                            aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                          <  8:275415/139‑1 (MQ=255)
                                            aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                          <  3:268264/139‑1 (MQ=255)
                                             aGGCTTGCAAACATCTTATTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACACCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                          <  5:420575/138‑1 (MQ=255)
                                              ggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGa                                                                        <  4:193701/139‑1 (MQ=255)
                                                    cAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCGTGCTGCTTTTCTACCGCAATCACAGCCATCCGCCATTCATAACCCGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGaa                                                                  <  5:298388/139‑1 (MQ=255)
                                                                   ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                      <  8:223983/72‑1 (MQ=255)
                                                                   ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                      >  7:223983/1‑72 (MQ=255)
                                                                    taccgcgcGGCGGCTTTCCCTTTCGGATTTTCTTCAACATCCACCCCCCACCGCCTTTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                  <  2:463066/135‑1 (MQ=255)
                                                                         cgcgGCGGCTTTCGCTTTCTGCTTTTCTTGCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTGGTTGGGTTGTTCGCGGCGGCGGATAGGCCCCGGCcgccgccggt                                             >  5:168702/1‑135 (MQ=255)
                                                                            ggcggcTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                          <  2:42360/139‑1 (MQ=255)
                                                                                 cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                     <  2:24819/139‑1 (MQ=255)
                                                                                              cttttCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCAcc                           <  6:179587/135‑1 (MQ=255)
                                                                                                 ttcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTTTGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGtt                      >  2:272758/1‑138 (MQ=255)
                                                                                                      ccTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTGTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCCCCGTTTTTGTGTCCGGGGTCACCCCGTtcatca                 >  2:302820/1‑138 (MQ=255)
                                                                                                            tCACAGCCATCCGCCATTCATAAAACGAAACAAGTCCTGCGGCCTGCACAGCCGCCGTTGTTTGGTATGTACGCGGAGCCGTATAGACCCCGGCAGCCGCAGGTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGAc          >  8:203636/1‑139 (MQ=255)
                                                                                                             cacaGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGAc          <  7:210528/138‑1 (MQ=255)
                                                                                                                  ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGt                                      <  7:252960/105‑1 (MQ=255)
                                                                                                                  ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGt                                      >  8:252960/1‑105 (MQ=255)
                                                                                                                  ccATCCGCCATTCAAAAAACCAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca    <  3:16974/139‑1 (MQ=255)
                                                                                                                  ccAGCCGCCATTCATTAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGATTCGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGGCACCCAGTTCATCAAAAAGACCATCca    <  8:198209/139‑1 (MQ=255)
                                                                                                                    aTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCcaca  <  8:453559/139‑1 (MQ=255)
                                                                                                                     tCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                     >  7:373901/1‑103 (MQ=255)
                                                                                                                     tCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                     <  8:373901/103‑1 (MQ=255)
                                                                                                                            |                                                                                                                                  
TGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACA  >  CP006881/3724886‑3725140

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: