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,0100CT95.3% 132.3 / NA 43E53E (GAG→GAA) U712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/33);  minor base A (2/0);  total (10/33)
Fisher's exact test for biased strand distribution p-value = 4.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCA  >  CP006881/3724883‑3725138
                                                                                                                               |                                                                                                                                
ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGGTTTATGGTTTTCTTGCGCCTTCACAGCCCTCCGCCATTAATAAAACg                                                                                                                        >  1:90607/1‑138 (MQ=255)
ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACg                                                                                                                        >  5:361410/1‑138 (MQ=255)
   tttAAAGATGCGGGCCCGCTCTCTCTTGCTTTGGGGATTCCTTGAAGGCTTTCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCATTCCGCCATTCATAAAACGAAAc                                                                                                                    <  4:93257/137‑1 (MQ=255)
    gTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTGGCTTTGTGGGTTTGTTGCACATTTACAGCCAAACGCCAATCATAAaaacaaacg                                                                                                                   >  8:141683/1‑138 (MQ=255)
     aaacgATCCGGGCCCGCTCAGTCCTGGTTTGGCGATTACTTGAAGGCTTCAAAAAATTTTATTTTATAAAGCGCGGCGGATTGCGCTTTCTGCTTTTCTTCCTAATACACCGCCAACCCCAATTAATAAAACCAAAcga                                                                                                                  <  4:195278/135‑3 (MQ=255)
          attcGGGCTCGCTCTCTGCTGCTTTGTCGATTCCTTGAAGGCTTGCAAACTTCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                             <  6:67908/136‑2 (MQ=255)
          atcggggCCCGCGCTCTCCTGCTTTGGCGATCCCTTGAAGGATTGCAAACATCTTATTTTATAACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAACCGAAACGAGTcct                                                                                                             <  1:172614/135‑2 (MQ=255)
          aggCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                                             <  5:104939/137‑2 (MQ=255)
          aggCGGGCCCGCCGTCCCCGGCTTGGCCGATGCCTTTAAGGCTGGAAAACATCTTCTTTTAAAACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACCGCCATCCGCCATTCAAAAAACCACACGAGTcct                                                                                                             <  2:372757/137‑2 (MQ=255)
             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                    >  7:179570/1‑129 (MQ=255)
             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                                                    <  8:179570/129‑1 (MQ=255)
             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                          <  6:163830/139‑1 (MQ=255)
             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTGGCGTTATGGGTTTGTTCCGCAATCAAAGCCCTCCGCCATTCATAAAACGAAAAGAGTCCtgcg                                                                                                          >  5:83426/1‑139 (MQ=255)
               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                                                        >  1:24813/1‑139 (MQ=255)
               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTATGCTTTTCTTACTCATTACCCGCCAAACGCCAATAATAAAACGAAACGAGTccataggc                                                                                                        >  4:155571/1‑133 (MQ=255)
                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                    <  6:361410/138‑1 (MQ=255)
                                     cGATTCCTTGAAGGCTAGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTgg                                                                                  <  8:350680/139‑1 (MQ=255)
                                         tCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTCTGCTTTTCTTCCTCATTAACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCCCAGCCGCCGTTATTTGGtttt                                                                              >  8:34701/1‑136 (MQ=255)
                                              gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                          <  4:156450/138‑1 (MQ=255)
                                              gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                          <  2:255427/138‑1 (MQ=255)
                                               agggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                        <  1:133392/137‑1 (MQ=255)
                                                  gctggaaaaCTTATTCTTTTATAATGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTAATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  1:186127/133‑1 (MQ=255)
                                                  gCTTGCAACCTGTTTCTTTTATAAAGCCGGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCTTTCCCAGCCATCCCCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  6:32921/139‑1 (MQ=255)
                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  5:269517/139‑1 (MQ=255)
                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  8:242938/139‑1 (MQ=255)
                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  5:264193/139‑1 (MQ=255)
                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  2:24813/139‑1 (MQ=255)
                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  7:34701/139‑1 (MQ=255)
                                                  ccTTCCAAACATCTCTTTTTATAACGCGCGGCGGCTTGCGCTTTCGGATTTTCTCCCTCATTCACCGCCATCCGCCATTCAAAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                     <  1:373151/138‑1 (MQ=255)
                                                     tGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCCCCGCCGCCAGTCTTTGGTATGTACGCTGAGcc                                                                   >  7:164225/1‑138 (MQ=255)
                                                                       taacccgcgGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCCCCCCCATCCGCCATTCAAAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                <  3:103901/134‑1 (MQ=255)
                                                                       taAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                <  5:351930/139‑1 (MQ=255)
                                                                            cgcgGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCCCAGCCGCCGTTGTTTGGTTTGTTCGGTGGGCCGGATAGGCCGCGGCCGGCGCAggt                                           >  5:345956/1‑137 (MQ=255)
                                                                               ggcggcTTTCGCTTTCTGCTTTTCTTCCGCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                        <  2:15790/139‑1 (MQ=255)
                                                                                         gCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCg                              <  7:141683/139‑1 (MQ=255)
                                                                                             tctgcttttCTTCCTCATTCACAGCCAACCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCa                           <  6:374533/133‑1 (MQ=255)
                                                                                                     tcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTtc                   <  6:83426/138‑1 (MQ=255)
                                                                                                     tcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTtc                   <  2:90607/138‑1 (MQ=255)
                                                                                                              ttCACAGCCAGCCGCCATTCATAAAAGGAAACTAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTTTGTACGCTGAGCCGAATAGACCGCGGCAGCGGCAGGTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa         <  5:207576/139‑1 (MQ=255)
                                                                                                                     ccATCCGCCATTCATAAAACGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca  <  2:98804/139‑1 (MQ=255)
                                                                                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca  <  7:269131/139‑1 (MQ=255)
                                                                                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca  <  3:155571/139‑1 (MQ=255)
                                                                                                                     ccATCCCCCATTCATAAAACGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca  <  2:275190/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                
CCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCA  >  CP006881/3724883‑3725138

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: ‑