Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,816,642	G→T	6.2%	T165N (ACT→AAT)	chbG ←	chito‑oligosaccharide deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,816,642	0	G	T	6.2%	214.2 / 3.6	64	T165N (ACT→AAT)	chbG	chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (37/23);  new base T (2/2);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

GGATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGCGTC  >  NC_000913/1816502‑1816789                                                                                                                                             |                                                                                                                                                    ggATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAgg                                                                                                                                             <  8:337351/149‑1 (MQ=255)  gATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTTCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGt                                                                                                                                            <  6:311736/149‑1 (MQ=255)  gATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGt                                                                                                                                            <  3:71959/149‑1 (MQ=255)         ttatCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCgg                                                                                                                                     >  8:236992/1‑149 (MQ=255)           atCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGc                                                                                                                                    >  6:371644/1‑148 (MQ=255)                        ggATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCAt                                                                                                                      <  2:237575/149‑1 (MQ=255)                           tGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGCCCGGCAAATCAAACGCCATCt                                                                                                                    >  1:17609/1‑148 (MQ=255)                           tGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCt                                                                                                                    >  2:340448/1‑148 (MQ=255)                             acacATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGa                                                                                                                  <  1:340448/148‑1 (MQ=255)                               acaTCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGa                                                                                                               <  7:279732/149‑1 (MQ=255)                                  tCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATcc                                                                                                            >  1:85417/1‑149 (MQ=255)                                      cTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATccgccc                                                                                                        >  4:227761/1‑149 (MQ=255)                                      cTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATccgccc                                                                                                        >  1:135343/1‑149 (MQ=255)                                       tCCAGTGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGGCCGGCAAATCAAACGCCATCTGACGATccgcccg                                                                                                       >  1:306441/1‑149 (MQ=255)                                         cAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATccgcccg                                                                                                       <  3:227761/147‑1 (MQ=255)                                               aCGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCcgcaacgcaa                                                                                               <  7:236992/149‑1 (MQ=255)                                               aCGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCcgcaacgcaa                                                                                               <  5:371644/149‑1 (MQ=255)                                                     aCCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATACACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCt                                                                                         >  7:297235/1‑149 (MQ=255)                                                          gATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGctct                                                                                    >  4:363461/1‑149 (MQ=255)                                                          gATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGctct                                                                                    >  3:286219/1‑149 (MQ=255)                                                                cTGCGTCATCCAGCACTTGCAGGCACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGATGAATCCCTGGGTAGTGCGCAGTTTTACCGGCAAATCAAACGCCATCTGACGCTCCGCCCGCAACGCAATGCCCTGCTCTGCCGc                                                                               >  7:287427/1‑148 (MQ=255)                                                                cTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGTTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGc                                                                               >  4:120589/1‑148 (MQ=255)                                                                        tCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  2:351473/1‑148 (MQ=255)                                                                        tCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  3:269690/1‑148 (MQ=255)                                                                         ccAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGc                                                                      >  8:56642/1‑148 (MQ=255)                                                                                ttGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCg                                                              >  6:138295/1‑149 (MQ=255)                                                                                ttGCAGGAACAGCGACTCACAGATCTCTTCACCGTAGACTGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTCGCAACAATCg                                                              >  6:248796/1‑149 (MQ=255)                                                                                   cAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCggg                                                            <  7:56642/148‑1 (MQ=255)                                                                                            cGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCAGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCg                                                  >  4:120891/1‑149 (MQ=255)                                                                                                tCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGaa                                              >  2:234943/1‑149 (MQ=255)                                                                                                  aCTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCCACGCAATGCCCTGCTCTGCCCCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACa                                            >  4:279581/1‑149 (MQ=255)                                                                                                  aCTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACa                                            <  3:363461/149‑1 (MQ=255)                                                                                                   cTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACat                                           <  2:85417/149‑1 (MQ=255)                                                                                                    tGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACata                                          >  1:301071/1‑149 (MQ=255)                                                                                                    tGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACata                                          >  2:259843/1‑149 (MQ=255)                                                                                                     gATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACatat                                         >  3:122293/1‑149 (MQ=255)                                                                                                      aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       <  3:78551/118‑1 (MQ=255)                                                                                                      aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  4:78551/1‑118 (MQ=255)                                                                                                      aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       >  4:83850/1‑118 (MQ=255)                                                                                                      aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTATTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTgg                                                                       <  3:83850/118‑1 (MQ=255)                                                                                                      aTCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATg                                        >  7:133624/1‑149 (MQ=255)                                                                                                             cACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatga                                 >  2:230412/1‑149 (MQ=255)                                                                                                              aCCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGGCCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatga                                 <  2:306441/148‑1 (MQ=255)                                                                                                                cGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACattatga                              <  4:122293/149‑1 (MQ=255)                                                                                                                cGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatga                              <  5:248796/149‑1 (MQ=255)                                                                                                                 gTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatgat                             <  4:269690/149‑1 (MQ=255)                                                                                                                 gTAGAATGCACTGCTGAATCACTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatgat                             <  3:6316/149‑1 (MQ=255)                                                                                                                  tAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACatgatgatg                            <  1:351473/149‑1 (MQ=255)                                                                                                                     aaTGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGc                          >  6:166421/1‑148 (MQ=255)                                                                                                                       tGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTAt                       <  8:287427/149‑1 (MQ=255)                                                                                                                          aCTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCCATGCCCTGCTCTGCCGCCAAACTGGCAACAATCGGGATAATCTGCGGGAACATTTGCACATGATGATGGCTATc                      >  4:122303/1‑147 (MQ=255)                                                                                                                          aCTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATc                      >  7:308741/1‑147 (MQ=255)                                                                                                                              ctgAATCCCTGGGGAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGCTCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATg                <  7:305766/149‑1 (MQ=255)                                                                                                                               tgAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               >  1:332205/1‑149 (MQ=255)                                                                                                                               tgAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               >  5:292880/1‑149 (MQ=255)                                                                                                                               tgAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               >  7:281653/1‑149 (MQ=255)                                                                                                                               tgAATCCCTGGGTAGTGCGCAGGTTCTCCGGCAAATCAACCGCCATCTGACCATCCGCCCGCAACGCAAAGCCCTGCACTGACGCAAACCTGCCAACAAGCGGGAAAAACTGCGGGACCATATGCACCTGATGATGGCTATCAACAt                 >  3:35305/1‑147 (MQ=255)                                                                                                                                 aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGTTCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGt             <  2:135343/149‑1 (MQ=255)                                                                                                                                 aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGc               <  4:286219/147‑1 (MQ=255)                                                                                                                                 aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCg              >  3:278978/1‑148 (MQ=255)                                                                                                                                 aaTCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCg              >  8:377873/1‑148 (MQ=255)                                                                                                                                      cTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGtt        >  3:257127/1‑149 (MQ=255)                                                                                                                                         ggTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGc     <  5:166421/149‑1 (MQ=255)                                                                                                                                            aGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGACCATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGCGTc  >  4:390720/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                    GGATAGTATTATCGATAAACGCCGGATGACACATCACCTCCAGCGAACGATCACCCCGATGGCCTGCGTCATCCAGCACTTGCAGGAACAGCGACTCACTGATCTCTTCACCGTAGAATGCACTGCTGAATCCCTGGGTAGTGCGCAGGTTGACCGGCAAATCAAACGCCATCTGACGATCCGCCCGCAACGCAATGCCCTGCTCTGCCGCAAACCTGGCAACAATCGGGAAAATCTGCGGGAACATATGCACATGATGATGGCTATCAAGATGCGTAGGTTTGCGTC  >  NC_000913/1816502‑1816789

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