Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,674,290	A→G	7.5%	L89P (CTC→CCC)	PP_3238 ←	transcriptional regulator PyrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,674,290	0	A	G	7.5%	148.7 / 2.5	53	L89P (CTC→CCC)	PP_3238	transcriptional regulator PyrR
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

GAGCACCGCGCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTGAATGCCTAC  >  NC_002947/3674142‑3674438                                                                                                                                                     |                                                                                                                                                       gaGCACCGCGCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtga                                                                                                                                                        <  3:249780/149‑1 (MQ=255)      ccGCGCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtgagttga                                                                                                                                                   <  5:53725/149‑1 (MQ=255)      ccGCGCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtga                                                                                                                                                        <  8:55660/144‑1 (MQ=255)      ccGCGCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtga                                                                                                                                                        >  7:55660/1‑144 (MQ=255)       cgcgCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtgagttgag                                                                                                                                                  >  3:310279/1‑149 (MQ=255)        gcgcTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtgagttgagt                                                                                                                                                 <  4:240212/149‑1 (MQ=255)           cTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATg                                                                                                                                              <  3:201235/149‑1 (MQ=255)              atgaGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATgc                                                                                                                                             <  7:139017/147‑1 (MQ=255)                   gTGGCGCGCCCCCAGTTCGGCGAGGTAGGCGCAGCTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGc                                                                                                                                       >  1:258748/1‑148 (MQ=255)                   gTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGc                                                                                                                                       >  5:105982/1‑148 (MQ=255)                        gcgcACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGTAAGAgg                                                                                                                                 <  4:310279/149‑1 (MQ=255)                                 gTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAag                                                                                                                                    >  4:320029/1‑137 (MQ=255)                                 gTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAag                                                                                                                                    <  3:320029/137‑1 (MQ=255)                                     ggCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCt                                                                                                                     <  2:198446/148‑1 (MQ=255)                                     ggCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTg                                                                                                                    <  8:272507/149‑1 (MQ=255)                                        gAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCa                                                                                                                  >  1:229800/1‑148 (MQ=255)                                             aGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAgcgc                                                                                                           <  6:105982/149‑1 (MQ=255)                                               gcgcAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCAGGTCATCGACCAGCAGCAAGGTGGTTCCGGCAAGTTGGAGCTTGGGGAGTTCCGGGTTTTCTGTTAGTTGAGTATGCGCATGCAAGATTCTCATTTCGTCTGCA‑TAGCGCtt                                                                                                         <  4:171298/149‑1 (MQ=255)                                               gcgcAGGTGCGTAGCAGAGAGTGCCCTTCGAAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCtt                                                                                                         <  6:210077/149‑1 (MQ=255)                                                       gCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtgagtt                                                                                                                                                     >  4:270239/1‑98 (MQ=255)                                                       gCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGtgagtt                                                                                                                                                     <  3:270239/98‑1 (MQ=255)                                                           aGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCgg                                                                                             >  3:184249/1‑149 (MQ=255)                                                            gCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCggg                                                                                            <  2:229800/149‑1 (MQ=255)                                                            gCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCggg                                                                                            <  5:118392/149‑1 (MQ=255)                                                              agagagTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTa                                                                                          <  5:248028/149‑1 (MQ=255)                                                                agagTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTa                                                                                          >  2:337727/1‑147 (MQ=255)                                                                   gTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCg                                                                                     <  3:198944/149‑1 (MQ=255)                                                                    tGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCgg                                                                                    >  3:3741/1‑149 (MQ=255)                                                                       ccTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCaa                                                                                 <  6:135389/149‑1 (MQ=255)                                                                        cTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTagag                                                                                       >  5:285595/1‑142 (MQ=255)                                                                        cTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTACGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTagag                                                                                       <  6:285595/142‑1 (MQ=255)                                                                           cGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGTGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATg‑g                                                                             <  2:258748/148‑1 (MQ=255)                                                                            gTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑gcg                                                                           >  4:346754/1‑149 (MQ=255)                                                                                  cACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCAc                                                                     >  7:97297/1‑149 (MQ=255)                                                                                       cATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGGGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGtt                                                                 >  1:273118/1‑148 (MQ=255)                                                                                       cATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGGGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGtt                                                                 >  1:267697/1‑148 (MQ=255)                                                                                       cATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCt                                                                                                                     >  5:380578/1‑98 (MQ=255)                                                                                       cATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCt                                                                                                                     <  6:380578/98‑1 (MQ=255)                                                                                                 caGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑gcg                                                                           <  1:21659/128‑1 (MQ=255)                                                                                                 caGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTg                                                        >  7:114440/1‑147 (MQ=255)                                                                                                 caGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTg                                                        >  3:239423/1‑147 (MQ=255)                                                                                                 caGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTg                                                        >  2:274805/1‑147 (MQ=255)                                                                                                 caGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTAAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑gcg                                                                           >  2:21659/1‑128 (MQ=255)                                                                                                  aGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAg                                                     >  2:100951/1‑149 (MQ=255)                                                                                                  aGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAg                                                     <  4:184249/149‑1 (MQ=255)                                                                                                    cAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCgg                                                                                             <  6:62125/108‑1 (MQ=255)                                                                                                     aaGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTg                                                  <  3:346754/149‑1 (MQ=255)                                                                                                               cGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGGGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAgg                                        <  2:273118/149‑1 (MQ=255)                                                                                                               cGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGGGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAgg                                        <  2:267697/149‑1 (MQ=255)                                                                                                                 gCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGgcg                                      >  5:322397/1‑149 (MQ=255)                                                                                                                     ggTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTgc                                  >  4:286557/1‑149 (MQ=255)                                                                                                                     ggTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTgc                                  <  3:286557/149‑1 (MQ=255)                                                                                                                               ggCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCgg                          <  4:3741/147‑1 (MQ=255)                                                                                                                               ggCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCt                        >  6:249948/1‑149 (MQ=255)                                                                                                                                gCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATGAGCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCt                        <  1:100951/149‑1 (MQ=255)                                                                                                                                        cGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCAcgcc               >  2:146431/1‑149 (MQ=255)                                                                                                                                         gggTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCAGTAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCAcgcc               >  1:291468/1‑149 (MQ=255)                                                                                                                                             tttCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTg           >  8:197396/1‑148 (MQ=255)                                                                                                                                                 ggTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTGAATGc      >  6:61960/1‑149 (MQ=255)                                                                                                                                                  gTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTGAATGcc     <  6:322397/149‑1 (MQ=255)                                                                                                                                                  gTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTGAATGcc     >  6:375142/1‑149 (MQ=255)                                                                                                                                                     agttgagtATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCACTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTGAATGCCTAc  <  1:274805/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GAGCACCGCGCTGATGAGGTGGCGCGCACCCAGTTCGGCGAGGTAGGCGCAGGTGCGTAGCAGAGAGTGCCCTTCGTAGAGCACGTCATCGACCACCAGCAAGGTGGTTCCGGCAAGGTCGAGCTTGGCGAGTTCCGGGTTTTCGGTGAGTTGAGTATGCGCATGCAAGAGGCTCAGGTCGTCTGCA‑TAGCGCTTGACCTTGATCGGGTAGAGCGGCAATG‑GCGCAGCACCGTTCAGTGCTTGCAGGTGCAGTTGCAGGCGCTGCGCCAGCGGCTCACCACGCCGCTGAATGCCTAC  >  NC_002947/3674142‑3674438

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