Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	370,869	0	T	G	42.6%	104.9 / 35.2	68	L198W (TTG→TGG)	mhpB	2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (36/3);  new base G (0/29);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 3.61e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGC  >  NC_000913/370719‑371015                                                                                                                                                              |                                                                                                                                                   tCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGtt                                                                                                                                                    >  2:242825/1‑151 (MQ=255)  ccAGCGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTg                                                                                                                                                   <  1:7456/151‑1 (MQ=255)    aGCGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTggg                                                                                                                                                 >  1:1334409/1‑151 (MQ=255)      cGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTgggg                                                                                                                                                >  1:1056594/1‑150 (MQ=255)      cGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTgggg                                                                                                                                                >  1:872445/1‑150 (MQ=255)        tACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGc                                                                                                                                             >  2:519012/1‑151 (MQ=255)            cGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCggg                                                                                                                                          >  1:527084/1‑150 (MQ=255)                tGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGaa                                                                                                                                        >  2:868593/1‑148 (MQ=255)                   tGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGAtt                                                                                                                                   <  2:127024/150‑1 (MQ=255)                   tGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGAtt                                                                                                                                   <  1:954702/150‑1 (MQ=255)                                        tCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                             >  2:1436416/1‑151 (MQ=255)                                        tCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                             >  2:741644/1‑151 (MQ=255)                                         caccacCCCTCAAAAAAACCGGGTGCTGTTCC‑GGGTTTCGGGGGGTTTTCCCCAT‑CGGCCGCGGGG‑CCCC‑GAG‑CGGGGAAAA‑GCCGATCCCCAT‑ATGCGCGCCCGTTTGTGGGGGAGGGGGAAATATTTCCCCGCCAGTGAGCGCGAATTg                                                                                                            <  1:1344212/151‑1 (MQ=255)                                          aCCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGc                                                                                                           >  2:346977/1‑151 (MQ=255)                                          aCCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGc                                                                                                           >  1:966721/1‑151 (MQ=255)                                           ccacccctctAAATAACCGGGTGTTTTTCCTGGTTTCGGGTGGGTTTT‑CCCAT‑CGCCGGCCGGT‑GCCC‑GAA‑CGGGAGAAA‑GCCGATCCCAAT‑CTGCCCGACCTTTTTTGGGGGGGGGGAAAAGATTTCCCCGCCAGTGAGCGCGAATTg                                                                                                            <  1:1476136/143‑1 (MQ=255)                                           ccAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCg                                                                                                          >  1:1402311/1‑151 (MQ=255)                                             aGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTc                                                                                                        >  1:540840/1‑151 (MQ=255)                                                   caaaacaaccgcgTGTTGTGGCGGGGTTCCGG‑GGGCTTT‑CCCCTGCGGCGGCGGGC‑GCCC‑GAA‑GGGGAAAGG‑GCGGAGCCCCAT‑AGGCGCGCCCTTTGGTGGGGGGGGGGGAAATATTGTCCCGCCAGTGAGCGCGAATTGCGTCAGCaa                                                                                                   <  2:1202074/141‑1 (MQ=255)                                                    tcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACg                                                                                                 >  2:1017974/1‑151 (MQ=255)                                                    tcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACg                                                                                                 >  1:119453/1‑151 (MQ=255)                                                      aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                               >  1:859743/1‑151 (MQ=255)                                                         aaaCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTg                                                                                             >  2:1390003/1‑150 (MQ=255)                                                           aCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                          >  2:1224844/1‑151 (MQ=255)                                                            cctctgcgtttCCTGGTTTGCGGGGGTCTTC‑CCCTT‑CCGCCGCC‑GT‑TCCC‑GAACCTGGAG‑AA‑GCCCCGGCCCATAACGCGCGCCCGTTTGTGGGGGGGGGGAAAATATTTTCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAt                                                                                           <  1:541900/140‑1 (MQ=255)                                                                  ttGTTCCTGGGTCCGGGTGGGTTTC‑CCAAC‑CAGCCGCGGGG‑GCCC‑GAA‑CGGGCGAAA‑GCCGACGCCCAT‑AGGCGCGCCCGTTTGTGGGGGGGCGGAAAAATTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                          <  1:310194/143‑1 (MQ=255)                                                                       ccTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                >  1:881945/1‑149 (MQ=255)                                                                        cTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                >  1:25197/1‑148 (MQ=255)                                                                        cTGGGCTCCGGGTGGCCTT‑CCCCA‑CCGCCCCGGGAGGCCC‑GAG‑CTGGAGAAA‑GCTGATGCCCTT‑GGGCCCGACCGTTTGTGGGGGGGGGAAAAATTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTg                                                                                <  2:944735/149‑1 (MQ=255)                                                                          gggTCCGGGGGGGTTTC‑CCCAC‑CGGCCGCCGGG‑GCCG‑GAC‑GTGGGAAAG‑GCGGGGGCCGAT‑CGGCGCGCCGGGTTGTGGGGGGGGGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTgag                                                                              <  2:173811/148‑1 (MQ=255)                                                                           ggTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGt                                                                          >  2:717683/1‑151 (MQ=255)                                                                              tCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttg                                                                       >  2:324665/1‑151 (MQ=255)                                                                                cGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                     >  1:1301322/1‑151 (MQ=255)                                                                                     ggCTTC‑CCCAG‑CGGCCGCGGGT‑CCCC‑AAC‑CTGGCAAAA‑GCCGATGCCCTT‑GTGCGCGACCGTGTTTGGGGGGGGGGAAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGa                                                                <  2:632778/151‑1 (MQ=255)                                                                                        ttt‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                              >  1:193590/1‑150 (MQ=255)                                                                                        ttt‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                              >  1:881797/1‑150 (MQ=255)                                                                                           ccccAT‑CGCCCGCGGGT‑CCCA‑AAA‑GTGGAAAAC‑GCCGATGCCAAT‑TTGCGCGCCCGTTTGTGGGGGGGGGGAAAATATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATc                                                              <  1:1217422/147‑1 (MQ=255)                                                                                            cccAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                            >  1:1211945/1‑149 (MQ=255)                                                                                             ccAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                            >  2:866915/1‑148 (MQ=255)                                                                                              cccaccgcccgcgGG‑GCCCGGAA‑CGGGCAAAA‑GCCGATCCCAAT‑ATGCGCGCCCGTTTGTGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCaga                                                           <  1:1139037/139‑1 (MQ=255)                                                                                                 ccAGCCGCCGGT‑GCCC‑AAC‑CTGGCAAAA‑CCCGACCCCCAT‑ATGCCCGACCGTTTGTGGGGGAGGGGAAAAATTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAgaga                                                         <  1:1055662/147‑1 (MQ=255)                                                                                                     ccgccgGT‑GCCC‑GAA‑CTGGAGAAA‑GCCGATGCCAAT‑ATGCGCGACCGTTTGTGGGGGGGGGGGAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCt                                                    <  1:487623/149‑1 (MQ=255)                                                                                                     ccgccgGG‑GCCA‑G‑A‑CGGGGGAAAGGCCGATCCCCAT‑ATGCGCGCCCGTCTGTGGGGGGGGGGAAAAGATTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGc                                                  <  1:129236/151‑1 (MQ=255)                                                                                                         cgGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATcc                                              >  2:978431/1‑151 (MQ=255)                                                                                                                 c‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAc                                       >  1:832206/1‑151 (MQ=255)                                                                                                                    aa‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAccc                                     >  2:311932/1‑151 (MQ=255)                                                                                                                       ggggagAAA‑GCCGACCCCCAT‑GGGCGCGACGTTTTGTGGGGGAGCGGGAAAGATTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGa                                   <  2:369900/146‑1 (MQ=255)                                                                                                                         ggaaaaa‑GCCGATCCCCAT‑ATGCGCGCCCGTTTGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAt                                  <  1:868609/146‑1 (MQ=255)                                                                                                                             aaa‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTggg                             >  1:261316/1‑151 (MQ=255)                                                                                                                             aaa‑GCCGATCCCCAT‑GGGCGCGCCCGTCTGTGGGGGGGCGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAttt                                <  1:1386329/148‑1 (MQ=255)                                                                                                                              aaaccggtgccccAT‑AGGCGCCCCCGTTTGTGGGGGAGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTg                               <  1:358987/139‑1 (MQ=255)                                                                                                                               acgcggATGCCCAT‑ATGGGGGACGGTTTGTGGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTg                               <  1:710051/143‑1 (MQ=255)                                                                                                                                  ccgttccctat‑agGCGCGACCGTTGGTGGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATaa                         <  2:980281/142‑1 (MQ=255)                                                                                                                                  ccgacccccAT‑ACGCGCGCCGGTTTGGGGGGGAGGGGAAAATTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATaa                         <  2:649518/145‑1 (MQ=255)                                                                                                                                  ccGTTGCCGAT‑ATGCGCGCCGTTCTGTGGGGGTGCGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTg                               <  1:205912/145‑1 (MQ=255)                                                                                                                                  ccGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATa                          >  2:91503/1‑150 (MQ=255)                                                                                                                                   ccgttccccagttGCGCGCCCGTTTGGGGGGGAGGGAAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATa                          <  2:643726/138‑1 (MQ=255)                                                                                                                                    gATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAAcc                       >  2:719699/1‑151 (MQ=255)                                                                                                                                      atcccaatgtGCGCGACCGTTTGTGGGGGGGGGGAAAATTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCa                      <  1:467762/142‑1 (MQ=255)                                                                                                                                        cccAT‑AGGCGCGCCCGTTTGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCa                      <  1:1032094/148‑1 (MQ=255)                                                                                                                                           at‑atGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                >  1:1120703/1‑151 (MQ=255)                                                                                                                                           at‑agGCGCGCCCGTCTGTGGGGGGGGGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                <  1:1143929/151‑1 (MQ=255)                                                                                                                                                gcgcgcCGGTGTGTGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCa                 <  2:1489677/146‑1 (MQ=255)                                                                                                                                                 cggacccGTGTGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACtt           <  1:596370/146‑1 (MQ=255)                                                                                                                                                   cgcccGTTTGTGGGGGGGGGGAAAATTTTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATg               <  1:810250/142‑1 (MQ=255)                                                                                                                                                         tCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTgg     >  1:1212963/1‑149 (MQ=255)                                                                                                                                                         tCTGTGGGGGAGGGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGc        <  1:1305664/146‑1 (MQ=255)                                                                                                                                                          cTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc  >  1:458743/1‑151 (MQ=255)                                                                                                                                                           tGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc  >  2:379041/1‑150 (MQ=255)                                                                                                                                                              |                                                                                                                                                   TCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTT‑CCCAT‑CAGCCGCCGGT‑GCCC‑GAA‑CTGGCGAAA‑GCCGATGCCCAT‑ATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGC  >  NC_000913/370719‑371015

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