Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,418,818	G→T	13.1%	P85T (CCT→ACT)	PP_4758 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,418,818	0	G	T	13.1%	171.1 / 15.7	61	P85T (CCT→ACT)	PP_4758	MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base T (4/4);  total (31/30)
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.69e-01

CTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACATGAT  >  NC_002947/5418671‑5418959                                                                                                                                                    |                                                                                                                                               cTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGt                                                                                                                                               <  8:159175/149‑1 (MQ=255)    cGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑Ac                                                                                                                                            <  8:85599/148‑1 (MQ=255)    cGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑Ac                                                                                                                                            <  8:269809/148‑1 (MQ=255)       ccACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCg                                                                                                                                        <  4:132639/149‑1 (MQ=255)         aCCCCGAGCAGAAAACGCAGCACGTAGAACTGGGGTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGGGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTA‑CCTTCGaa                                                                                                                                      <  2:105381/149‑1 (MQ=255)           cccGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACa                                                                                                                                    <  4:580/149‑1 (MQ=255)            ccgcccAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCCTGGCCGTGGCCACCATCCCCCAGGTGATCATGATACGGGTCAGCCACACCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTA‑CCTTCGAACat                                                                                                                                   <  8:117355/144‑2 (MQ=255)                           agcaCGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATgaa                                                                                                                    <  8:171920/149‑1 (MQ=255)                                   gAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGc                                                                                                             <  1:19807/148‑1 (MQ=255)                                          gTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCAc                                                                                                      >  7:45546/1‑148 (MQ=255)                                           tttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCACGGCAcc                                                                                                     >  7:117481/1‑148 (MQ=255)                                           tttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGAGGAAGAACAGCCCGGCAcc                                                                                                     >  7:140423/1‑148 (MQ=255)                                                ttCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTA‑CCTTCGAACat                                                                                                                                   >  3:293485/1‑112 (MQ=255)                                                ttCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTA‑CCTTCGAACat                                                                                                                                   <  4:293485/113‑2 (MQ=255)                                                     cacaAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACAGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCAt                                                                                             <  1:68749/146‑1 (MQ=255)                                                      acaAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGc                                                                                         >  2:157198/1‑149 (MQ=255)                                                        aaaGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCAACGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCg                                                                                       >  8:21700/1‑149 (MQ=255)                                                          aGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCa                                                                                     <  1:157198/149‑1 (MQ=255)                                                          aGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAACCCATAGGCCGc                                                                                      >  4:8607/1‑148 (MQ=255)                                                             ccATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAgcg                                                                                  >  5:113932/1‑149 (MQ=255)                                                              cATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACAGATGAAGAACAGCCCGGCCCCGAAGCCATCGGCCGCAgcgc                                                                                 >  3:197889/1‑149 (MQ=255)                                                                      tGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGccc                                                                         <  6:113932/149‑1 (MQ=255)                                                                             caGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCTGCAGCGCTGATGCCCAGGt                                                                     >  4:239816/1‑146 (MQ=255)                                                                             caGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCTGCAGCGCTGATGCCCAGGt                                                                     <  3:239816/146‑1 (MQ=255)                                                                                   ccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACGGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTc                                                              >  5:276227/1‑147 (MQ=255)                                                                                       gTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTg                                                        <  7:129266/149‑1 (MQ=255)                                                                                       gTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCAAAGGT‑ACATCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTAGCTTTCCACGTg                                                        >  1:252167/1‑149 (MQ=255)                                                                                        tGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgc                                                       <  3:45587/149‑1 (MQ=255)                                                                                             aTGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcg                                                  <  8:45546/149‑1 (MQ=255)                                                                                               gATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGa                                                >  3:252978/1‑149 (MQ=255)                                                                                                aTACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGa                                                <  7:225131/148‑1 (MQ=255)                                                                                                aTACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGAc                                               >  5:39336/1‑149 (MQ=255)                                                                                                aTACGGGGCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGa                                                <  3:77304/148‑1 (MQ=255)                                                                                                    gggTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAg                                           >  6:196159/1‑149 (MQ=255)                                                                                                      gTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGc                                          >  2:141494/1‑148 (MQ=255)                                                                                                        cAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑AATTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCAAGGTGCGGGCGGACGAAGCCaa                                       >  4:84584/1‑149 (MQ=255)                                                                                                            cacaGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGtt                                   <  1:141494/149‑1 (MQ=255)                                                                                                             acaGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCCCCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGtt                                   >  7:296556/1‑148 (MQ=255)                                                                                                               aGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGAc                                >  3:272943/1‑149 (MQ=255)                                                                                                                   gggCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTcc                                                             >  4:98472/1‑116 (MQ=255)                                                                                                                   gggCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTcc                                                             <  3:98472/116‑1 (MQ=255)                                                                                                                    ggCACCGACCTTTTGCAGCAGCATGTTCGAAGGTAACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCgg                            >  1:212688/1‑149 (MQ=255)                                                                                                                    ggCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCa                                                            >  4:228735/1‑116 (MQ=255)                                                                                                                    ggCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCa                                                            <  3:228735/116‑1 (MQ=255)                                                                                                                        ccGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAg                       <  4:252978/149‑1 (MQ=255)                                                                                                                         cGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAgg                      >  8:218557/1‑149 (MQ=255)                                                                                                                         cGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCAATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAgg                      >  5:164131/1‑149 (MQ=255)                                                                                                                          gACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGt                     <  1:128473/149‑1 (MQ=255)                                                                                                                                ttGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGtt                 <  3:8607/147‑1 (MQ=255)                                                                                                                                 tGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGa               <  2:211620/148‑1 (MQ=255)                                                                                                                                  gcagcagcaTGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGAc              >  3:257255/1‑148 (MQ=255)                                                                                                                                      cagcaTGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGt                     >  4:91871/1‑137 (MQ=255)                                                                                                                                      cagcaTGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGt                     <  3:91871/137‑1 (MQ=255)                                                                                                                                       agcaTGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGaa        <  6:276227/149‑1 (MQ=255)                                                                                                                                           tGTTCGAAGTA‑CCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTa                    >  7:285226/1‑133 (MQ=255)                                                                                                                                           tGTTCGAAGTA‑CCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTa                    <  8:285226/133‑1 (MQ=255)                                                                                                                                           tGTTCGAAGTA‑CCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATTAACat     >  5:253538/1‑148 (MQ=255)                                                                                                                                             ttCGAAGTA‑CCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatgat  >  2:229075/1‑149 (MQ=255)                                                                                                                                             ttCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatgat  <  7:37433/149‑1 (MQ=255)                                                                                                                                             ttCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatgat  >  8:37433/1‑149 (MQ=255)                                                                                                                                             ttCGAAGGT‑AATTCGAACAGCGCGTAACCGATGAAGAACAGCCCGCAACCGAAGCCATAGGCCGCAGCCCTAATGCCCAGGTCGCTTTCCAGGTGCGGCCGGACGAAGCCAATGTTGACGCGGTGGAGGAAGTTGACGATGAACatgat  >  7:184144/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               CTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGT‑ACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACATGAT  >  NC_002947/5418671‑5418959

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