Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	371,120	A→T	9.5%	K201* (AAA→TAA) ‡	dgcA →	dimethylglycine dehydrogenase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	371,120	0	A	T	9.5%	176.2 / 8.1	63	K201* (AAA→TAA) ‡	dgcA	dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (4/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAC  >  NC_002947/370973‑371265                                                                                                                                                    |                                                                                                                                                  tACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGa                                                                                                                                                   >  8:306886/1‑148 (MQ=255) tACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGa                                                                                                                                                   >  3:17062/1‑148 (MQ=255)      gcAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCg                                                                                                                                             >  5:205548/1‑149 (MQ=255)        aGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCTTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGc                                                                                                                                           >  8:345938/1‑149 (MQ=255)             gcgcGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGa                                                                                                                                      <  7:306886/149‑1 (MQ=255)              cgcgTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTg                                                                                                                                       >  8:23977/1‑147 (MQ=255)                   gttgcCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTAtt                                                                                                                                <  4:347328/146‑1 (MQ=255)                       ccAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGa                                                                                                                              >  8:31199/1‑147 (MQ=255)                        cAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAg                                                                                                                           >  5:183590/1‑149 (MQ=255)                              gggCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGAATTTCATGCGGCCTTCaaa                                                                                                                                    >  5:137250/1‑133 (MQ=255)                              gggCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGAATTTCATGCGGCCTTCaaa                                                                                                                                    <  6:137250/134‑2 (MQ=255)                                gcggcCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTAcgc                                                                                                                   >  3:240521/1‑149 (MQ=255)                                     cTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCg                                                                                                               >  2:207409/1‑148 (MQ=255)                                        gACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGt                                                                                                           >  5:120838/1‑149 (MQ=255)                                               tGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTgac                                                                                                    <  4:240521/149‑1 (MQ=255)                                               gggAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTgac                                                                                                    <  5:157100/148‑1 (MQ=255)                                                 taaTGGTCGCCAGTGCACCAGCACATTATCTACCTGTGCTGTAGCCAACGGTTTCAGAAGCGTACACATATTTGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTgacga                                                                                                  <  1:305407/148‑1 (MQ=255)                                                 gAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACTAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGaaa                                                                                                                            <  3:114707/123‑1 (MQ=255)                                                 gAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGaaa                                                                                                                            >  4:114707/1‑123 (MQ=255)                                                       tCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCt                                                                                             <  6:161128/148‑1 (MQ=255)                                                       tCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTg                                                                                            >  4:78304/1‑149 (MQ=255)                                                          gCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGt                                                                                         <  1:207409/149‑1 (MQ=255)                                                            aGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTgg                                                                                       <  7:345938/149‑1 (MQ=255)                                                              tGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTggg                                                                                      >  3:82791/1‑148 (MQ=255)                                                                  ccAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGc                                                                                 <  7:86411/149‑1 (MQ=255)                                                                     gCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTa                                                                              <  6:183590/149‑1 (MQ=255)                                                                      cacaTGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTa                                                                                                                      >  8:13414/1‑108 (MQ=255)                                                                      cacaTGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTa                                                                                                                      <  7:13414/108‑1 (MQ=255)                                                                         aTGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTg                                                                          >  5:277681/1‑149 (MQ=255)                                                                            aTCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCgg                                                                       <  7:23977/149‑1 (MQ=255)                                                                             tCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCAt                                                                                   >  5:130863/1‑136 (MQ=255)                                                                             tCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCAt                                                                                   <  6:130863/136‑1 (MQ=255)                                                                             tCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGt                                                                      >  2:236689/1‑149 (MQ=255)                                                                               gACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGa                                                                    >  8:258569/1‑149 (MQ=255)                                                                                aCCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGAc                                                                   <  7:31199/149‑1 (MQ=255)                                                                                aCCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGAc                                                                   >  4:153714/1‑149 (MQ=255)                                                                                aCCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAACCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGAc                                                                   <  1:288640/149‑1 (MQ=255)                                                                                   aGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAg                                                                >  6:320202/1‑149 (MQ=255)                                                                                    gTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGt                                                               >  6:6247/1‑149 (MQ=255)                                                                                    gTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGt                                                               <  3:167569/149‑1 (MQ=255)                                                                                       cTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTc                                                             >  8:157493/1‑148 (MQ=255)                                                                                        tGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCa                                                           >  5:89170/1‑149 (MQ=255)                                                                                          gaccccTCGGGTCACCAAGCGTACCTACTCAGGGGGCGCCACATTTGCAGGCGGCATGAAATTCGGCCTTTACGTTTTGAAAGCCGTACGCGCCGTCGTGGGTGACGACTTCTGGGTGGGCATGCGTATCTGCGGTGACGAGTTcccc                                                          <  6:277681/145‑3 (MQ=255)                                                                                           aGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCAccc                                                        <  6:120838/149‑1 (MQ=255)                                                                                               cACGGGTCAACAATCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTTAAAGCCGTACGAGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCCACCTGAt                                                    >  3:209513/1‑149 (MQ=255)                                                                                                aCGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATg                                                   <  3:78304/149‑1 (MQ=255)                                                                                                  gggTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGc                                                 <  4:82791/149‑1 (MQ=255)                                                                                                       aacaaGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCa                                             <  7:258569/148‑1 (MQ=255)                                                                                                       aacaaGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCa                                             <  6:205548/148‑1 (MQ=255)                                                                                                       aacaaGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTAGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAg                                            >  8:334123/1‑149 (MQ=255)                                                                                                             cGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGa                                      <  3:153714/149‑1 (MQ=255)                                                                                                                aCCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCa                                                           >  3:301098/1‑125 (MQ=255)                                                                                                                aCCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCa                                                           <  4:301098/125‑1 (MQ=255)                                                                                                                  cgacgaATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACa                                 <  4:17062/149‑1 (MQ=255)                                                                                                                  agaagaATGGGGCGGAAACTTTGAAGGCCTCATTAAATTCGGCCTTGAAGTATTGAAAGCAGTACGCGCCGAATTGTGTGACGACTTCTGCGTGGGAATGAGTATCTGCGGTGACGAGTTCCACCATGATTGCCTCAGCCACGAGGACa                                 >  2:65494/5‑149 (MQ=255)                                                                                                                           gggCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGa                        <  6:89170/149‑1 (MQ=255)                                                                                                                                 cACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTg                  >  7:289839/1‑149 (MQ=255)                                                                                                                                   tacttcaaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGcc                >  1:79357/7‑149 (MQ=255)                                                                                                                                   ccTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGcc                >  5:353219/1‑149 (MQ=255)                                                                                                                                      ttcaaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGAt                                                    <  7:221326/107‑1 (MQ=255)                                                                                                                                      ttcaaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGAt                                                    >  8:221326/4‑110 (MQ=255)                                                                                                                                         aaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCtactac          >  7:307098/1‑149 (MQ=255)                                                                                                                                                 aTGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAc  <  8:289839/149‑1 (MQ=255)                                                                                                                                                   gAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCctacta           >  3:162936/1‑138 (MQ=255)                                                                                                                                                   gAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCctacta           <  4:162936/138‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAC  >  NC_002947/370973‑371265

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: ‑