Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,305,201	C→G	A68G (GCC→GGC)	LY180_20645 →	rhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,305,201	0	C	G	95.8%	56.5 / NA	24	A68G (GCC→GGC)	LY180_20645	rhsA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (8/15);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTG  >  CP006584/4304910‑4305341                                                                                                                                                                                                                                                                                                    |                                                                                                                                             aaaagTTCCGGGGTGTTAACTTATTTTCATTTCCTTTTGTTGACTGGCGCTCGGCTTCTGGTTCTTATTTATTTCAAAGGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCt                                                                                                                                        <  2:123544/298‑1 (MQ=255)    aaTTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCCGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGCTCGGTGTGCCCCGGCGGCGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTCCCCGCCCCGCTGcc                                                                                                                                     >  2:484950/1‑298 (MQ=255)                       aTTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCCCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCCCCCTGCCCGGCCCGCTGCCGTTCATCCTCCCCCGcccc                                                                                                                  >  2:250368/1‑298 (MQ=255)                          aTCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCCGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTAc                                                                                                               >  1:322677/1‑298 (MQ=255)                               ttCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTAcagcag                                                                                                          <  1:415872/298‑1 (MQ=255)                                tCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGt                                                                                                         >  1:232803/1‑298 (MQ=255)                                    gttgttGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTAcc                                                                                                     <  1:157350/298‑1 (MQ=255)                                      tgttgACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCgg                                                                                                   <  1:165474/298‑1 (MQ=255)                                                  cTAGGCTGTGGGTTCTTATTTATTTAAAAAGGATTTTACGATGAGCGGAAACCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGccc                                                                                       <  1:354251/298‑1 (MQ=255)                                                   tAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcg                                                                                      <  1:113420/298‑1 (MQ=255)                                                   tAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCCCTGCCGTTCATCCTCTCCCCGACCTACCGCCGTTACCGGCCCAAAACGcccc                                                                                      >  2:306063/1‑297 (MQ=14)                                                            ggTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTgg                                                                             <  1:110184/298‑1 (MQ=255)                                                                tttttttttttCAAAAGGTTTTGACGATGGGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAg                                                                         <  2:254643/294‑1 (MQ=255)                                                                   atttttttAAAAGGGTTTTGACGGGGAGGGGAAAGCCGGGGGCGCGGCAGGGGGAGAGGCGGCAGTATGGGGGTAGCATTGGCCCGGGTTCAGCGGGGGTACGTGTTGGTGCCCCCACCGGGGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCt                                                                      <  2:412485/298‑1 (MQ=11)                                                                   atttatttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCATGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTCGGGAGccc                                                                      >  1:18191/1‑297 (MQ=255)                                                                              aaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGACGGGGGGGAGCCTCGGCCCCGGca                                                           >  1:414896/1‑297 (MQ=255)                                                                                          aTGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCCCCCGCCCCTACCGCAGTTCCCGGACCAAACCCACCCCCCCGGGGGGGGCCCCCGGCCCCGGCCGGAAAATGcccc                                               >  2:458701/1‑296 (MQ=11)                                                                                           tGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGTGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGc                                              >  1:58768/1‑298 (MQ=14)                                                                                              gCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGa                                           <  1:565261/298‑1 (MQ=14)                                                                                                ggAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGata                                         <  1:250368/298‑1 (MQ=14)                                                                                                   aagccGGCGGCGCGGCGGGGCGACATGACACAGGATGTCGGGGGCATTGGCCAGGGTTCAGCCGGGGTACGTATTGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCGGCCAGCCGGTCAAGCCCCTGCTCGGTGCAACGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc                                      <  2:104807/295‑1 (MQ=11)                                                                                                   aaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc                                      <  1:165777/298‑1 (MQ=14)                                                                                                          cggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa                               <  1:383063/298‑1 (MQ=11)                                                                                                                                       ggCGGTAGCCTTGTCCGGGGTGCAGCCGGGGTCCGTATTGGTGCCCCCGCCGGCGTGGCGTGTTCGGGGTGCCCCGACGGGGTGACGTCCGGCGATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTg  <  2:143491/298‑1 (MQ=11)                                                                                                                                                                                                                                                                                                    |                                                                                                                                             AAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTG  >  CP006584/4304910‑4305341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑