Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,971,956	A→C	10.0%	intergenic (‑391/+26)	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,956	0	A	C	10.0%	174.6 / 10.2	60	intergenic (‑391/+26)	PP_4379/flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base A (32/22);  new base C (3/3);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

GGCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGGC  >  NC_002947/4971808‑4972104                                                                                                                                                     |                                                                                                                                                      ggcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgca                                                                                                                                                       <  4:119483/149‑1 (MQ=255)    aCTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                                      >  5:133642/1‑147 (MQ=255)      tGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCagcagc                                                                                                                                                  >  8:286056/1‑149 (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                                                                                                        >  1:264874/1‑149 (MQ=255)                                               ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCt                                                                                                        <  8:210647/149‑1 (MQ=255)                                               ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCt                                                                                                        >  7:210647/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                                                                                >  4:358042/1‑148 (MQ=255)                                                                        aaGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                                >  3:280981/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                               <  1:235692/129‑1 (MQ=255)                                                                                                                                       tcctggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtg                               >  2:235692/6‑134 (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          >  3:352083/1‑149 (MQ=255)                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCg          >  2:152282/1‑149 (MQ=255)                                                                                                                                              gCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCg          >  1:79147/1‑149 (MQ=255)                                                                                                                                                    cagcagCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATgg   <  4:352083/149‑1 (MQ=255)                                                                                                                                                     agcagcCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGGc  <  6:190557/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTC‑GAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGGC  >  NC_002947/4971808‑4972104

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