Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,025 A→G 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,0250AG97.4% 136.2 / NA 38L48L (CTT→CTCU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (9/28);  minor base C (1/0);  total (10/28)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAG  >  CP006881/3724896‑3725149
                                                                                                                                  |                                                                                                                             
cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTATGCTTTTCTTGCTCATTGACCGCCCTCCGCCATTCATAAAACGAAACGAGTcctgag                                                                                                                      >  8:229833/1‑135 (MQ=255)
  gacaggCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATATAGGGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                                                                    <  7:403301/134‑1 (MQ=255)
   gCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTGTGCGTTTCTTCCTGATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                                                                      >  5:425690/1‑136 (MQ=255)
     ccgCTCTCTCCTGCTTTGGCGATTCCTTTATGGCTTGCAAACTTCTTATTTTTTAATGCGCGGGGGCTTTCGCTTTCTGCTTTTATACCTCATTCACAGCCATCGGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                                                                    <  4:42257/134‑1 (MQ=255)
     ccgCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTGGCGTTGTGCGGTTCTTCCTCATTGCCCGCCCGCCGCCCTTCATAAAACGCAACCCGTCCAGCGGc                                                                                                                    >  3:46887/3‑136 (MQ=255)
     ccgCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGGTTTTCTTCCTCATTCAAAGCCATCCGCCATTCATTAAACGAAACGAGTCCTGCGGc                                                                                                                    >  3:345412/3‑136 (MQ=255)
             tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagc                                                                                                           <  1:129988/137‑1 (MQ=255)
               cTGCTTTGGCGATTCCTTGAAGGATTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                                                                       <  2:299038/139‑2 (MQ=255)
                  cTTTGGCGATTCCTTGAAGGCTTGCAAACATCAACTTTAATAAAGCGCGGCGGCTTTCGCGTTAATGTTTTATTACTCAATCACACCCCTCCCCCAGTCATAAAACGAAACGAGTCCTGCGGCCCGCAAcaccgcat                                                                                                      >  8:403301/1‑135 (MQ=255)
                       gCGATTCCTTGAAGGGTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCCTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                <  7:416567/138‑1 (MQ=255)
                                 gAAGGCTTGAAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGt                                                                                         <  3:256480/135‑1 (MQ=255)
                                  aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                    <  7:21900/139‑1 (MQ=255)
                                  aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCCTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                    <  8:43493/139‑1 (MQ=255)
                                    cgCTTGAAACCTTCTTCTTTATTAACGGGGGGCGGCTTTGGATTTGTGCTTTGCTTCCTCATTCCCCCCCACCCCCATTTCATAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                 <  1:101369/138‑1 (MQ=255)
                                     gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                    <  5:61767/136‑1 (MQ=255)
                                     gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGa                                                                                  <  2:69459/138‑1 (MQ=255)
                                          cAAACAGCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGcc                                                                               <  5:359578/136‑1 (MQ=255)
                                                          taACGCGCGGGGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCTTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACTGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAg                                                               <  5:435830/136‑1 (MQ=255)
                                                                       cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                              <  2:48684/139‑1 (MQ=255)
                                                                       cTTTCGCTTTCTGCTTTTCTTCCTCAATCACAGCCATCCGCCATTCATAAACCGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCt                                                  <  5:257692/136‑1 (MQ=255)
                                                                                   gcttttCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCCCCGCCGCCGTTGGTTTGGATGTTCGCGGAGGCGGATAGACCCCGGCCGCCCCCGGTTTGTGGTCCGGCggccc                                     >  5:435630/3‑132 (MQ=255)
                                                                                          ttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTATTTGGTTTTTTCGCTGAGGCGCATATAACCCGGCAGCCCCAGTTTTTCTGTGCAGCGGCACCCAGtt                               >  6:248871/1‑135 (MQ=255)
                                                                                                aTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaaa                      <  1:114203/138‑1 (MQ=255)
                                                                                                aTTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaa                       <  1:108639/137‑1 (MQ=255)
                                                                                                 ttCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCaaaa                       <  6:425690/136‑1 (MQ=255)
                                                                                                      aGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGAcc                  <  3:104669/136‑1 (MQ=255)
                                                                                                      aGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGAcc                  <  5:270749/136‑1 (MQ=255)
                                                                                                        cactCCGCCATTCATAAAACGAACCGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca             <  8:373078/136‑1 (MQ=255)
                                                                                                         cctCCCCCATTAATAAAACCAAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca             <  1:352659/137‑1 (MQ=255)
                                                                                                           tCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                    >  1:61959/1‑129 (MQ=255)
                                                                                                           tCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                    <  2:61959/129‑1 (MQ=255)
                                                                                                            ccGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTGTTTGGTTTGTACGCTGAGCCGCATAGACCGCGGCAGCCGCAGTTTTTCTGTGCCGGGGGCCCCAGGTGATCCAAAACACCAACCCCAcc         >  7:149195/1‑139 (MQ=255)
                                                                                                                cattCATAAAACGAAACGAGTCCTGCGGCCTGGACAGCCGGCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCTCCCAGTTCATCAAAAAGACCATCCTCACCCGTa     <  8:146262/136‑1 (MQ=255)
                                                                                                                cattCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACAccc        <  4:345412/133‑1 (MQ=255)
                                                                                                                cattCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTa     <  1:232107/136‑1 (MQ=255)
                                                                                                                 attCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCg       <  4:113986/134‑1 (MQ=255)
                                                                                                                      tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAg  <  4:164361/136‑1 (MQ=255)
                                                                                                                       aaaaCGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTTTGTACGCGGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCCGCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAg  >  3:185941/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                             
CGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAG  >  CP006881/3724896‑3725149

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: