Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	3,724,977	T→C	Q64Q (CAA→CAG)	U712_19205 ←	Uncharacterized protein ywcI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,724,977	0	T	C	97.8%	131.3 / NA	46	Q64Q (CAA→CAG)	U712_19205	Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/36);  minor base A (1/0);  total (10/36)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724843‑3725105                                                                                                                                        |                                                                                                                                 ccAATTTGCTTTTGACTGGGATTTTACATCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTTGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTgct                                                                                                                              >  8:91940/1‑137 (MQ=255)    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                                         <  3:91949/138‑1 (MQ=255)                                                                                       gAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                         <  1:337114/138‑1 (MQ=255)                                                                                        aaGGCTTGC‑AACCATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  2:40568/139‑1 (MQ=255)                                                                                         aGGCTTGC‑AAACATCTTCTTTTATAACGCGCGGCGGGTTTCGCTTTCTGCTTTTCTTCCTCATTCACTGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                       <  1:329508/138‑1 (MQ=255)                                                                                           gCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  7:257657/139‑1 (MQ=255)                                                                                           gCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  6:304222/139‑1 (MQ=255)                                                                                           gCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  2:321461/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)                                                                                                                                        |                                                                                                                                 CCAATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGC‑AAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724843‑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: ‑