Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,586	C→G	12.3%	L31V (CTG→GTG) ‡	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,586	0	C	G	12.3%	149.6 / 12.4	57	L31V (CTG→GTG) ‡	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base G (4/3);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

GGGGCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAGG  >  NC_002947/5187438‑5187734                                                                                                                                                     |                                                                                                                                                     ggggCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGcc                                                                                                                                                      <  3:360201/149‑1 (MQ=255)  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                                                                                                                                            >  1:341867/1‑139 (MQ=255)                     gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGa                                                                                                                                            <  2:341867/139‑1 (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              >  2:77626/1‑149 (MQ=255)                                                                                                                                         gTCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcg              >  4:4269/1‑149 (MQ=255)                                                                                                                                          tCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcgc             >  4:113915/1‑149 (MQ=255)                                                                                                                                          tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                 <  2:103994/129‑1 (MQ=255)                                                                                                                                          tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                 >  1:103994/1‑129 (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)                                                                                                                                                     |                                                                                                                                                     GGGGCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAGG  >  NC_002947/5187438‑5187734

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