Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	3,724,980	T→C	K63K (AAA→AAG)	U712_19205 ←	Uncharacterized protein ywcI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,724,980	0	T	C	88.6%	129.2 / ‑3.6	44	K63K (AAA→AAG)	U712_19205	Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (5/34);  minor base G (4/1);  total (9/35)
Fisher's exact test for biased strand distribution p-value = 4.18e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724846‑3725105                                                                                                                                        |                                                                                                                              aTTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGctttt                                                                                                                           <  8:381768/139‑1 (MQ=255)      cTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTttcttc                                                                                                                       <  5:297308/138‑1 (MQ=255)      cTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTcc                                                                                                                      <  8:26091/139‑1 (MQ=255)         ttGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCATCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCt                                                                                                                     <  2:141281/137‑1 (MQ=255)             cTGTGATTTTACGTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGGTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTca                                                                                                               <  5:304782/139‑1 (MQ=255)                         ctcctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGcc                                                                                                   <  4:259453/139‑1 (MQ=255)                            ctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTGGCTTTCTTGTTTTCTTGCTCATTCAAAGCCATACGGCagt                                                                                                >  3:374760/1‑137 (MQ=255)                            ctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCat                                                                                                 >  8:257657/1‑137 (MQ=255)                                  aCCGCCATGTAACGTTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATaaa                                                                                          <  5:297972/139‑1 (MQ=255)                                  aCCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGCGTTATGGTTTTGATACTCATTAAACGCCCACCGCCCATCATaa                                                                                           >  5:368300/1‑138 (MQ=255)                                    agCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAAc                                                                                        <  7:16525/138‑1 (MQ=255)                                      ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGa                                                                                      <  6:273991/139‑1 (MQ=255)                                      ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCGGGTTTTGTTGCTCATTCCCAGCCAGCCGCCATTCATAAAACg                                                                                       >  2:89714/1‑138 (MQ=255)                                      ccATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCGTTAGCTTTGAGCTTTTATTACTCATTAAAACCCAACCCCCAATAaaaaaaaa                                                                                       >  6:297972/1‑136 (MQ=255)                                                 tGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTCTTCTTTTGTTGCTCATTAACAGCCATCCCCCAATCATAAAAACAAACCAGTcc                                                                             >  8:16525/1‑137 (MQ=255)                                                  gCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  4:108210/138‑2 (MQ=255)                                                   cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGGTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                         >  1:35259/1‑139 (MQ=255)                                                     ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                       <  6:76013/139‑1 (MQ=255)                                                           ttctctcTTGCTTTGGGGATTCTTTGATGTTTTGC‑AAAAATCTTCTTTTATAATGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATCCACTGCCATCCGCCTTTCATAAAACGAAACGAGTCCTGCGGCCTGCAc                                                                 <  4:92805/138‑1 (MQ=255)                                                                tcCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagccg                                                            >  5:373973/1‑139 (MQ=255)                                                                  cTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCATCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                          <  8:89524/139‑2 (MQ=255)                                                                  cTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTATTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACagccgcc                                                          <  1:129698/139‑2 (MQ=255)                                                                  cTGCTTTGGAGATTCCTTTAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGATTTCTGCTTTTCTTCCTCCTTCACAGCCATCAGCCATTCATAAAACGAAACGATTCCTGCGGCCTGCACagccgcc                                                          <  1:55722/139‑2 (MQ=255)                                                                        tGGCGATTCCTTGACGGCTTGC‑AAACATCTTCTTTTAAAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCCCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                     <  1:161658/138‑1 (MQ=255)                                                                        tGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                     <  5:120922/138‑1 (MQ=255)                                                                         ggggATTCCTTGAAGGCTTGC‑AAACATTTTCTTTTAAAACGCGCGGCGGCTTTCGCTTTCTGGTTTTCTTCCTCATTCAAAGCCATCCGCAATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  2:122894/136‑1 (MQ=255)                                                                         ggCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  3:32726/139‑1 (MQ=255)                                                                         ggCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                     <  5:276299/137‑1 (MQ=255)                                                                          gCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  6:373973/138‑1 (MQ=255)                                                                          gCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTAACAGCCCTCCGCCATTAATAAAACGAAACGAGTACAGCGGCCTGCACACACGCCGTTCTTTg                                                  >  4:175954/1‑139 (MQ=255)                                                                          agacagccTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTg                                                  <  3:132888/133‑1 (MQ=255)                                                                               tCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCTGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTAt                                              <  7:91940/138‑1 (MQ=255)                                                                                    gacggCTTGC‑AAACTGATCATTTTTTCACGCGGGGCGGCTTTCGCTTTCTGCTTTTCTACCTCAATCACAGCCACCCGCCATACATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                         <  2:391580/135‑1 (MQ=255)                                                                                    gAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                         <  4:112540/138‑1 (MQ=255)                                                                                    gAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                         <  1:337114/138‑1 (MQ=255)                                                                                    gAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                         <  3:91949/138‑1 (MQ=255)                                                                                     aaGGCTTGC‑AACCATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  2:40568/139‑1 (MQ=255)                                                                                      aGGCTTGC‑AAACATCTTCTTTTATAACGCGCGGCGGGTTTCGCTTTCTGCTTTTCTTCCTCATTCACTGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  1:329508/138‑1 (MQ=255)                                                                                        gCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  6:304222/139‑1 (MQ=255)                                                                                        gCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  2:321461/139‑1 (MQ=255)                                                                                        gCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  7:257657/139‑1 (MQ=255)                                                                                                             ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc                 <  1:89714/138‑1 (MQ=255)                                                                                                             ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc                 <  3:397709/138‑1 (MQ=255)                                                                                                                      ggcggcTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt       <  6:160304/139‑1 (MQ=255)                                                                                                                           cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAACCGCCAATCATAAACCGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc  <  5:212173/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                              ATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724846‑3725105

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