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	94.4%	42.2 / NA	18	A68G (GCC→GGC)	LY180_20645	rhsA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (10/7);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGC  >  CP006584/4304906‑4305307                                                                                                                                                                                                                                                                                                        |                                                                                                           tCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGcccggccc                                                                                                         <  1:311234/299‑1 (MQ=255)  cTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggcccg                                                                                                        <  1:453546/299‑1 (MQ=255)     aaaaaTTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGCTGCAAAGGTCCTTCCCGGTGACACCGACATCGCCCTGCCCGGCCCGCTg                                                                                                     >  1:425874/1‑299 (MQ=255)        aaTTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCCTTGTCCACGGTTCAGCCGCGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGCTCGGGGTGCCCCGGCGCAGTTACGTCCGCCCCTCCGGCCAACCCCCCGCTCTGTGCAAAGCTCCTTCCCCGTGAAACCCACACCCCCCCGCCCGCCCCCCCGccc                                                                                                  >  2:557454/1‑298 (MQ=21)             gggggTGTTAACTGTTTATCATTTCCTGTGTTTGACTGGCGTAAGGCTTCTGGTTGTTATTTATTTCAAAAGGATTTTAGGGTGGGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAt                                                                                             <  2:122617/298‑1 (MQ=17)                 gtgtTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGCCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGCCCCTCCCGGTGAAACCGACACCGCCCcccccgcccccccgcccgtcctccct                                                                                         >  2:408811/1‑289 (MQ=18)                  tgtTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGATGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct                                                                                        >  1:96533/1‑299 (MQ=255)                     tAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTccc                                                                                     <  1:239361/299‑1 (MQ=255)                          cATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCAcc                                                                                <  1:390973/299‑1 (MQ=255)                                    tCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGtt                                                                      >  1:132769/1‑299 (MQ=255)                                        gttgttGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTAATACTCTCCCGCACCTAAAGCAGTTAcag                                                                  >  1:429005/1‑297 (MQ=255)                                          tgttgACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTTAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGa                                                                >  1:144917/1‑299 (MQ=255)                                              gagtGGCGCTGGGCTTCGGGTTCTTATTTATTTCGAAAGGATTTGAGGAGGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACaaa                                                            <  2:132769/296‑1 (MQ=14)                                                     gcTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCACAGGTCCTTCCCGTTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACCGCAGTTACCGGACAAAAACGcac                                                     >  1:254719/1‑297 (MQ=255)                                                                               cAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCgg                           >  1:213595/1‑299 (MQ=25)                                                                                   aGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTgg                       <  1:171529/299‑1 (MQ=21)                                                                                        tttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAAt                  >  1:573241/1‑299 (MQ=17)                                                                                             gATGGGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTg             <  1:74124/299‑1 (MQ=14)                                                                                                   cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGata       >  1:510524/1‑299 (MQ=14)                                                                                                      aaaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGTGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc    <  1:123058/299‑1 (MQ=14)                                                                                                        aaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGc  >  1:172593/1‑299 (MQ=14)                                                                                                                                                                                                                                                                                                        |                                                                                                           TCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGC  >  CP006584/4304906‑4305307

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