Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,445,063	T→C	10.9%	intergenic (‑438/‑221)	yveA ← / → PP_1260	aspartate‑proton symporter/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,445,063	0	T	C	10.9%	172.6 / 7.7	64	intergenic (‑438/‑221)	yveA/PP_1260	aspartate‑proton symporter/transposase
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (2/5);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGG  >  NC_002947/1444915‑1445091                                                                                                                                                      |                             gTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTa                             >  2:5274677/1‑150 (MQ=255) gTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTa                             >  2:6045608/1‑150 (MQ=255) gTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTa                             >  2:7288556/1‑150 (MQ=255)  tACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCGTCAGTTCCGCAGCGACGCGGACATTTTGAAAAAGAAAAGGCGGTTTGCCCGGTAg                            >  1:1231469/1‑150 (MQ=255)   aCGGCGCAGAATTCGTAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAgg                           <  1:6677224/150‑1 (MQ=255)   aCGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAgg                           <  2:1232846/150‑1 (MQ=255)   aCGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAgg                           <  1:4218632/150‑1 (MQ=255)     ggCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACACTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGc                         <  1:2756439/150‑1 (MQ=255)      gcgcAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGcc                        >  1:7181315/1‑150 (MQ=255)       cgcAGAATTCGCCATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCa                       <  2:6693669/150‑1 (MQ=255)       cgcAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCa                       >  1:7444026/1‑150 (MQ=255)       cgcAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGAAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCa                       <  2:1996671/150‑1 (MQ=255)         cAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgc                     <  1:231519/150‑1 (MQ=255)         cAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGGAGGGCCAgc                     >  1:6312822/1‑150 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGCGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    <  2:4407573/150‑1 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAATGGCGGTTGGCCCGGTAGGGCCAgcc                    >  1:2397742/1‑150 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    <  1:2624010/150‑1 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    >  2:295557/1‑150 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    >  2:295491/1‑150 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    >  2:1440631/1‑150 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    <  2:6111597/150‑1 (MQ=255)          aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    <  1:5067093/150‑1 (MQ=255)           gAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCAGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAcccc                   >  1:2797613/1‑149 (MQ=255)            aaTTCGACATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGTGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAaccgc                  <  2:6502249/150‑1 (MQ=14)             aTTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCTACCGGGCCAaccgcc                 <  2:4501473/150‑1 (MQ=14)             aTTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAaccgcc                 <  2:3293665/150‑1 (MQ=14)              ttCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCt                >  1:4855799/1‑150 (MQ=255)              ttCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACAAGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCt                >  1:3502389/1‑150 (MQ=255)                cGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAaccgcc                 <  1:5136381/147‑1 (MQ=14)                cGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAaccgcc                 >  2:5136381/1‑147 (MQ=14)                 gAAATCTGTGGGGGCGCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGTCATTTTGTAATAGAAAAGGGGGTTGGCCCGGTAGGGCCAGCCGGCtttt             >  2:3257416/1‑150 (MQ=14)                 gAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGGTGTCATTTTGTAATAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCtttt             >  2:6866959/1‑150 (MQ=255)                 gAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACAGGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCtgtt             >  2:1164328/1‑150 (MQ=255)                 gAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCACCAAACATGAGGTGTACGGCTTCAGTTCCGCACCGACGCTGACATATTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCtttt             >  1:5167888/1‑150 (MQ=255)                  aaaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCCAGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCtttt             <  2:4792210/150‑1 (MQ=255)                  aaaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTg            >  1:1769537/1‑150 (MQ=39)                  aaaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTg            <  1:4121244/150‑1 (MQ=39)                  aaaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACCCTGACAGTTTGAAAAATAAAAGGCGGTTGGCCCGGTAGGGCCAGCCCCCGTTTg            >  1:3832506/1‑150 (MQ=14)                   aaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGc           <  2:7335848/150‑1 (MQ=38)                   aaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGc           <  2:216378/150‑1 (MQ=38)                   aaTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCGTTTGc           >  1:254809/1‑150 (MQ=38)                    aTCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCt          <  1:4273690/150‑1 (MQ=38)                     tCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTAAGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAATAAAAGGCGGTTGGCCCGGTAGGGCCACCCGCCTTTTGCtt         >  1:3783777/1‑150 (MQ=255)                      cTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCCTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCttt        <  2:5760495/150‑1 (MQ=25)                       tgtgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCCGCCGCCTTGTGCTTTg       >  2:785261/1‑150 (MQ=25)                       tgtgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTg       <  2:3267908/150‑1 (MQ=35)                       tgtgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGGACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTg       <  1:2329764/150‑1 (MQ=25)                        gtgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGt      <  2:5425114/150‑1 (MQ=34)                        gtgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACAAGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTTAAAAATAAAAGTCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGt      >  1:2506093/1‑150 (MQ=255)                        gtgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACAAGAGGTGTACGGCGTCAGTTCCGCAGCGACCCGGACATTTTGAAAAAGAAAAGGCGGTTGGCCCCGTAGGGCCAGCCGCCGTTTGCGTTGt      >  1:1591600/1‑150 (MQ=11)                         tgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGtt     <  2:3030066/150‑1 (MQ=33)                         tgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGtt     <  1:755784/150‑1 (MQ=33)                         tgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGtt     >  1:6019386/1‑150 (MQ=33)                         tgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGtt     <  2:7181315/150‑1 (MQ=33)                         tgGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACAGTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGGTTTGtt     >  1:1875305/1‑150 (MQ=18)                          gggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTa    <  2:1498504/150‑1 (MQ=32)                          gggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTa    <  1:3556965/150‑1 (MQ=32)                          gggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTa    <  1:1886110/150‑1 (MQ=32)                           ggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGTAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCGTTTGCTTTGTTAg   >  1:5668890/1‑150 (MQ=17)                           ggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAATAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAg   >  1:5668603/1‑150 (MQ=17)                           ggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGGTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAg   >  2:2477120/1‑150 (MQ=17)                           ggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCtt               <  2:2402750/138‑1 (MQ=11)                           ggggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCtt               >  1:2402750/1‑138 (MQ=11)                            gggCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAgg  <  2:7336433/150‑1 (MQ=31)                                                                                                                                                      |                             GTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGG  >  NC_002947/1444915‑1445091

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑