Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	991,557	T→C	7.6%	I60T (ATC→ACC)	der →	GTPase Der

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	991,557	0	T	C	7.6%	191.6 / 5.1	67	I60T (ATC→ACC)	der	GTPase Der
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (2/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

TCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAG  >  NC_002947/991424‑991690                                                                                                                                      |                                                                                                                                      tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGATACCACCGGTGTCGatca                                                                                                                        <  6:63387/149‑4 (MQ=255) tCCACCAGGTTCAACCGCCTGACCACGACGCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTAGGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCTGTGATACCACCGGTGTCGatca                                                                                                                        <  7:203614/149‑4 (MQ=255)    accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                     >  5:61955/1‑149 (MQ=255)        tGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTAt                                                                                                                                       >  6:157143/1‑127 (MQ=255)        tGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTAt                                                                                                                                       <  5:157143/127‑1 (MQ=255)               ccGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATgg                                                                                                                >  8:56771/1‑143 (MQ=255)               ccGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATgg                                                                                                                <  7:56771/143‑1 (MQ=255)                  ccTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaaaaa                                                                                                       >  4:141963/1‑149 (MQ=255)                        cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATg                                                                                                     <  1:175026/145‑1 (MQ=255)                        cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATg                                                                                                     >  2:175026/1‑145 (MQ=255)                        cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATAACCGGTGACGAAGTGGGAATGGACGAAAAAATTGCCg                                                                                                 >  4:222251/1‑149 (MQ=255)                          aGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                               >  6:19673/1‑149 (MQ=255)                             cccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAg                                                                                            >  2:1521/1‑149 (MQ=255)                              ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGt                                                                                           >  3:190757/1‑149 (MQ=255)                               cgcgATGCCATCGTCGCTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTAACACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGt                                                                                           >  8:266111/1‑148 (MQ=255)                                      ccATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATg                                                                                   <  7:129005/149‑1 (MQ=255)                                      ccATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATg                                                                                   <  3:141963/149‑1 (MQ=255)                                       cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATgg                                                                                                    >  5:38993/1‑131 (MQ=255)                                       cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATgg                                                                                                    <  6:38993/131‑1 (MQ=255)                                       cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATgg                                                                                  <  4:190757/149‑1 (MQ=255)                                        aTCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                 >  4:45676/1‑149 (MQ=255)                                           gtcgGTGACCTGTCGGGCCGGACCCGTGACCGCCAGTATGGTGATGCCGGCTGGCATGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCAt                                                                              <  6:89004/149‑1 (MQ=255)                                             cgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTACTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGCACAtt                                                                             >  7:21134/1‑148 (MQ=255)                                              ggTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTg                                                                            <  5:234010/148‑1 (MQ=255)                                               gTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCCGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaa                                                                          <  5:143308/149‑1 (MQ=255)                                                tGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaag                                                                         <  3:45676/149‑1 (MQ=255)                                                  aCCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGTCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaag                                                                         >  1:233399/1‑147 (MQ=255)                                                    cTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                 >  3:118641/1‑137 (MQ=255)                                                    cTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                 <  4:118641/137‑1 (MQ=255)                                                       tCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCg                                                                   <  1:1521/148‑1 (MQ=255)                                                         gggCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGAc                                                                 >  5:163750/1‑148 (MQ=255)                                                          ggCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTa                                                               >  3:144256/1‑149 (MQ=255)                                                           gCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAt                                                              >  6:182771/1‑149 (MQ=255)                                                                 cccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGTTg                                                        >  1:275790/1‑149 (MQ=255)                                                                  ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGt                                                       >  2:106880/1‑149 (MQ=255)                                                                     tGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGTAGAAGCCGACTATGTGCTGTTcc                                                    >  2:244720/1‑149 (MQ=255)                                                                     cctatctcAGTATGGTGATGGCAGCTGGTAGGGTCGTTCCTTCATGCTTATCGACACCGGTGGTATCACGGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTc                                                     <  8:60551/141‑1 (MQ=255)                                                                      gACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCt                                                   >  6:138885/1‑149 (MQ=255)                                                                      gACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCt                                                   >  3:225109/1‑149 (MQ=255)                                                                       aCCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCt                                                   <  5:138885/148‑1 (MQ=255)                                                                       aCCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTg                                                  <  1:106880/149‑1 (MQ=255)                                                                        ccGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTgg                                                 <  5:182771/149‑1 (MQ=255)                                                                           ccAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                               >  4:44148/1‑148 (MQ=255)                                                                            cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGa                                             <  2:275790/149‑1 (MQ=255)                                                                             aGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTc                                                     >  7:60551/1‑140 (MQ=255)                                                                                  ggTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCg                                       >  1:218165/1‑149 (MQ=255)                                                                                      aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGTAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGc                                    <  1:244720/148‑1 (MQ=255)                                                                                      aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTc                                                     <  2:14777/131‑1 (MQ=255)                                                                                      aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTc                                                     >  1:14777/1‑131 (MQ=255)                                                                                      aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCg                                   <  5:19673/149‑1 (MQ=255)                                                                                      aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCg                                   >  8:305261/1‑149 (MQ=255)                                                                                         ccAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGt                                <  4:144256/149‑1 (MQ=255)                                                                                         ccAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGt                                <  2:218165/149‑1 (MQ=255)                                                                                          cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTgtc                                                                                                                             >  3:215765/1‑53 (MQ=21)                                                                                          cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTgtc                                                                                                                             <  4:215765/55‑3 (MQ=21)                                                                                          cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGt                                                >  6:59784/1‑132 (MQ=255)                                                                                          cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGt                                                <  5:59784/132‑1 (MQ=255)                                                                                           aGCTTGCAGGGTCGTTACTTCATCCTGATCGACACCTGTGGTATCACCGGTGACTAAGTGGGCATGGACGATAAAATGGCCGAGCAGTCGCTCTTGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATTCCCGTGCGGGTAt                              <  7:266111/149‑1 (MQ=255)                                                                                               ggCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACt                          >  7:285620/1‑149 (MQ=255)                                                                                                   gggTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                  >  5:325179/1‑137 (MQ=255)                                                                                                   gggTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                  <  6:325179/137‑1 (MQ=255)                                                                                                   gggTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCgg                      <  6:163750/149‑1 (MQ=255)                                                                                                    ggTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaa                                                                          <  1:312832/96‑1 (MQ=255)                                                                                                    ggTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaa                                                                          >  2:312832/1‑96 (MQ=255)                                                                                                            cTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAga             <  2:233399/149‑1 (MQ=255)                                                                                                                   cTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGc      >  4:304685/1‑149 (MQ=255)                                                                                                                    tGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGcc     >  6:158419/1‑149 (MQ=255)                                                                                                                       tCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAg  <  8:21134/149‑1 (MQ=255)                                                                                                                        cgtcacCGGTGATACCTCCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGTGTATTACTGCGGCAGACCAGATGATCGCCGAg  >  8:267775/4‑148 (MQ=255)                                                                                                                                      |                                                                                                                                      TCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAG  >  NC_002947/991424‑991690

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