Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,587	T→C	11.3%	L31P (CTG→CCG) ‡	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,587	0	T	C	11.3%	161.7 / 10.2	62	L31P (CTG→CCG) ‡	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (4/3);  total (33/29)
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.55e-01

GGGCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAGGC  >  NC_002947/5187439‑5187735                                                                                                                                                     |                                                                                                                                                     ggtcATGGCTGCTGAATGAGCGGTTGTCGATGGCGCATCTATACGGCGGACTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATCAACCGCGCGTACAAAAAGGTCCTGTTCCGCCt                                                                                                                                                      <  6:323606/146‑1 (MQ=255)       ggCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGaa                                                                                                                                                  >  4:359294/1‑147 (MQ=255)         ctgctgAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAgg                                                                                                                                              <  1:71651/149‑1 (MQ=255)          tgctgAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGa                                                                                                                                             <  4:362960/149‑1 (MQ=255)             tgAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACg                                                                                                                                           >  3:239377/1‑148 (MQ=255)              gAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGAc                                                                                                                                            >  7:68748/1‑146 (MQ=255)                 tGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGACCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGAccttttt                                                                                                                                      <  7:225813/149‑7 (MQ=255)                   aGCCGTTGTCGATGGCTCATCTATACGGCGGCCTTCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATTAACACCGCGTACAAAAAGGTACTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                     <  7:43934/148‑1 (MQ=255)                    gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                     <  4:36018/147‑1 (MQ=255)                    gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                     >  3:36018/1‑147 (MQ=255)                    gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGa                                                                                                                                             <  2:341867/139‑1 (MQ=255)                    gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGa                                                                                                                                             >  1:341867/1‑139 (MQ=255)                      cGTTGTCGATGGCCCATCTATACGGCGGCCTCCTGATTGGACTGGCCGTGTGGCTGGTGCTGCGCCCGGCGCGAAGTAGAAGCGGGGAACTGCATGAACAGCGCGTACAAAACGGTCCTGTTCCGCCTGGGACAGGCCTAGAACGGCTa                                                                                                                                 >  4:333760/1‑149 (MQ=255)                        ttGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTAcc                                                                                                                               >  2:33905/1‑149 (MQ=255)                              aTGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCgg                                                                                                                         <  8:238327/149‑1 (MQ=255)                                 gcgcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCt                                                                                                                      >  2:335882/1‑149 (MQ=255)                                   gcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCg                                                                                                                    >  8:20439/1‑149 (MQ=255)                                         ataCGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAg                                                                                                              >  5:9885/1‑149 (MQ=255)                                         ataCGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAg                                                                                                              >  3:79744/1‑149 (MQ=255)                                          taCGGCGGCCTGCTGATTGGCCTGGCCGTGTGGCGGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAg                                                                                                              <  4:119392/148‑1 (MQ=255)                                             ggcggcCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGcc                                                                                                          >  6:316302/1‑149 (MQ=255)                                                 gcCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCtgtg                                                                                                      <  6:363015/149‑1 (MQ=255)                                                       tgATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGcc                                                                                                  >  5:259291/1‑147 (MQ=255)                                                         aTTGGCCTGGCCTTGGGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAg                                                                                               <  3:136391/148‑1 (MQ=255)                                                              cctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTg                                                                                         <  6:259291/149‑1 (MQ=255)                                                               ctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTg                                                                                         >  4:45770/1‑148 (MQ=255)                                                                    ctTGTGGCTGGTGCTGCGGCCGGGGCGGAGTTGATGCGGGGTACTGCATGAAGAGGGCGTAAAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAAGGGGTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAggc                                                                                   <  4:239377/149‑1 (MQ=255)                                                                       gtgGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAggcgg                                                                                 <  7:175851/148‑1 (MQ=255)                                                                          gCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTAt                                                                             <  8:68748/149‑1 (MQ=255)                                                                            ttcaactgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGc                                                                           >  5:351803/6‑149 (MQ=255)                                                                            ttcaactgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGc                                                                           >  6:80257/6‑149 (MQ=255)                                                                                actgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATc                                                                       >  4:19002/2‑149 (MQ=255)                                                                                actgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATc                                                                       >  1:365206/2‑149 (MQ=255)                                                                                    cgcCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGac                                                                    <  3:359294/148‑1 (MQ=255)                                                                                                               cTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATc                                         <  1:33905/148‑1 (MQ=255)                                                                                                                  cATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGccc                                     >  5:112478/1‑149 (MQ=255)                                                                                                                     gAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                  >  7:335666/1‑149 (MQ=255)                                                                                                                      aaCAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCggg                                 <  6:351803/149‑1 (MQ=255)                                                                                                                      aaCAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                  <  3:19002/148‑1 (MQ=255)                                                                                                                        cAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGCCATCGCCCcgcgcg                               >  6:117093/1‑149 (MQ=255)                                                                                                                           cgcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGAc                             <  4:295844/148‑1 (MQ=255)                                                                                                                           cgcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACg                            >  2:222604/1‑149 (MQ=255)                                                                                                                            gcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGa                                                                     <  1:62831/107‑1 (MQ=255)                                                                                                                            gcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGa                                                                     >  2:62831/1‑107 (MQ=255)                                                                                                                              gcaCAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGa                         >  3:401765/3‑149 (MQ=255)                                                                                                                               tACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGTCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGAt                        <  1:335882/149‑1 (MQ=255)                                                                                                                                  aaaaaGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGa                     >  3:88747/1‑149 (MQ=255)                                                                                                                                    aaaGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTGGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGaca                   <  3:126436/149‑1 (MQ=255)                                                                                                                                        gTCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcg               >  4:4269/1‑149 (MQ=255)                                                                                                                                        gTCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcg               >  2:77626/1‑149 (MQ=255)                                                                                                                                         tCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcgc              >  4:113915/1‑149 (MQ=255)                                                                                                                                         tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                  >  1:103994/1‑129 (MQ=255)                                                                                                                                         tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                  <  2:103994/129‑1 (MQ=255)                                                                                                                                           cTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGa            <  2:365206/149‑1 (MQ=255)                                                                                                                                             gTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAg           >  8:110996/1‑148 (MQ=255)                                                                                                                                                ccGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAg       >  8:112361/1‑149 (MQ=255)                                                                                                                                                 cGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAg       <  3:45770/148‑1 (MQ=255)                                                                                                                                                   ccTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAg    <  4:401765/149‑1 (MQ=255)                                                                                                                                                    cTGGGACAGGTCGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAgg   <  3:333760/149‑1 (MQ=255)                                                                                                                                                    cTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATg                                              <  3:153001/106‑1 (MQ=255)                                                                                                                                                    cTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATg                                              >  4:153001/1‑106 (MQ=255)                                                                                                                                                    cTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                  <  4:46188/118‑1 (MQ=255)                                                                                                                                                    cTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                  >  3:46188/1‑118 (MQ=255)                                                                                                                                                     tGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGTCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAGGc  >  2:167644/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGGCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAGGC  >  NC_002947/5187439‑5187735

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