Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,006,474	A→C	15.0%	K1019Q (AAG→CAG)	PP_2627 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,006,474	0	A	C	15.0%	149.5 / 21.4	60	K1019Q (AAG→CAG)	PP_2627	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/19);  new base C (6/3);  total (38/22)
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.91e-01

TACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAGCC  >  NC_002947/3006330‑3006622                                                                                                                                                 |                                                                                                                                                      tACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAgg                                                                                                                                                    <  5:182929/148‑1 (MQ=255) tACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGa                                                                                                                                                   >  5:110140/1‑149 (MQ=255)    aGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACaa                                                                                                                                                >  6:316309/1‑149 (MQ=255)             cAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGc                                                                                                                                       >  4:129413/1‑149 (MQ=255)                cAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCtt                                                                                                                                    <  5:316309/149‑1 (MQ=255)                   cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                                      >  6:24378/1‑144 (MQ=255)                   cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                                      <  5:24378/144‑1 (MQ=255)                                     cccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGc                                                                                                                >  6:18186/1‑148 (MQ=255)                                       cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAAt                                                                                                             >  3:68047/1‑149 (MQ=255)                                               tgAAGGCAACCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑tt                                                                                                    <  3:283133/149‑1 (MQ=255)                                                gAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTc                                                                                                   >  2:52741/1‑149 (MQ=255)                                                   ggCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTg                                                                                                <  5:18186/149‑1 (MQ=255)                                                     cACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGcc                                                                                              >  4:59940/1‑149 (MQ=255)                                                      aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGcc                                                                                              >  4:59226/1‑148 (MQ=255)                                                      aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGcc                                                                                              >  1:306026/1‑148 (MQ=255)                                                         cGATGCTTCTTTGGATGACTTCGAGGCCTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGAACGGGCAATCCATGA‑TTCGTGCCGACg                                                                                          >  7:160145/1‑149 (MQ=255)                                                             gCTTATTTGGATGACTTCGTGGCGTTCTTCGGGCCTAATGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGa                                                                                       <  2:306026/148‑1 (MQ=255)                                                                  tttGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGc                                                                                               <  7:1223/135‑1 (MQ=255)                                                                  tttGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGc                                                                                               >  8:1223/1‑135 (MQ=255)                                                                    tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCt                                                                                                         <  8:214450/124‑1 (MQ=255)                                                                    tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCt                                                                                                         >  7:214450/1‑124 (MQ=255)                                                                    tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGagc                                                                                 <  2:144181/147‑1 (MQ=255)                                                                      gATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCt                                                                             >  5:145723/1‑149 (MQ=255)                                                                      gATGACTTCGAGGCGTTCTTCCGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCt                                                                             >  8:101279/1‑149 (MQ=255)                                                                       aTGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTc                                                                            >  5:239769/1‑149 (MQ=255)                                                                        tGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTc                                                                            >  3:124672/1‑148 (MQ=255)                                                                            ttCGATGCGTTCTTCGGGCCGAAAGGTCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTGCTGAAGTACAACCTTGATGCCTTGCAGTCCGCTCAGCACTGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCt                                                                       <  8:160145/149‑1 (MQ=255)                                                                            ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCt                                                                       >  3:299500/1‑149 (MQ=255)                                                                                aGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCg                                                                   <  4:124672/149‑1 (MQ=255)                                                                                   cGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGa                                                                 <  1:52741/149‑1 (MQ=255)                                                                                    gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                               >  2:104690/1‑149 (MQ=255)                                                                                        ttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGcgcggag                                                           >  6:206904/1‑146 (MQ=255)                                                                                        ttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCAt                                                           >  2:177565/1‑149 (MQ=255)                                                                                           gggCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATcc                                                         <  6:239769/148‑1 (MQ=255)                                                                                                      gcgcCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACtttct                                             <  4:299500/149‑1 (MQ=255)                                                                                                                cagTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCg                                   >  4:243222/1‑149 (MQ=255)                                                                                                                      aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                               <  4:199796/115‑1 (MQ=255)                                                                                                                      aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                               >  3:199796/1‑115 (MQ=255)                                                                                                                      aaTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGc                               >  2:193093/1‑147 (MQ=255)                                                                                                                       aTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAg                            <  5:206904/149‑1 (MQ=255)                                                                                                                        tGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCg                                                                   >  2:40861/1‑109 (MQ=255)                                                                                                                        tGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCg                                                                   <  1:40861/109‑1 (MQ=255)                                                                                                                          aTCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCt                          <  6:110140/148‑1 (MQ=255)                                                                                                                                  ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCg                                                                           <  5:174878/91‑1 (MQ=38)                                                                                                                                  ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCg                                                                           >  6:174878/1‑91 (MQ=38)                                                                                                                                  ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg                   >  2:68595/1‑147 (MQ=255)                                                                                                                                   tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg                >  1:307809/1‑149 (MQ=255)                                                                                                                                   tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg                >  8:110628/1‑149 (MQ=255)                                                                                                                                   tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                            <  6:218507/121‑1 (MQ=255)                                                                                                                                   tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                            >  5:218507/1‑121 (MQ=255)                                                                                                                                   tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg                >  3:114634/1‑149 (MQ=255)                                                                                                                                   tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGGGATCGAGCAGCTCGCGCTGGCGGAGCGCATCAGCGAAACCTTCTGCGCTCAGCGTGGCAAGCTAAGTGTGCAg                >  8:283981/1‑149 (MQ=255)                                                                                                                                    cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt               >  1:191833/1‑149 (MQ=255)                                                                                                                                    cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt               >  1:178306/1‑149 (MQ=255)                                                                                                                                     aaGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt               >  7:96686/1‑148 (MQ=255)                                                                                                                                       ggTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCa            <  4:114634/149‑1 (MQ=255)                                                                                                                                            tCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATc       >  3:290315/1‑149 (MQ=255)                                                                                                                                             ccTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGCAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATc       >  8:209220/1‑148 (MQ=255)                                                                                                                                               tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCa                 <  3:253525/136‑1 (MQ=255)                                                                                                                                               tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGCCAAGCTAAGTGTGCa                 >  4:253525/1‑136 (MQ=255)                                                                                                                                                 aaGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCACTGA‑TTAGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAAACTTCTTCGATCAGCGTCGCAAGCTAAGTGTGCAGTACAGCATCGAGcc  >  2:115025/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGA‑TTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAGCC  >  NC_002947/3006330‑3006622

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