Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,173,009	T→A	8.9%	Q165H (CAA→CAT) ‡	PP_1925 ←	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,173,009	0	T	A	8.9%	152.5 / 5.8	56	Q165H (CAA→CAT) ‡	PP_1925	monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base A (2/3);  total (26/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 = 7.99e-01

CCCGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAGC  >  NC_002947/2172861‑2173157                                                                                                                                                      |                                                                                                                                                      cccGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCt                                                                                                                                                       >  8:155204/1‑149 (MQ=255) cccGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCt                                                                                                                                                       <  5:112591/149‑1 (MQ=255)  ccGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCtt                                                                                                                                                      <  4:296476/149‑1 (MQ=255)    gTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTgg                                                                                                                                                    <  1:357227/149‑1 (MQ=255)        gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCTTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGaaa                                                                                                                                                 <  4:65573/148‑1 (MQ=255)        gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGaaa                                                                                                                                                 <  7:155204/148‑1 (MQ=255)           cAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCt                                                                                                                                             >  2:220532/1‑149 (MQ=255)           cAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                             <  5:28554/149‑3 (MQ=255)           cAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                             <  2:115745/149‑3 (MQ=255)              cacaCGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCgg                                                                                                                                          <  1:355999/149‑1 (MQ=255)              cacaCGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTGGGAAAGCCTCgg                                                                                                                                          >  8:240213/1‑149 (MQ=255)               acacGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCggg                                                                                                                                         >  2:300119/1‑149 (MQ=255)                   ggCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGtgat                                                                                                                                     <  1:295086/149‑1 (MQ=255)                      cATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGCTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCGTTCCTTGGAAAGCCTCGGGggattga                                                                                                                                  >  7:250934/1‑149 (MQ=255)                          aaCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAAt                                                                                                                              >  4:20588/1‑149 (MQ=255)                            cGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATgg                                                                                                                            <  4:323616/149‑1 (MQ=255)                              tcatcaGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGc                                                                                                                          <  8:34486/149‑1 (MQ=255)                                      cAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGTGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATTGGC‑GGAGTGaa                                                                                                                 >  7:252421/1‑149 (MQ=255)                                          atcGCCGGTTCCTCCGTGGTTACGCACGTTGCGTCGGCGACCTGGGTCAGCTCCGCATAGATTTGCGCCCGCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGTTAATGGGC‑GGAGTGAATctg                                                                                                             <  2:88037/147‑1 (MQ=255)                                            cGCCGGTTCGTCTCTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGAAAATGGGC‑GGAGTGAATctgct                                                                                                           <  1:321178/149‑1 (MQ=255)                                               cGGTTCCACCGTCGTAACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGAATTGCTCCCCCGAGTTGACAACACCAAC‑GACCAGCACGCCCTTTCCTTGGAAAGCCTCGGGAGATTGATAATGGGC‑GGAGTGAATCTGCTGTc                                                                                                        >  2:210183/1‑149 (MQ=255)                                                ggTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGGACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTcc                                                                                                       <  3:291178/149‑1 (MQ=255)                                                        cGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGACAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACa                                                                                               <  5:21054/149‑1 (MQ=255)                                                          ttggtACCCACGTTGCGTGGGTGACCTCGGAAAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGc                                                                                             <  8:252421/145‑1 (MQ=255)                                                                             ggCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCg                                                                          <  1:14714/149‑1 (MQ=255)                                                                                gACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGaa                                                                       <  4:162487/149‑1 (MQ=255)                                                                                   cTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACa                                                                     >  7:69875/1‑148 (MQ=255)                                                                                      ggACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATa                                                                   <  1:300119/147‑1 (MQ=255)                                                                                       gACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                                <  6:169379/149‑1 (MQ=255)                                                                                        aCAGCTCCACATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAgggg                                                               >  1:204532/1‑149 (MQ=255)                                                                                               cGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTc                                                         >  6:50435/1‑148 (MQ=255)                                                                                                gCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGcc                                                                           <  5:28421/129‑1 (MQ=255)                                                                                                gCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGcc                                                                           >  6:28421/1‑129 (MQ=255)                                                                                                  aTAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAgg                                                     >  4:203973/1‑149 (MQ=255)                                                                                                   tAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAAT‑GGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTc                                                   >  8:183037/1‑149 (MQ=255)                                                                                                    aGATTTGCGCCCCCGAGTTGCCACCACCAACGGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGATGTCCCTGGAACAGATCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGt                                                    >  1:231410/1‑149 (MQ=255)                                                                                                           cgcCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTg                                            <  1:220532/149‑1 (MQ=255)                                                                                                              ccccGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTgcgc                                         <  3:203973/149‑1 (MQ=255)                                                                                                              ccccGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTgcgc                                         >  2:346871/1‑149 (MQ=255)                                                                                                               cgcGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCt                                        >  6:29794/3‑149 (MQ=255)                                                                                                               ccagAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCGCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCt                                        <  1:150816/146‑1 (MQ=255)                                                                                                                        ccaccaccaAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATgacgc                                >  6:199760/1‑147 (MQ=255)                                                                                                                          accaccaAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCCGATGACGacac                             >  7:209566/1‑149 (MQ=255)                                                                                                                            caccaAC‑GACCAGCACGCGAGTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGa                           <  8:131596/149‑1 (MQ=255)                                                                                                                             accaAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAg                          <  1:346871/149‑1 (MQ=255)                                                                                                                                aaC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCt                       >  7:144605/1‑149 (MQ=255)                                                                                                                                       cAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGt                 >  7:146814/1‑149 (MQ=255)                                                                                                                                         gCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTc                >  1:78351/1‑149 (MQ=255)                                                                                                                                           aCGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTcgcg             >  2:314258/1‑149 (MQ=255)                                                                                                                                           aCGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATATGATTTCCCTGGAACAGATCTATCCCCGAATAGGCCGGAACATAGGGGTTGATCCAGGTCCCCGTTGCGATGATTACGACACGAGCCTCACAGTcgcg             >  6:369678/1‑149 (MQ=255)                                                                                                                                           aCGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAAT‑GGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            <  7:183037/149‑1 (MQ=255)                                                                                                                                            gaggcTTTCCAAGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            >  7:333416/4‑149 (MQ=255)                                                                                                                                            gaggcTTTCCAAGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            <  8:333416/146‑1 (MQ=255)                                                                                                                                            gaggcTTTCCAAGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            >  5:64090/4‑149 (MQ=255)                                                                                                                                            cgcgCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            <  5:29794/149‑1 (MQ=255)                                                                                                                                            cgcgCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            <  8:69875/149‑1 (MQ=255)                                                                                                                                               gcTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTcc         <  3:20588/149‑1 (MQ=255)                                                                                                                                                 tttCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGCCACGAGCCTCCCAGTCTCGGACCt        >  3:75879/1‑148 (MQ=255)                                                                                                                                                     cTTTGAAAGCCTCGGGTGATTGATAATGGGC‑TGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGa    <  5:369678/148‑1 (MQ=255)                                                                                                                                                      ttGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAGc  >  3:200309/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCCGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAAC‑GACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGC‑GGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAGC  >  NC_002947/2172861‑2173157

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