Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	371,119	G→T	11.1%	M200I (ATG→ATT)	dgcA →	dimethylglycine dehydrogenase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	371,119	0	G	T	11.1%	190.0 / 11.0	63	M200I (ATG→ATT)	dgcA	dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base T (5/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAC  >  NC_002947/370971‑371265                                                                                                                                                     |                                                                                                                                                   aCTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATg                                                                                                                                                    <  3:11040/149‑1 (MQ=255)   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)                                                                                                                                                     |                                                                                                                                                   ACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAC  >  NC_002947/370971‑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: ‑