Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	812,048	T→C	6.2%	Q68Q (CAA→CAG)	PP_0698 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	812,048	0	T	C	6.2%	188.3 / 2.4	64	Q68Q (CAA→CAG)	PP_0698	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TGGCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTC  >  NC_002947/811903‑812195                                                                                                                                                  |                                                                                                                                                      tggCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCGAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTg                                                                                                                                                     >  4:1480/1‑148 (MQ=255)   gCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAg                                                                                                                                                  >  3:151340/1‑149 (MQ=255)    cATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGATCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGGAGATCGCCTGGGTGGTGCCTTGCAgg                                                                                                                                                 <  4:258188/149‑1 (MQ=255)    cATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAgg                                                                                                                                                 <  3:220505/149‑1 (MQ=255)     aTGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAggg                                                                                                                                                <  2:366965/149‑1 (MQ=255)     aTGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCAAGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAggg                                                                                                                                                >  1:266961/1‑149 (MQ=255)      tGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGc                                                                                                                                               <  3:156321/149‑1 (MQ=255)           gcgcGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTGCAAGGcaccac                                                                                                                                          <  1:204705/149‑1 (MQ=255)            cgcgGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGAcc                                                                                                                                         <  6:171490/149‑1 (MQ=255)             gcgGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCAGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACc‑‑g                                                                                                                                      >  6:183444/1‑149 (MQ=255)              cgGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccaccca                                                                                                                                       <  7:99623/149‑3 (MQ=255)                gCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACc‑‑ggcg                                                                                                                                   >  1:189873/1‑149 (MQ=255)                   ggCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCGAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCgggg                                                                                                                                <  3:1480/149‑1 (MQ=255)                    gcgcGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCggggg                                                                                                                               >  6:319694/1‑149 (MQ=255)                        gggATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCGCAGCC‑‑AGCGGGGGTGAg                                                                                                                           >  1:344629/1‑149 (MQ=255)                        gggATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAg                                                                                                                           >  6:350930/1‑149 (MQ=255)                        gggATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAg                                                                                                                           <  5:350930/149‑1 (MQ=255)                                ccagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGCGGAAAATgg                                                                                                                   >  3:101033/1‑149 (MQ=255)                                 cagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGaa                                                                                                                        >  6:59216/1‑143 (MQ=255)                                 cagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGaa                                                                                                                        <  5:59216/143‑1 (MQ=255)                                    cgcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCg                                                                                                               <  7:356810/149‑1 (MQ=255)                                       aTGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGtt                                                                                                            >  6:94622/1‑149 (MQ=255)                                          gCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGGGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTgcg                                                                                                         <  5:89022/149‑1 (MQ=255)                                            gATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTgcggc                                                                                                       >  1:362220/1‑149 (MQ=255)                                                   ggTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGa                                                                                                 <  3:326067/148‑1 (MQ=255)                                                    gTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                >  1:146270/1‑148 (MQ=255)                                                        gCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGt                                                                                            <  5:319694/148‑1 (MQ=255)                                                         cAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTc                                                                                          >  3:200862/1‑149 (MQ=255)                                                                        acaTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttc                                                                           >  8:216942/1‑149 (MQ=255)                                                                          aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttctt                                                                         >  2:130646/1‑149 (MQ=255)                                                                           ttGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                        >  2:28894/1‑149 (MQ=255)                                                                           ttGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                        >  4:289500/1‑149 (MQ=255)                                                                            tGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTg                                                                                                           >  6:392065/1‑113 (MQ=255)                                                                            tGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTg                                                                                                           <  5:392065/113‑1 (MQ=255)                                                                                gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAg                                                                                 <  3:292140/135‑1 (MQ=255)                                                                                gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAg                                                                                 >  4:292140/1‑135 (MQ=255)                                                                                             tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcg                                                      >  5:393107/1‑149 (MQ=255)                                                                                               cAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcggc                                                    >  8:51428/1‑149 (MQ=255)                                                                                               cAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcggc                                                    >  1:96606/1‑149 (MQ=255)                                                                                               cAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcggc                                                    >  4:41565/1‑149 (MQ=255)                                                                                               cAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                     <  4:151340/148‑1 (MQ=255)                                                                                               cAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcg                                                      <  7:216942/147‑1 (MQ=255)                                                                                                   cgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCGGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGAc                                                >  4:189399/1‑149 (MQ=255)                                                                                                      tGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACcgc                                             >  8:339743/1‑149 (MQ=255)                                                                                                        cACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCgcgc                                           <  2:146270/149‑1 (MQ=255)                                                                                                          cGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCg                                         <  2:189873/149‑1 (MQ=255)                                                                                                           ggCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgt                                        <  5:183444/149‑1 (MQ=255)                                                                                                           ggCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgt                                        <  2:96606/149‑1 (MQ=255)                                                                                                             cgcCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                      <  1:28894/149‑1 (MQ=255)                                                                                                                   cGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGAc                                 <  2:362220/148‑1 (MQ=255)                                                                                                                    gAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGt                               <  5:94622/149‑1 (MQ=255)                                                                                                                      gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGc                             <  2:266961/149‑1 (MQ=255)                                                                                                                          ttCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGc                                                                <  4:69783/110‑1 (MQ=255)                                                                                                                          ttCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGc                                                                >  3:69783/1‑110 (MQ=255)                                                                                                                                  ccccgccggtcgtGCCCTGCAAGGCACCACCCAGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTCACGTGCAAGTGCTcgg                   >  4:265837/8‑149 (MQ=38)                                                                                                                                     ccTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                     <  2:223404/142‑1 (MQ=255)                                                                                                                                     ccTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                     >  1:223404/1‑142 (MQ=255)                                                                                                                                          gtggtgCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCa          >  4:310342/1‑148 (MQ=255)                                                                                                                                          gtggtgCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCa          <  3:310342/148‑1 (MQ=255)                                                                                                                                          gtcgtgCCCTGCAAGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCAc         >  7:25897/4‑149 (MQ=255)                                                                                                                                              tgCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAg     <  7:339743/149‑1 (MQ=255)                                                                                                                                              tgCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGGATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGATCGGCGGCGCGCACGAAg     >  8:141198/1‑149 (MQ=255)                                                                                                                                                 cTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGctc  <  7:51428/149‑1 (MQ=255)                                                                                                                                                 cTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGAAGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCCCACGAAGct   >  1:134842/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                      TGGCATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTC  >  NC_002947/811903‑812195

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