Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,971,955	C→G	9.8%	intergenic (‑390/+27)	PP_4379 ← / ← flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,971,955	0	C	G	9.8%	183.7 / 10.0	61	intergenic (‑390/+27)	PP_4379/flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base C (33/22);  new base G (3/3);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GGCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGG  >  NC_002947/4971808‑4972103                                                                                                                                                    |                                                                                                                                                      ggcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgca                                                                                                                                                      <  4:119483/149‑1 (MQ=255)    aCTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                                     >  5:133642/1‑147 (MQ=255)      tGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCagcagc                                                                                                                                                 >  8:286056/1‑149 (MQ=255)           cccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTgcgc                                                                                                                                                       >  1:58059/1‑138 (MQ=255)           cccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTgcgc                                                                                                                                                       <  2:58059/138‑1 (MQ=255)           cccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAg                                                                                                                                            <  2:267600/149‑1 (MQ=255)            ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCCGGAGAGAGCACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAGccctatc                                                                                                                                           <  1:236539/149‑5 (MQ=39)                aTGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       <  6:356263/149‑1 (MQ=255)                    tCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATcagg                                                                                                                                     >  3:210409/1‑147 (MQ=255)                      aGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATcaggcag                                                                                                                                  <  7:286056/148‑1 (MQ=255)                              agcagATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGAt                                                                                                                            >  6:309005/1‑146 (MQ=255)                                        gaaagggCCGCGATGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑g                                                                                                               >  5:157912/5‑148 (MQ=255)                                         aTGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑Gaa                                                                                                             >  4:3888/1‑149 (MQ=255)                                              gCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGc                                                                                                        >  1:72111/1‑149 (MQ=255)                                               ccGGATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGAt                                                                                                                                         >  5:80387/1‑116 (MQ=255)                                               ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGAt                                                                                                                                         <  6:80387/116‑1 (MQ=255)                                               ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCt                                                                                                       >  7:210647/1‑149 (MQ=255)                                               ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCt                                                                                                       <  8:210647/149‑1 (MQ=255)                                               ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCt                                                                                                       >  1:264874/1‑149 (MQ=255)                                                      gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTga                                                                                                <  4:210409/149‑1 (MQ=255)                                                       cgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       >  1:177279/1‑110 (MQ=255)                                                       cgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       <  2:177279/110‑1 (MQ=255)                                                                cgcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATgcgc                                                                                       <  6:157912/148‑1 (MQ=255)                                                                 gcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGAAGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGGAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATgcgct                                                                                      >  6:371783/1‑147 (MQ=255)                                                                   gACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCaa                                                                                   <  2:72111/149‑1 (MQ=255)                                                                      gCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATcag                                                                                                                                      >  6:373834/1‑96 (MQ=255)                                                                        aaGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                               >  3:280981/1‑148 (MQ=255)                                                                        aaGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                               >  4:358042/1‑148 (MQ=255)                                                                         aGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCa                                                                             <  6:133642/149‑1 (MQ=255)                                                                           gCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGCTCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGt                                                                           >  3:20647/1‑149 (MQ=255)                                                                            ccGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCGGAGAAATGCGCGCAAACGCCAGtt                                                                          >  3:154530/1‑149 (MQ=255)                                                                                tCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                                       >  5:24655/1‑148 (MQ=255)                                                                                        aTACGCCACCTCCAGGAGAGAGGACGCCCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTc                                                                >  5:125422/1‑147 (MQ=255)                                                                                        aTACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              <  2:264874/149‑1 (MQ=255)                                                                                            gCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCAt                                                          >  5:390771/1‑149 (MQ=255)                                                                                              cACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTA‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATgg                                                        >  3:59285/1‑149 (MQ=255)                                                                                                    cAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctgct                                                  <  3:3888/149‑1 (MQ=255)                                                                                                       gagagaGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctgctg                                                 >  2:113463/1‑147 (MQ=255)                                                                                                       gagagaGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCa                                               <  6:24655/149‑1 (MQ=255)                                                                                                        agagagGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAg                                              <  3:358042/149‑1 (MQ=255)                                                                                                                cgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCt                                      <  6:390771/149‑1 (MQ=255)                                                                                                                 gacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTg                                                                                                                   >  5:316896/1‑72 (MQ=255)                                                                                                                 gacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTg                                                                                                                   <  6:316896/72‑1 (MQ=255)                                                                                                                      gtgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAGCGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtgg                                >  7:200946/1‑149 (MQ=255)                                                                                                                       tgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtg                                 <  4:280981/147‑1 (MQ=255)                                                                                                                         ggggCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtgg                             >  4:239845/1‑149 (MQ=255)                                                                                                                          gggCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGc                            <  4:59285/149‑1 (MQ=255)                                                                                                                                   gtCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCAt                   >  5:190557/1‑149 (MQ=255)                                                                                                                                    tCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCtgctgc                                                >  8:126557/1‑119 (MQ=255)                                                                                                                                    tCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCtgctgc                                                <  7:126557/119‑1 (MQ=255)                                                                                                                                    tCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATg                  >  7:132247/1‑149 (MQ=255)                                                                                                                                      gcgATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGt                 >  4:234515/1‑148 (MQ=255)                                                                                                                                       tcctggggCTGCTGCGCAGCCCGAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              <  7:23668/97‑1 (MQ=255)                                                                                                                                       tcctggggCTGCTGCGCAGCCCGAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              >  8:23668/6‑102 (MQ=255)                                                                                                                                       tcctggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtg                              >  2:235692/6‑134 (MQ=255)                                                                                                                                       tcctggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtg                              <  1:235692/129‑1 (MQ=255)                                                                                                                                          tggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa            >  8:125438/3‑149 (MQ=255)                                                                                                                                          tAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa            >  6:36290/1‑149 (MQ=255)                                                                                                                                            gggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAGCGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTAt           <  8:200946/148‑1 (MQ=255)                                                                                                                                            gggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCCGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATc          <  4:154530/149‑1 (MQ=255)                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGt             <  6:336513/145‑1 (MQ=255)                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGt             >  5:336513/1‑145 (MQ=255)                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCg         >  2:152282/1‑149 (MQ=255)                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCg         >  3:352083/1‑149 (MQ=255)                                                                                                                                              gCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCg         >  1:79147/1‑149 (MQ=255)                                                                                                                                                    cagcagCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATgg  <  4:352083/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GGCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGG  >  NC_002947/4971808‑4972103

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