Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	932,603	A→G	6.3%	V4340A (GTG→GCG) ‡	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	932,603	0	A	G	6.3%	178.8 / 2.3	63	V4340A (GTG→GCG) ‡	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGG  >  NC_002947/932458‑932744                                                                                                                                                  |                                                                                                                                              cACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGt                                                                                                                                            <  3:70150/149‑1 (MQ=255)       aCCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCgg                                                                                                                                      >  4:236712/1‑149 (MQ=255)        ccGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCgg                                                                                                                                      <  3:98804/148‑1 (MQ=255)          gTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCAACGTGTCGCCGGccg                                                                                                                                   <  2:205273/149‑3 (MQ=255)          gTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCAACGTGTCGCCGGcc                                                                                                                                    >  1:210528/1‑147 (MQ=255)          gTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCAACGTGTCGCCGGcc                                                                                                                                    <  2:210528/148‑2 (MQ=255)           tCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGc                                                                                                                                  <  3:236712/149‑1 (MQ=255)               gTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGACTGCGTCg                                                                                                                              <  2:71884/149‑1 (MQ=255)                  cTCGTACCGTCTCTGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGTCTGCCAGGTTGGCAACTGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGt                                                                                                                            <  3:101821/148‑1 (MQ=255)                   tCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTgg                                                                                                                          >  8:330198/1‑149 (MQ=255)                     gTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCg                                                                                                                        <  1:92248/149‑1 (MQ=255)                       aCCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGt                                                                                                                      >  6:97245/1‑149 (MQ=255)                        ccGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGt                                                                                                                      >  1:173947/1‑148 (MQ=255)                         cGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCa                                                                                                                    >  1:240164/1‑149 (MQ=255)                               ggTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCa                                                                                                              <  2:240164/149‑1 (MQ=255)                               ggTTGCGCCTGCTTCGCCATCACCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCa                                                                                                              <  2:212607/149‑1 (MQ=255)                                gTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGc                                                                                                                                  <  7:54578/128‑1 (MQ=255)                                gTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGc                                                                                                                                  >  8:54578/1‑128 (MQ=255)                                     gcCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCt                                                                                                        >  5:36443/1‑149 (MQ=255)                                      ccTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTg                                                                                                       <  4:58638/149‑1 (MQ=255)                                       cTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTgg                                                                                                      <  3:340461/149‑1 (MQ=255)                                        tGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGtt                                                                                                                                           <  6:197238/111‑1 (MQ=255)                                        tGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGtt                                                                                                                                           >  5:197238/1‑111 (MQ=255)                                          cTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGc                                                                                                     <  6:130382/147‑1 (MQ=255)                                              gCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTgg                                                                                               >  2:246994/1‑149 (MQ=255)                                               ccATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGt                                                                                              <  7:330198/149‑1 (MQ=255)                                                cATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGGCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTc                                                                                             >  1:252944/1‑149 (MQ=255)                                                cATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTg                                                                                                                >  8:260572/1‑130 (MQ=255)                                                cATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTg                                                                                                                <  7:260572/130‑1 (MQ=255)                                                     ccTGTCAGGGTGGCGCCATTGGCGGAGAAAGCCAGGTTGGCCACCGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGAGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTCGTCAggg                                                                                         >  8:43716/1‑148 (MQ=255)                                                           aGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGtt                                                                                  >  5:290553/1‑149 (MQ=255)                                                                       ttGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTg                                                                      >  5:65558/1‑149 (MQ=255)                                                                          gCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGaa                                                                   <  2:64198/149‑1 (MQ=255)                                                                            ggAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTg                                                                                                       >  5:144594/1‑111 (MQ=255)                                                                            ggAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTg                                                                                                       <  6:144594/111‑1 (MQ=255)                                                                                   gCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGtt                                                          <  5:10064/149‑1 (MQ=255)                                                                                   gCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGtt                                                          <  6:290553/149‑1 (MQ=255)                                                                                     cAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGc                                                              <  6:46967/143‑1 (MQ=255)                                                                                     cAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGc                                                              >  5:46967/1‑143 (MQ=255)                                                                                     cAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGc                                                        <  1:246994/149‑1 (MQ=255)                                                                                           ggCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCgg                                                  >  8:63556/1‑149 (MQ=255)                                                                                             ccACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGc                                                 >  7:216060/1‑148 (MQ=255)                                                                                              cACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTgg                                                                                          <  3:16015/106‑1 (MQ=255)                                                                                              cACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTa                                               >  6:356863/1‑149 (MQ=255)                                                                                              cACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTAGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTgg                                                                                          >  4:16015/1‑106 (MQ=255)                                                                                                  gcggccggcggcGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGc                                           >  6:196125/1‑149 (MQ=255)                                                                                                     gccggcggcGTTGTGTCCGCTGCTGTGTAAGGTGCGGCGGGCGACAGGTTGCCGGCTGGGTCTGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCtgt                                        <  7:43716/149‑1 (MQ=255)                                                                                                            gcgTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGc                                 <  5:356863/149‑1 (MQ=255)                                                                                                                ggtgTCCGCTGCGGTGTAAGGTGAGGCCGGCGACAAGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTTGTCTGGGCACTGTTCAAGGTGACCTGGAAGTTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGt                             <  6:58521/149‑1 (MQ=255)                                                                                                                 gtgtCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGtt                            >  3:10932/1‑149 (MQ=255)                                                                                                                       gCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCt                      <  2:252944/149‑1 (MQ=255)                                                                                                                         tGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGaa                                                                   <  2:102019/102‑1 (MQ=255)                                                                                                                         tGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGaa                                                                   >  1:102019/1‑102 (MQ=255)                                                                                                                            ggTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCAGTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCAt                 >  4:92914/1‑149 (MQ=255)                                                                                                                             gtgtAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATc                >  1:152623/1‑149 (MQ=255)                                                                                                                             gtgtAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGGTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATc                >  7:318619/1‑149 (MQ=255)                                                                                                                              tgtAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATcg               >  8:129338/1‑149 (MQ=255)                                                                                                                                tAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATcgc              <  6:65558/148‑1 (MQ=255)                                                                                                                                   ggTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcac          <  2:173947/149‑1 (MQ=255)                                                                                                                                        ggccggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtca     >  3:236292/1‑149 (MQ=255)                                                                                                                                        agccggcAACGTGTCGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtca     >  2:239263/2‑149 (MQ=255)                                                                                                                                          ccggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtcacg   <  5:97245/149‑1 (MQ=255)                                                                                                                                           cggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtcacgg  <  4:10932/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              CACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGG  >  NC_002947/932458‑932744

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