Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,950,972	T→G	18.8%	intergenic (‑122/+35)	fliL ← / ← fliK	flagellar protein FliL/flagellar hook‑length control protein FliK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,950,972	0	T	G	18.8%	90.9 / 21.5	53	intergenic (‑122/+35)	fliL/fliK	flagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base G (5/5);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

CTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAC  >  NC_002947/4950832‑4951114                                                                                                                                                      |                                                                                                                                               cTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCcgc                                                                                                                                         <  8:185239/148‑1 (MQ=255)    cTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGGGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCcgcg                                                                                                                                        >  1:51384/1‑146 (MQ=255)    cTTCGCTCTTCGCCAAGCCAATAATCCGTCGTCCATCGGGGGTTCAAATGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTGCTT‑‑‑‑‑CGCGGGCTCGCCCGCgaag                                                                                                                                     <  7:303302/149‑1 (MQ=255)     ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCgaag                                                                                                                                     >  6:227378/1‑148 (MQ=255)     ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCgaag                                                                                                                                     >  6:156275/1‑148 (MQ=255)     ttCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCgaag                                                                                                                                     >  1:140370/1‑148 (MQ=255)       cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCCGCgaagaa                                                                                                                                   <  8:64521/148‑1 (MQ=255)       atctctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCCGCGaagaag                                                                                                                                  <  3:298672/147‑1 (MQ=255)        gCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGaagaag                                                                                                                                  <  3:246713/148‑1 (MQ=255)        gCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAagc                                                                                                                                 <  8:156398/149‑1 (MQ=255)          tctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAagc                                                                                                                                 >  3:109455/1‑147 (MQ=255)            ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGTATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAagcca                                                                                                                               <  6:43581/147‑1 (MQ=255)            ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAagccagc                                                                                                                             >  4:268574/1‑149 (MQ=255)            ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAagcca                                                                                                                               <  2:140370/147‑1 (MQ=255)            ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCCGCGaagaag                                                                                                                                  >  4:298672/1‑144 (MQ=255)                ccATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAg                                                                                                                         <  5:227378/149‑1 (MQ=255)                ccATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAg                                                                                                                         >  6:236532/1‑149 (MQ=255)                  aTGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCCGCgaa                                                                                                                                      >  3:313197/1‑134 (MQ=255)                  aTGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCGAGCCCGCgaa                                                                                                                                      <  4:313197/134‑1 (MQ=255)                      cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                    >  5:369271/1‑148 (MQ=255)                       aatactCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACa                                                                                                                  <  7:33444/149‑1 (MQ=255)                            tCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑AGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTAt                                                                                                             >  2:17188/1‑149 (MQ=255)                                   tcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgt                                                                                                      <  5:156275/149‑1 (MQ=255)                                     cATCGGGGTTTCAATGCATAATGGGTTGGCTTGAGCATGTGTTATGCAAGATTTTGCTGGGTATGTATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgtag                                                                                                    <  1:180761/149‑1 (MQ=255)                                     cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                     <  5:116339/148‑1 (MQ=255)                                     cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                     <  6:369271/148‑1 (MQ=255)                                                   aGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGTTTGGGGGCTTCTTGCTAGCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGAGTAGTAATc                                                                                               >  1:300071/1‑149 (MQ=255)                                                                 gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAac                                                                          >  6:68132/1‑147 (MQ=255)                                                                 gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAac                                                                          <  5:68132/147‑1 (MQ=255)                                                                  cttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAaccac                                                                       <  5:236532/149‑1 (MQ=255)                                                                   ttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGCCCGCACTGTCGCCGCCCaccacc                                                                      >  8:202891/1‑149 (MQ=255)                                                                     gAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctg                                                                    >  8:43798/1‑149 (MQ=255)                                                                                  gCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGc                                                       >  4:276546/1‑149 (MQ=255)                                                                                   ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                      >  3:159817/1‑149 (MQ=255)                                                                                   ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑AGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGc                                                       <  1:17188/148‑1 (MQ=255)                                                                                             catgggATGTATTCCATTTTGTTTTGGGTTGGT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCGCAGGCCAACAGTTATCATGAGTAGTATTCGAGCGCTCTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCTGGCGCTATc                                            <  4:268857/146‑1 (MQ=255)                                                                                             cAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                            <  3:268574/149‑1 (MQ=255)                                                                                             cAGCGGAGGGATTCCTGCTGGTTTTGCGTTTCT‑‑‑‑GGCTGATT‑‑‑‑‑CGGGGGCTCCCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCACGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTAt                                             <  8:253590/148‑1 (MQ=255)                                                                                                 ggatggatTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAc                                         <  1:37343/148‑1 (MQ=255)                                                                                                 ggatggatTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAc                                         >  2:37343/1‑148 (MQ=255)                                                                                                       atTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCg                                  <  2:300071/149‑1 (MQ=255)                                                                                                           cagcagGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATAGACCGCACTGTCGCCCACAACCACCAGCTTCTACAAAGGCCGGGCGCTATAGACCGCTTCCGccc                              >  7:278605/1‑149 (MQ=255)                                                                                                                     tGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCt                     <  3:276546/148‑1 (MQ=255)                                                                                                                       cGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCg                   <  7:202891/148‑1 (MQ=255)                                                                                                                                  ggCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGc                                >  4:71293/1‑128 (MQ=255)                                                                                                                                  ggCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGc                                <  3:71293/128‑1 (MQ=255)                                                                                                                                   gCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgc          >  3:71418/1‑149 (MQ=255)                                                                                                                                        tt‑‑‑‑‑cGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCCCTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTACTCGGCCCGGCGcgccgc     >  1:31846/1‑149 (MQ=255)                                                                                                                                        tt‑‑‑‑‑cGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccgc     >  1:231466/1‑149 (MQ=255)                                                                                                                                        tt‑‑‑‑‑cGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc       <  4:59068/147‑1 (MQ=255)                                                                                                                                        tt‑‑‑‑‑cGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc       >  3:59068/1‑147 (MQ=255)                                                                                                                                         t‑‑‑‑‑cGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc    >  3:17962/1‑149 (MQ=255)                                                                                                                                                gcgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCCCCCTGCTCGGCCCGGCGCGCCGCCAc  >  4:46431/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                               CTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCT‑‑‑‑GGCTTCTT‑‑‑‑‑CGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAC  >  NC_002947/4950832‑4951114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑