Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,186,053	T→G	7.9%	E136A (GAG→GCG) ‡	aceK ←	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,186,053	0	T	G	7.9%	214.1 / 9.8	76	E136A (GAG→GCG) ‡	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base G (3/3);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

GCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTGGATATAGGC  >  NC_002947/5185908‑5186201                                                                                                                                                   |                                                                                                                                                      gcaTTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTcgc                                                                                                                                                    >  7:141611/1‑149 (MQ=255)       gcgcAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCTCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCcgc                                                                                                                                             <  8:287903/149‑1 (MQ=255)       gcgcAGTTGTTCTTCCAGCCGCGCCAGGGCGCCATCCAGGTCTCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCcgc                                                                                                                                             <  8:282754/149‑1 (MQ=255)          cAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCa                                                                                                                                          <  6:191669/149‑1 (MQ=255)           aGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCAt                                                                                                                                         >  8:298830/1‑149 (MQ=255)               gttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                           >  5:162982/1‑143 (MQ=255)               gttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                           <  6:162982/143‑1 (MQ=255)                ttcttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGaa                                                                                                                                     <  5:313440/148‑1 (MQ=255)                   ttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGcc                                                                                                                                 >  2:161214/1‑149 (MQ=255)                   ttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGcc                                                                                                                                 >  4:110112/1‑149 (MQ=255)                     ccAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCg                                                                                                                                <  1:161214/148‑1 (MQ=255)                      cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCggg                                                                                                                              >  6:386021/1‑149 (MQ=255)                      cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCgg                                                                                                                               <  5:19081/148‑1 (MQ=255)                        gCCGCGCCAGGTCGCCTTCCCGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                           <  1:23783/134‑1 (MQ=255)                        gCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                           >  2:23783/1‑134 (MQ=255)                        gCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTa                                                                                                                            <  3:110112/149‑1 (MQ=255)                         ccGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAg                                                                                                                           <  5:84576/149‑1 (MQ=255)                               cAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGa                                                                                                                     <  2:282819/149‑1 (MQ=255)                                ctgcagCCCTCCCGCTCGCCCTAGCGCATGTCCAACGGGGACTCGGCGCACCCCGCGCGCAGCACGTTCCTGCCACCGCCGTCCAGCCGATCGGTCCGGGTGTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAACGCCGGGTAGTGCGGAt                                                                                                                    <  8:348848/144‑1 (MQ=255)                                atgTCGCCGTCCTGGTCGCCATAGGGCATGTCGAACGGGTAATCTGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGAt                                                                                                                    <  4:55078/147‑1 (MQ=255)                                  gTCGCCTTCCAGGTAGCCATAGGGCATGTCGTACGGGTAATCGGCGAACAGCGCGCGCAGCAAGTTCCTGAGACTGCCGTCTAGCCGATAGGTACGGGTTTGCGGGG‑CGCGCTCGCG‑GCCGCGGATGGAAGGCCGGGTAGTGTGGATGa                                                                                                                  <  4:145931/149‑1 (MQ=255)                                   tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                                   <  6:231420/131‑1 (MQ=255)                                   tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                                   >  5:231420/1‑131 (MQ=255)                                   tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTATTGTGGAt                                                                                                                    >  7:348848/1‑146 (MQ=255)                                     gCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCgg                                                                                                                               <  3:372437/133‑1 (MQ=255)                                     gCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCgg                                                                                                                               >  4:372437/1‑133 (MQ=255)                                     gCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAAc                                                                                                                <  3:145566/148‑1 (MQ=255)                                            aGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑C‑CGCGAGCGCGCCGCGCaaaccc                                                                                                                                     <  8:186985/120‑6 (MQ=38)                                            aGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑C‑CGCGAGCGCGCCGCGCaaaccc                                                                                                                                     >  7:186985/1‑115 (MQ=38)                                                          cATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAgg                                                                                           >  7:190729/1‑148 (MQ=255)                                                           aagTCGAACGTGTAATCGTCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCTTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGt                                                                                          <  8:141611/146‑1 (MQ=255)                                                                gAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                                     <  7:298830/148‑1 (MQ=255)                                                                gAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGc                                                                                    >  7:347792/1‑149 (MQ=255)                                                                  aCGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTg                                                                                  <  8:190729/149‑1 (MQ=255)                                                                        aaTCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCACATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATc                                                                                              <  7:281352/131‑1 (MQ=255)                                                                        aaTCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCACATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATc                                                                                              >  8:281352/1‑131 (MQ=255)                                                                          tCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTa                                                                                                                            >  4:31133/1‑99 (MQ=255)                                                                          tCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCt                                                                          >  4:95285/1‑149 (MQ=255)                                                                           cGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                                   >  3:83963/1‑91 (MQ=255)                                                                           cGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                                   <  4:83963/91‑1 (MQ=255)                                                                             gCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGc                                                                        >  1:8328/1‑148 (MQ=255)                                                                             gCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCAGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCTTCGCTGATCAGGTCATGGCTGAAGAGGCTGc                                                                        >  8:173549/1‑148 (MQ=255)                                                                               gAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaa                                                                     <  5:386021/149‑1 (MQ=255)                                                                                aaCACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaaa                                                                    >  1:15265/1‑149 (MQ=255)                                                                                 acacCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaaa                                                                    >  7:108891/1‑148 (MQ=255)                                                                                  cacCGCGCGCAGCAAGTTCCTGAGACAGCCGTCAAGCCGATAGCTACGGGTTTGCGCGG‑CGCGCTCGAG‑GCCGCGCATCCAAGGCCCGGTAGAGTGGAAGAACAAGCAGCCGTCGCTGAACAGGTCAAGGCTGAAGACGCTGCAAAACa                                                                  >  6:69482/1‑149 (MQ=255)                                                                                   acCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                                     <  6:258069/129‑1 (MQ=255)                                                                                   acCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAg                                                                 >  4:211704/1‑149 (MQ=255)                                                                                     cgcgcgCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGAAAAACAGc                                                                >  3:91391/1‑148 (MQ=255)                                                                                              caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAcc                                                      <  6:388913/149‑1 (MQ=255)                                                                                              caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                       >  4:31716/1‑148 (MQ=255)                                                                                              caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                       >  2:56595/1‑148 (MQ=255)                                                                                               aaGTTCCTGAGACTGCCGTCAAGCCGATAGGTTCGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                                     <  5:69482/149‑1 (MQ=255)                                                                                               aaGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                                     <  2:129919/149‑1 (MQ=255)                                                                                                aGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAg                                                    <  5:127973/149‑1 (MQ=255)                                                                                                  ttCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGt                                                            <  4:22105/139‑1 (MQ=255)                                                                                                  ttCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGt                                                            >  3:22105/1‑139 (MQ=255)                                                                                                      tGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTctct                                              >  1:249755/1‑149 (MQ=255)                                                                                                        cgaCTGCCGTCATGGCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GGCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGc                                            <  4:56008/148‑1 (MQ=255)                                                                                                        agaCTTCCGTCAAGCCTAGAGTTACGTGTTTGCGCGG‑CGAGCTCGCG‑TCCGCTCATGGAAGGCCGGGTAGTGTGTATTAACATGCAGCCGTCGCTGATCAGGTCATGTCTGAAGAGGCTGCAAAACAGCTAGTTGTACCAGGTctct                                              <  4:91391/147‑1 (MQ=255)                                                                                                        agaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTctct                                              <  3:95285/147‑1 (MQ=255)                                                                                                         gaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGcc                                           <  3:211704/149‑1 (MQ=255)                                                                                                          aCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGcc                                           >  5:257426/1‑148 (MQ=255)                                                                                                              ccGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtc                                      <  2:8328/149‑1 (MQ=255)                                                                                                              ccGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtc                                      <  8:347792/149‑1 (MQ=255)                                                                                                               cGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcg                                     <  2:15265/149‑1 (MQ=255)                                                                                                                      ccGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCa                              <  2:249755/149‑1 (MQ=255)                                                                                                                        gATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAGACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAgg                            >  5:342879/1‑149 (MQ=255)                                                                                                                               aCGCGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGACGGGTAGTGTGGAAGAACATGCAGCCGTCGCTGATCATGTCATGGATGAAGAGGCTACAAAACAGCGAATTGTACGAGGTATCTGACAACTCGACGTCCAGGCGGGGGt                     >  2:268989/1‑149 (MQ=255)                                                                                                                                cGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGt                     <  3:31716/148‑1 (MQ=255)                                                                                                                                    tttGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAATACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATc                <  6:342879/149‑1 (MQ=255)                                                                                                                                     atGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtt                                                           <  2:135566/104‑1 (MQ=255)                                                                                                                                     atGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtt                                                           >  1:135566/2‑105 (MQ=255)                                                                                                                                      tGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg                       <  3:193199/140‑1 (MQ=255)                                                                                                                                      tGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg                       >  4:193199/1‑140 (MQ=255)                                                                                                                                         gcgg‑cgcgCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTg           <  6:257426/149‑1 (MQ=255)                                                                                                                                                   tCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTGGATATAg    >  2:254876/1‑147 (MQ=255)                                                                                                                                                   tCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCGGAAGAGGCTGCAAACCCGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGAGCAGGTGGATATAGGc  >  2:373473/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTGGATATAGGC  >  NC_002947/5185908‑5186201

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