Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,400,124	A→G	5.8%	S32G (AGC→GGC)	ybgF →	TolQRA transport system periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,400,124	0	A	G	5.8%	205.5 / 2.1	69	S32G (AGC→GGC)	ybgF	TolQRA transport system periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base G (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCAAGGCATCATCGA  >  NC_002947/1399977‑1400263                                                                                                                                                    |                                                                                                                                              gacgaGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACgcaggcag                                                                                                                                              <  5:137273/149‑1 (MQ=255)       cAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATcc                                                                                                                                        <  4:274982/149‑1 (MQ=255)       cAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATcc                                                                                                                                        <  5:196114/149‑1 (MQ=255)         gTCCTGGGCTCAGAGCCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑Gc                                                                                                                                     >  7:316747/1‑149 (MQ=255)            cTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTg                                                                                                                                  >  3:190143/1‑149 (MQ=255)            cTGGGCTCAGAACCGTCGCGTAGAACTGCGAAAGTAAGTTCTTATGCGTATGTGCCGCCGAGTAGTAACCGTCCACGAACACAGCCTGCCGCGTGCGGCCTGGGCGGAGGACCCGGAAGGAGATGACAACGCAGGCAGCAAACC‑Gccgg                                                                                                                                  >  8:57149/1‑147 (MQ=255)            cTGGGCTCAGAACCATCACGTAAAACTGCGTAAATAAAATCTTATACATATGTGCCACCGTATAATAACCGTCATCGAACACAGCCTGCCGATTGCGGCATGGGCTGAGGTCCCTGTAGTAGATGACAACCCAGGCAGCGAACC‑Gcagg                                                                                                                                  >  8:56022/1‑146 (MQ=255)             tGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGc                                                                                                                                 >  8:136958/1‑149 (MQ=255)               ggCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCgg                                                                                                                               >  6:200030/1‑149 (MQ=255)               ggCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGACTGAAGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCgg                                                                                                                               <  1:244598/149‑1 (MQ=255)                gCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCggg                                                                                                                              <  8:316747/149‑1 (MQ=255)                 cTCAGAACCGTCGCGTAGAACTGTGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGt                                                                                                                             >  2:284703/1‑149 (MQ=255)                  tCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCGGATAGCT‑GCCTGCgt                                                                                                                               >  2:303288/1‑145 (MQ=255)                  tCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCGGATAGCT‑GCCTGCgt                                                                                                                               <  1:303288/146‑2 (MQ=255)                  tCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGtt                                                                                                                            <  3:253124/149‑1 (MQ=255)                  tCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGt                                                                                                                             >  4:307328/1‑148 (MQ=255)                         cGTCGCGTAGAACTGTGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCAc                                                                                                                     <  1:284703/149‑1 (MQ=255)                         cGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCAc                                                                                                                     <  5:74432/149‑1 (MQ=255)                               gTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAgcgg                                                                                                               <  5:200030/149‑1 (MQ=255)                               gTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAgcgg                                                                                                               >  6:17854/1‑149 (MQ=255)                                    aCTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGcc                                                                                                           >  7:153568/1‑148 (MQ=255)                                      tGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTAt                                                                                                        >  8:20657/1‑149 (MQ=255)                                             gtaagtTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGt                                                                                                 >  2:92106/1‑149 (MQ=255)                                                  ttCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAggg                                                                                            <  7:180012/149‑1 (MQ=255)                                                  ttCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAggg                                                                                            >  2:35565/1‑149 (MQ=255)                                                   tCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAgggg                                                                                           >  8:28518/1‑149 (MQ=255)                                                     ttATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCt                                                                                         >  4:156183/1‑149 (MQ=255)                                                     ttATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCt                                                                                         >  2:42302/1‑149 (MQ=255)                                                     ttATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCCGGTCAGGGGCt                                                                                         >  6:224002/1‑149 (MQ=255)                                                        tGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCg                                                                                      >  6:144723/1‑149 (MQ=255)                                                             aTGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCG‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGcccc                                                                                 >  2:282080/1‑149 (MQ=255)                                                                tgCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGc                                                                              <  7:346434/149‑1 (MQ=255)                                                                   cgccgTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCctc                                                                           <  1:188608/149‑1 (MQ=255)                                                                    gccgTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCctct                                                                          <  6:103593/149‑1 (MQ=255)                                                                    gccgTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCctct                                                                          <  5:17854/149‑1 (MQ=255)                                                                    gccgTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCctct                                                                          >  5:205563/1‑149 (MQ=255)                                                                    gccgTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCctct                                                                          <  7:56022/149‑1 (MQ=255)                                                                    gccgGGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCctct                                                                          <  7:57149/149‑1 (MQ=255)                                                                           agtaACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCG‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAg                                                                   <  1:282080/148‑1 (MQ=255)                                                                             taACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGc                                                                >  5:56164/1‑149 (MQ=255)                                                                             taACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGc                                                                <  4:190143/149‑1 (MQ=255)                                                                               aCCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCa                                                              <  6:355660/149‑1 (MQ=255)                                                                                 cGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGc                                                            >  4:50740/1‑149 (MQ=255)                                                                                     cTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCggg                                                                                                  <  5:293972/108‑1 (MQ=255)                                                                                     cTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCggg                                                                                                  >  6:293972/1‑108 (MQ=255)                                                                                     cTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGt                                                         <  7:136958/148‑1 (MQ=255)                                                                                      tCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGtt                                                        >  8:295411/1‑148 (MQ=255)                                                                                       cGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCa                                                      >  5:30654/1‑149 (MQ=255)                                                                                       cGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGAAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCa                                                      >  4:237631/1‑149 (MQ=255)                                                                                          aCTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGc                                                   <  3:307328/149‑1 (MQ=255)                                                                                           cTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATg                                                    >  4:350267/1‑147 (MQ=255)                                                                                            tCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAg                                                 >  5:279390/1‑149 (MQ=255)                                                                                                   gCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTc                                                       >  3:40833/1‑136 (MQ=255)                                                                                                   gCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTc                                                       <  4:40833/136‑1 (MQ=255)                                                                                                   gCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAAc                                          <  5:144723/149‑1 (MQ=255)                                                                                                     cGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCCGGCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCaa                                           >  1:101424/1‑147 (MQ=255)                                                                                                                  ggCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAg                            <  6:30654/148‑1 (MQ=255)                                                                                                                    cTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCtt                         >  6:324838/1‑149 (MQ=255)                                                                                                                       aGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCTCGCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTc                       <  1:92106/149‑1 (MQ=255)                                                                                                                       aGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCA‑GCAAAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTcc                      >  1:113243/1‑149 (MQ=255)                                                                                                                        ggTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTccc                     <  1:42302/149‑1 (MQ=255)                                                                                                                            ccTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGc                   >  8:327174/1‑147 (MQ=255)                                                                                                                            ccTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCa                 <  3:156183/149‑1 (MQ=255)                                                                                                                                 agtagATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCaa             <  8:127741/148‑1 (MQ=255)                                                                                                                                      cataacccGCAGGCGGATAGCT‑GCCTGCGTGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAg                >  6:329045/8‑140 (MQ=255)                                                                                                                                      cataacccGCAGGCGGATAGCT‑GCCTGCGTGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAg                <  5:329045/133‑1 (MQ=255)                                                                                                                                      aTGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCAAGGcatc       <  8:153568/149‑1 (MQ=255)                                                                                                                                       tGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCAAGGcatca      >  5:376726/1‑149 (MQ=255)                                                                                                                                         aCAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCAAGGcatcat     <  1:116532/148‑1 (MQ=255)                                                                                                                                           aaCGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCAAGGCATCATCGa  <  5:224002/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              GACGAGCAGTCCTGGGCTCAGAACCGTCGCGTAGAACTGCGTAAGTAAGTTCTTATGCGTATGTGCCGCCGTGTAGTAACCGTCCTCGCACTCAGCCTGCCGCTTGCGGCCTGGGCTGAGGTCCCTGTAGTAGATGACAACGCAGGCAGCTATCC‑GCCTGCGGGTTATGGCACGAGCGGCGCCTATGCCGGGTCAGGGGCTTCGGCCCCTGCCTCT‑GCACAGGGCCAGCTGTTCATGCAGCTGCAACAGATGCAGGATCAGCTTTCCCGCCAGCAAGGCATCATCGA  >  NC_002947/1399977‑1400263

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