Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,445,062	G→A	10.9%	intergenic (‑437/‑222)	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,062	0	G	A	10.9%	174.3 / 10.6	64	intergenic (‑437/‑222)	yveA/PP_1260	aspartate‑proton symporter/transposase
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (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

AGGTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGG  >  NC_002947/1444913‑1445091                                                                                                                                                       |                              cggTACGGCGCAGAATTCGAAAGCTGTGCGGGCCCGCGACGAGATCACCCAGCCCTTCC‑GGCCTGCGCTGTCGCGCAGCCAACATGAGGTCTCCCGCTTCAGTGCCGCAGCGACGCTGCCAGAGTGAAAAAGAAAAGGCCGATGGCCCgg                               >  2:5365353/2‑150 (MQ=14) aGGTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCgg                               <  2:1188477/150‑1 (MQ=255) aGGTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTACGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCgg                               >  1:3176434/1‑150 (MQ=255) aGGTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCGTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCgg                               >  1:3176270/1‑150 (MQ=255)   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:6111597/150‑1 (MQ=255)            aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    >  2:1440631/1‑150 (MQ=255)            aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    <  1:5067093/150‑1 (MQ=255)            aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    >  2:295491/1‑150 (MQ=255)            aGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAgcc                    >  2:295557/1‑150 (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                 >  2:5136381/1‑147 (MQ=14)                  cGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAaccgcc                 <  1:5136381/147‑1 (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)                                                                                                                                                       |                              AGGTACGGCGCAGAATTCGAAATCTGTGGGGGCCCTGTAGGAGATCACCCAGCCCTTCC‑GGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGG  >  NC_002947/1444913‑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: ‑