Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	713,097	T→A	7.5%	I29N (ATC→AAC)	lspA →	lipoprotein signal peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	713,097	0	T	A	7.5%	188.7 / 3.2	63	I29N (ATC→AAC)	lspA	lipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (2/3);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTG  >  NC_002947/712952‑713244                                                                                                                                                  |                                                                                                                                                     gCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGa                                                                                                                                                   >  4:41383/1‑149 (MQ=255)    gCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCa                                                                                                                                                <  5:38720/149‑1 (MQ=255)    gCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCGGAGCTTGCTGGTCCTGGGCAGCgcc                                                                                                                                                  >  7:208509/1‑145 (MQ=255)     ccATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAg                                                                                                                                               >  6:82856/1‑149 (MQ=255)     ccATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAg                                                                                                                                               <  6:54810/149‑1 (MQ=255)     ccATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAg                                                                                                                                               <  3:84109/149‑1 (MQ=255)          cGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑c                                                                                                                                            >  8:194219/1‑146 (MQ=255)          cGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑c                                                                                                                                            <  7:194219/146‑1 (MQ=255)              gCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAg                                                                                                                                     <  6:51804/149‑1 (MQ=255)               cACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGG‑CCTGGTCGATGACCAGGACCAGCAAGc                                                                                                                                    <  6:258330/149‑1 (MQ=38)                aCCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGAc                                                                                                                                                  <  1:170808/135‑1 (MQ=255)                aCCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGAc                                                                                                                                                  >  2:170808/1‑135 (MQ=255)                  ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGG‑CCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                                 <  8:44224/149‑3 (MQ=38)                  ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTg                                                                                                                                  >  7:307030/1‑148 (MQ=255)                     gAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACt                                                                                                                              >  1:139569/1‑149 (MQ=255)                        aTCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTca                                                                                                                           <  6:198690/149‑1 (MQ=255)                              ggCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACg                                                                                                                     <  8:252784/149‑1 (MQ=255)                               gCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGc                                                                                                                    >  5:20766/1‑149 (MQ=255)                                   ttGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTAc                                                                                                                               >  3:263546/1‑134 (MQ=255)                                   ttGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTAc                                                                                                                               <  4:263546/134‑1 (MQ=255)                                       gTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCGTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCa                                                                                                            <  8:288082/149‑1 (MQ=255)                                          gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGt                                                                                                         <  5:99341/149‑1 (MQ=255)                                                cTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagc                                                                                                   <  1:185282/149‑1 (MQ=255)                                                              tGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCAtt                                                                                     >  1:116834/1‑149 (MQ=255)                                                              tGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCAtt                                                                                     >  6:83307/1‑149 (MQ=255)                                                               gAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                                    >  4:99410/1‑149 (MQ=255)                                                                aGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                                       >  3:178352/1‑129 (MQ=255)                                                                aGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                                       <  4:178352/129‑1 (MQ=255)                                                                aGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTcc                                                                                   <  8:307030/149‑1 (MQ=255)                                                                       ccACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGactac                                                                            >  4:89839/1‑149 (MQ=255)                                                                           tATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACtt                                                                          >  1:9179/1‑147 (MQ=255)                                                                             tGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCa                                                                        <  2:9179/147‑1 (MQ=255)                                                                             tGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAg                                                                       <  5:61267/148‑1 (MQ=255)                                                                               ccTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                           >  6:139207/1‑126 (MQ=255)                                                                               ccTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCt                                                                     <  6:16255/148‑1 (MQ=255)                                                                                 tAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCAACGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTgaa                                                                  >  5:319030/1‑147 (MQ=255)                                                                                   aCCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGa                                                                  <  2:139569/147‑1 (MQ=255)                                                                                     ccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCt                                                              >  1:172977/1‑149 (MQ=255)                                                                                     ccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCt                                                              <  3:251734/149‑1 (MQ=255)                                                                                        gCAGCGTGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCACCAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGc                                                           <  7:217439/149‑1 (MQ=255)                                                                                             ggggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTac                                                       >  3:36754/1‑148 (MQ=255)                                                                                                         gcgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcg                                          >  4:9480/1‑149 (MQ=255)                                                                                                          cgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGTTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTgg                                                                                          >  6:161676/1‑100 (MQ=255)                                                                                                          cgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGTTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTgg                                                                                          <  5:161676/100‑1 (MQ=255)                                                                                                          cgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcgc                                         <  5:83307/149‑1 (MQ=255)                                                                                                          cgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcgc                                         >  1:207249/1‑149 (MQ=255)                                                                                                             cTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGtt                                      <  2:116834/149‑1 (MQ=255)                                                                                                                gCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACTACACCGGTGCCGCGTTCa                                    <  5:82856/148‑1 (MQ=255)                                                                                                                      ctgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCt                             >  4:321990/1‑149 (MQ=255)                                                                                                                         tggctgAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTg                            <  2:231579/147‑1 (MQ=255)                                                                                                                         tggctgAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGc                          <  1:60040/149‑1 (MQ=255)                                                                                                                              caGCTTGGTGG‑CCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCt                             >  7:18064/2‑141 (MQ=38)                                                                                                                              caGCTTGGTGG‑CCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCt                             <  8:18064/140‑1 (MQ=38)                                                                                                                                  ttGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGc                 >  3:117168/1‑149 (MQ=255)                                                                                                                                    gCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGc                 <  3:99410/147‑1 (MQ=255)                                                                                                                                      tggtccttgtcATCGACCAGG‑CCAACAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTacaa                                                     <  7:88447/109‑1 (MQ=255)                                                                                                                                      tggtcctggtcATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGt             >  5:179636/1‑149 (MQ=255)                                                                                                                                      tggtcctggtcATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTACTGACTACTTCAGCTGGACCCTGGATTaaaa                                                     >  8:88447/1‑109 (MQ=255)                                                                                                                                       ggtcctggtcATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGtt            >  8:197705/1‑149 (MQ=255)                                                                                                                                          cctggtcATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTgg          >  2:267686/1‑148 (MQ=255)                                                                                                                                                 aTCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCt                         >  4:217764/1‑126 (MQ=255)                                                                                                                                                 aTCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCt                         <  3:217764/126‑1 (MQ=255)                                                                                                                                                 aTCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTg  >  8:119176/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGG‑CCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTG  >  NC_002947/712952‑713244

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