Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,031 C→T A46A (GCG→GCA)  U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,0310CT97.1% 118.8 / NA 35A46A (GCG→GCA) U712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (1/33);  minor base A (1/0);  total (2/33)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATA  >  CP006881/3724896‑3725152
                                                                                                                                       |                                                                                                                         
cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGGTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                                        >  1:35259/1‑139 (MQ=255)
  ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                                                                      <  6:76013/139‑1 (MQ=255)
        ttctctcTTGCTTTGGGGATTCTTTGATGTTTTGCAAAAATCTTCTTTTATAATGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATCCACTGCCATCCGCCTTTCATAAAACGAAACGAGTCCTGCGGCCTGCAc                                                                                                                <  4:92805/138‑1 (MQ=255)
             tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagccg                                                                                                           >  5:373973/1‑139 (MQ=255)
               cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCATCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                                                                         <  8:89524/139‑2 (MQ=255)
               cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTATTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                                                                         <  1:129698/139‑2 (MQ=255)
               cTGCTTTGGAGATTCCTTTAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGATTTCTGCTTTTCTTCCTCCTTCACAGCCATCAGCCATTCATAAAACGAAACGATTCCTGCGGCCTGCACagccgcc                                                                                                         <  1:55722/139‑2 (MQ=255)
                     tGGCGATTCCTTGACGGCTTGCAAACATCTTCTTTTAAAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCCCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                    <  1:161658/138‑1 (MQ=255)
                     tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                    <  5:120922/138‑1 (MQ=255)
                      ggggATTCCTTGAAGGCTTGCAAACATTTTCTTTTAAAACGCGCGGCGGCTTTCGCTTTCTGGTTTTCTTCCTCATTCAAAGCCATCCGCAATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                  <  2:122894/136‑1 (MQ=255)
                      ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                  <  3:32726/139‑1 (MQ=255)
                      ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                    <  5:276299/137‑1 (MQ=255)
                       gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                  <  6:373973/138‑1 (MQ=255)
                       gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTAACAGCCCTCCGCCATTAATAAAACGAAACGAGTACAGCGGCCTGCACACACGCCGTTCTTTg                                                                                                 >  4:175954/1‑139 (MQ=255)
                       agacagccTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTg                                                                                                 <  3:132888/133‑1 (MQ=255)
                            tCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCTGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTAt                                                                                             <  7:91940/138‑1 (MQ=255)
                                 gacggCTTGCAAACTGATCATTTTTTCACGCGGGGCGGCTTTCGCTTTCTGCTTTTCTACCTCAATCACAGCCACCCGCCATACATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                        <  2:391580/135‑1 (MQ=255)
                                 gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                        <  1:337114/138‑1 (MQ=255)
                                 gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                        <  4:112540/138‑1 (MQ=255)
                                 gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                        <  3:91949/138‑1 (MQ=255)
                                  aaGGCTTGCAACCATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                      <  2:40568/139‑1 (MQ=255)
                                   aGGCTTGCAAACATCTTCTTTTATAACGCGCGGCGGGTTTCGCTTTCTGCTTTTCTTCCTCATTCACTGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                      <  1:329508/138‑1 (MQ=255)
                                     gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                   <  7:257657/139‑1 (MQ=255)
                                     gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                   <  2:321461/139‑1 (MQ=255)
                                     gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                   <  6:304222/139‑1 (MQ=255)
                                                         ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc                                                                <  3:397709/138‑1 (MQ=255)
                                                         ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc                                                                <  1:89714/138‑1 (MQ=255)
                                                                  ggcggcTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                                      <  6:160304/139‑1 (MQ=255)
                                                                       cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAACCGCCAATCATAAACCGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                                 <  5:212173/139‑1 (MQ=255)
                                                                                 ctgcttttCTTCCTCATTCACAGCCATCCGCAATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCAc                                        <  1:105502/134‑1 (MQ=255)
                                                                                        tcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTACTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTtc                                 <  5:371913/138‑1 (MQ=255)
                                                                                                aTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAg                        <  4:209330/139‑1 (MQ=255)
                                                                                                  tCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCTCAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                       <  7:93616/138‑1 (MQ=255)
                                                                                                        ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                <  7:301646/139‑1 (MQ=255)
                                                                                                        ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                <  1:56701/139‑1 (MQ=255)
                                                                                                                      tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATa  <  3:175954/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                         
CGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATA  >  CP006881/3724896‑3725152

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