Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	790,620	T→C	6.6%	Q17R (CAG→CGG)	PP_0676 ←	transcription elongation factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	790,620	0	T	C	6.6%	174.1 / 2.9	62	Q17R (CAG→CGG)	PP_0676	transcription elongation factor
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/2);  total (30/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 = 4.36e-01

GGCCAGCCGTTGTTTGTCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTGGCGATAGGCCTTGG  >  NC_002947/790475‑790759                                                                                                                                                   |                                                                                                                                             ggCCAGCCGTTGTTTGTCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGAt                                                                                                                                           >  5:258855/1‑149 (MQ=255) ggCCAGCCGTTGTTTGTCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGAt                                                                                                                                           <  7:60246/149‑1 (MQ=255)  gCCAGCCGTTGTTTGTCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTG‑ACAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCg                                                                                                                                         >  1:78826/1‑149 (MQ=255)     aGCCGTTGTTTGTCACCCCGTTCGCCCTGCCCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGc                                                                                                                                       >  6:106527/1‑149 (MQ=255)      gCCGTTGTTTGTCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGcc                                                                                                                                      <  5:82977/149‑1 (MQ=255)             tttgtCACCCGGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCTCATAGTTTGGCTGTTCGCTGATGCGCCGCTCGTCC‑GGCTGATCGGCCTGGTcg‑g                                                                                                                              <  5:259158/149‑1 (MQ=255)             tttgtCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTcg                                                                                                                                <  2:83222/148‑1 (MQ=255)              ttgtCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTcg‑gc                                                                                                                             >  4:91170/1‑149 (MQ=255)                gtCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTcg‑gcgg                                                                                                                           >  2:215315/1‑149 (MQ=255)                  cACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCt                                                                                                                         >  7:260947/1‑149 (MQ=255)                    ccccGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTg                                                                                                                       <  7:275155/149‑1 (MQ=255)                    ccccGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTg                                                                                                                       <  7:33550/149‑1 (MQ=255)                       cGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTc                                                                                                                    >  1:76882/1‑149 (MQ=255)                              ccTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTg                                                                                                                       <  4:9655/139‑1 (MQ=255)                              ccTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTg                                                                                                                       >  3:9655/1‑139 (MQ=255)                              ccTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGC‑CGACCAGGCCGATCAGCC‑GGTCGAgcggc                                                                                                                          <  2:237893/136‑1 (MQ=38)                              ccTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGC‑CGACCAGGCCGATCAGCC‑GGTCGAgcggc                                                                                                                          <  2:314711/136‑1 (MQ=38)                              ccTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGC‑CGACCAGGCCGATCAGCC‑GGTCGAgcggc                                                                                                                          >  1:314711/1‑136 (MQ=38)                              ccTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGC‑CGACCAGGCCGATCAGCC‑GGTCGAgcggc                                                                                                                          >  1:237893/1‑136 (MQ=38)                               cTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGa                                                                                                            <  8:21713/149‑1 (MQ=255)                                   gcCTGCAGTTCGCTGCGCAAGGCATTCAGCGCGGCCACCCGCGTTTGCAACTGCCGCATGCCGCTGGCAGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTGCGTTGACGaa                                                                                                        <  3:26475/149‑1 (MQ=255)                                    ccTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAAt                                                                                                       <  3:250664/149‑1 (MQ=255)                                     cTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATg                                                                                                      >  3:147635/1‑149 (MQ=255)                                     cTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATg                                                                                                      >  2:224476/1‑149 (MQ=255)                                         aGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAAAGCTCg                                                                                                  >  3:239974/1‑149 (MQ=255)                                          gTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCgg                                                                                                 >  6:175615/1‑149 (MQ=255)                                           ttCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGc                                                                                                >  7:198620/1‑149 (MQ=255)                                               cTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGCCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCAt                                                                                            <  4:141610/149‑1 (MQ=255)                                                    cAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAgg                                                                                       <  2:78826/149‑1 (MQ=255)                                                     aaGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAggg                                                                                      >  1:316831/1‑149 (MQ=255)                                                      aGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGc                                                                                     <  7:80072/149‑1 (MQ=255)                                                                  cGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCt                                                                                                                   <  5:74774/107‑1 (MQ=255)                                                                  cGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCt                                                                                                                   >  6:74774/1‑107 (MQ=255)                                                                    gCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATga                                                                       <  8:260947/149‑1 (MQ=255)                                                                      cagccGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGGCTGGTCG‑TCGGCTTGGTCTTCGTTGACGAATGCTCTGCTCATGTAGGGCCTCCTTGTGCATgaga                                                                     <  7:191253/146‑1 (MQ=255)                                                                      cACCCGCGTTTGGAACTGGCGCAGGCCTCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGAGCATgaga                                                                     <  5:106527/149‑1 (MQ=255)                                                                      cACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATgaga                                                                     <  6:169792/149‑1 (MQ=255)                                                                      cACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATgaga                                                                     >  7:34530/1‑149 (MQ=255)                                                                        cccGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCgg                                                                                                 <  7:29179/119‑1 (MQ=255)                                                                        cccGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCgg                                                                                                 >  8:29179/1‑119 (MQ=255)                                                                           gcgTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACcat                                                                 >  6:132484/1‑148 (MQ=255)                                                                              tttGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATg                                                             <  1:224476/149‑1 (MQ=255)                                                                              tttGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATg                                                             <  6:258855/149‑1 (MQ=255)                                                                              tttGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGACCTCCTTGTGCATGAGACCATCATg                                                             >  6:27415/1‑149 (MQ=255)                                                                                  cAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCa                                                          <  1:215315/148‑1 (MQ=255)                                                                                   aaCTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCgg                                                                                                 <  6:310306/108‑1 (MQ=255)                                                                                   aaCTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCgg                                                                                                 >  5:310306/1‑108 (MQ=255)                                                                                      tGGCGCAGGCCGCTGGACGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTAGGCTCATGTAGGGCATCATTGTGAATGAGACCATCATAGAAGgccg                                                     >  8:246814/1‑149 (MQ=255)                                                                                          gcAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGTCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCgg                                                 >  3:48749/1‑149 (MQ=255)                                                                                           cAGGCCGCTGGACGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATAGCAGGCCGCCGGt                                                <  7:246814/149‑1 (MQ=255)                                                                                                 gCTGGCCGTCTCTTAGTTGGGCTGTTCGGTGATGGGCCGCTCGACC‑GGCTGATCGGCCTTGTGG‑GCGGCTTTGTGTTCGTTGACGAATGCTCGGCTCATGTAGGTCCTCCTTGTGCATGAGAACATCATGGCAGGCCGCCGGTTTAAAg                                          <  4:147635/149‑1 (MQ=255)                                                                                                         tCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTcc                                  <  3:91170/149‑1 (MQ=255)                                                                                                             aTAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGACTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGttt                              >  8:165811/1‑149 (MQ=255)                                                                                                                      cTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCaa                     <  2:316831/149‑1 (MQ=255)                                                                                                                            gCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTgg               <  8:34530/149‑1 (MQ=255)                                                                                                                                    gccgcTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGTTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTGGCGATAGGc       >  2:250791/1‑149 (MQ=255)                                                                                                                                      cgcTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTGGCGATAGGCCt     >  6:248110/1‑149 (MQ=255)                                                                                                                                       gcTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTTAAGGAATGTCCGTTTGCGACTCAATTGTGGAGATAGGCCtt    >  5:89919/1‑149 (MQ=255)                                                                                                                                         tCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTaaa                                           <  2:167663/108‑1 (MQ=255)                                                                                                                                         tCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTaaa                                           >  1:167663/1‑108 (MQ=255)                                                                                                                                         tCGACC‑GGATGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGGCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTGGCGATAGGCCTTgg  >  6:124875/1‑149 (MQ=255)                                                                                                                                         gcGACC‑GGCGGATCGGCCTGGGCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTGGCGATAGGCCTTgg  <  2:76882/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             GGCCAGCCGTTGTTTGTCACCCCGTTCGCCCTGCGCCTGCAGTTCGCTGCGCAAGGCATTCAGCTCGGCCACCCGCGTTTGCAACTGGCGCAGGCCGCTGGCCGTCACATAGTTGGGCTGTTCGCTGATGCGCCGCTCGACC‑GGCTGATCGGCCTGGTCG‑GCGGCTTGGTCTTCGTTGACGAATGCTCGGCTCATGTAGGGCCTCCTTGTGCATGAGACCATCATGGCAGGCCGCCGGTTTAAAGGAATGTCCGTTTGCGACTCAATTGTGGCGATAGGCCTTGG  >  NC_002947/790475‑790759

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