Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,452,122	Δ1 bp	28.8%	coding (284/1542 nt)	PP_1271 →	putative Multidrug efflux MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,452,122	0	T	.	28.8%	134.1 / 58.7	59	coding (284/1542 nt)	PP_1271	putative Multidrug efflux MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base . (9/8);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

AACGAGTCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATGCTGATGAGCGTG  >  NC_002947/1451974‑1452267                                                                                                                                                     |                                                                                                                                                     aaCGAGTCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCt                                                                                                                                                      >  5:1639/1‑149 (MQ=255) aaCGAGTCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTAGCt                                                                                                                                                      >  1:243741/1‑149 (MQ=255)   cGAGTCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTgc                                                                                                                                                    >  6:300412/1‑149 (MQ=255)     aGTCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTgc‑gg                                                                                                                                                 <  8:15969/149‑1 (MQ=255)       tCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTgc‑ggcg                                                                                                                                               <  5:300412/149‑1 (MQ=255)          gTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCtt                                                                                                                                            >  3:112821/1‑149 (MQ=255)          gTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCt                                                                                                                                             >  2:74544/1‑148 (MQ=255)          gTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCtt                                                                                                                                            >  1:118380/1‑148 (MQ=255)            aaCCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCAc                                                                                                                                         >  1:241706/1‑149 (MQ=255)              ccAAGATTTCCCTGGCCGAAATCAGCGCCGCCATCCCCATCGGTGCCGACGAAGGAGCCTGGATGCTGGCGGTAAACAGCGCCGCTACGCACTCGTCCAAGGACATTGCTCCGTGGCTGGCAAAGACCTTCTCGCAGC‑GTCGCTTCAc                                                                                                                                         >  6:97050/1‑148 (MQ=255)               cAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCAt                                                                                                                                      <  8:32913/149‑1 (MQ=255)                 aGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATg                                                                                                                                     <  8:23213/149‑1 (MQ=255)                 aGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGa                                                                                                                                    <  1:22883/149‑1 (MQ=255)                    tttCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGAGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAgcg                                                                                                                                 >  1:244092/1‑149 (MQ=255)                        ccTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAg                                                                                                                              >  8:296144/1‑149 (MQ=255)                         cTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGt                                                                                                                             <  8:135716/149‑1 (MQ=255)                          tGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTgg                                                                                                                           >  5:247308/1‑149 (MQ=255)                                   tCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCg                                                                                                                  <  7:77468/149‑1 (MQ=255)                                     cgcgGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGc                                                                                                                 <  2:135904/149‑1 (MQ=255)                                     cgcgGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGcc                                                                                                                <  4:316692/149‑1 (MQ=255)                                      gcggcgCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGGAGTGGGGGTGTTCGcc                                                                                                                >  4:252119/1‑149 (MQ=255)                                            ccATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCtt                                                                                                          >  4:201317/1‑149 (MQ=255)                                             cATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTATCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTg                                                                                                         >  6:185816/1‑149 (MQ=255)                                                  gCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCctg                                                                                                    >  4:96120/1‑149 (MQ=255)                                                   cATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCTCGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCctgc                                                                                                   >  2:92975/1‑149 (MQ=255)                                                           ccGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCg                                                                                           <  2:241706/148‑1 (MQ=255)                                                             gacgaAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGt‑‑t                                                                                       <  7:296144/149‑1 (MQ=255)                                                                gaAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGAGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGc                                                                                    >  2:194762/1‑149 (MQ=255)                                                                        cTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAAcc                                                                            >  6:293659/1‑149 (MQ=255)                                                                            ttGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGca                                                                        <  6:1639/149‑1 (MQ=255)                                                                            ttGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGca                                                                        <  3:96120/149‑1 (MQ=255)                                                                              gCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGcac                                                                       <  7:120911/148‑1 (MQ=255)                                                                                  gCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGGTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGc                                                                    <  5:243258/147‑1 (MQ=255)                                                                                  gCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGcc                                                                   <  4:112821/148‑1 (MQ=255)                                                                                  gCGGTAAACAGCGCCGCTTCGGCCGCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGGGCTTGGCGTGCTGCAGCGGT‑‑TCGCCCCCAACCTGCACAGCCt                                                                  <  2:243741/149‑1 (MQ=255)                                                                                    ggTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGAt                                                               >  6:86395/1‑149 (MQ=255)                                                                                      tataCAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGt                                                             >  5:164171/1‑149 (MQ=255)                                                                                         aCAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCctgctg                                                                                                 >  4:20578/1‑112 (MQ=255)                                                                                         aCAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCctgctg                                                                                                 <  3:20578/112‑1 (MQ=255)                                                                                           aGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑ACGCCCCCAACCTGCACAGCCTGATGTtgct                                                         <  1:194762/149‑1 (MQ=255)                                                                                            gcgcCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTtgttgc                                                       <  2:118380/149‑1 (MQ=255)                                                                                             cgccgcTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGc                                                                                                       <  2:45312/103‑1 (MQ=255)                                                                                             cgccgcTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGc                                                                                                       >  1:45312/1‑103 (MQ=255)                                                                                             cgccgcTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTtgctgc                                                       <  1:74544/149‑1 (MQ=255)                                                                                                     cGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGCGGGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACt                                                >  1:249063/1‑149 (MQ=255)                                                                                                       gTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTg                                               >  6:19362/1‑147 (MQ=255)                                                                                                         ctcGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAggg                                          >  3:247444/1‑149 (MQ=255)                                                                                                          tcGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGccc                                                                                  <  5:245480/109‑1 (MQ=255)                                                                                                          tcGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGccc                                                                                  >  6:245480/1‑109 (MQ=255)                                                                                                                    ccTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCt                                <  5:19362/149‑1 (MQ=255)                                                                                                                          cgcCGTGGCTGGCCTCGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGTCCTCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGgc                             <  1:92975/148‑1 (MQ=255)                                                                                                                           gcCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCt                         <  3:252119/149‑1 (MQ=255)                                                                                                                            ccGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCtt                        <  1:73578/149‑1 (MQ=255)                                                                                                                                   tggcCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTg                               >  1:168013/1‑134 (MQ=255)                                                                                                                                   tggcCACGACCTTCTCGC‑GC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTg                               <  2:168013/134‑1 (MQ=255)                                                                                                                                      ccACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATg              <  5:293659/149‑1 (MQ=255)                                                                                                                                         cGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATGCtg           <  7:213555/149‑1 (MQ=255)                                                                                                                                             cTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATGCtgat         >  2:244753/1‑147 (MQ=255)                                                                                                                                              ttCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TAGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGACGCCGATGCTGATGAg      >  5:45243/1‑149 (MQ=255)                                                                                                                                                  cGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATGCTGATGAGCGTg  >  6:146053/1‑149 (MQ=255)                                                                                                                                                   gCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATGCTGATGAGCGt   <  2:249063/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AACGAGTCGGTAACCAAGATTTCCCTGGCCGATATCCGCGGCGCCATGGGCATCGGTGCCGACGAAGGTGCCTGGTTGCTGGCGGTATACAGCGCCGCTTCGGTCTCGGCCATGGCCTTTGCGCCGTGGCTGGCCACGACCTTCTCGCTGC‑GGCGCTTCACCATGAGCGCAGTGGGTCTGTTCGCCGTGCTTGGCCTGCTGCAGCCGT‑‑TCGCCCCCAACCTGCACAGCCTGATGTTGCTGCGCGTACTGCAGGGCTTTGCCGCTGGCGCCTTGCCGCCGATGCTGATGAGCGTG  >  NC_002947/1451974‑1452267

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