Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,688,515	Δ1 bp	6.2%	coding (451/2061 nt)	pilJ ←	twitching motility protein PilJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,688,515	0	C	.	6.2%	204.2 / 6.8	64	coding (451/2061 nt)	pilJ	twitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base . (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCA  >  NC_002947/5688367‑5688657                                                                                                                                                     |                                                                                                                                                 ttGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGc                                                                                                                                                  <  6:221594/149‑1 (MQ=255)    aCCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCa                                                                                                                                                <  6:307712/148‑1 (MQ=255)    aCCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCAc                                                                                                                                               >  4:257493/1‑149 (MQ=255)      cGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCt                                                                                                                                             >  5:63947/1‑149 (MQ=255)      cGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCt                                                                                                                                             <  7:284907/149‑1 (MQ=255)          ccAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑g                                                                                                                                        >  3:189253/1‑149 (MQ=255)              ggATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGc                                                                                                                                    >  1:6299/1‑149 (MQ=255)              ggATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGc                                                                                                                                    <  5:291308/149‑1 (MQ=255)               gATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCg                                                                                                                                   <  4:189253/149‑1 (MQ=255)               gATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCg                                                                                                                                   <  3:112527/149‑1 (MQ=255)                aTGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGa                                                                                                                                  >  3:43568/1‑149 (MQ=255)                 tGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTG‑CACCAGGTGGCAGCGAc                                                                                                                                 <  3:124787/149‑1 (MQ=255)                    gTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAg                                                                                                                               >  7:33478/1‑148 (MQ=255)                    gTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAg                                                                                                                               >  8:167447/1‑148 (MQ=255)                    gTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGc                                                                                                                              >  3:42433/1‑149 (MQ=255)                     ttCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCa                                                                                                                             <  1:172271/149‑1 (MQ=255)                         gCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                         >  3:160971/1‑149 (MQ=255)                             gCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGc                                                                                                                     >  3:412/1‑149 (MQ=255)                              cAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCt                                                                                                                    >  2:199754/1‑149 (MQ=255)                              cAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTAGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCt                                                                                                                    >  6:1325/1‑149 (MQ=255)                                             ggCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGgatga                                                                                                     >  4:24458/1‑149 (MQ=255)                                             ggCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGgat                                                                                                       <  3:257493/147‑1 (MQ=255)                                                  ccGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCg                                                                                                <  4:160971/149‑1 (MQ=255)                                                        ccTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTc                                                                                          <  4:412/149‑1 (MQ=255)                                                        ccTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTc                                                                                          <  7:167447/149‑1 (MQ=255)                                                         cTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAgg                                                                                                         >  3:219348/1‑133 (MQ=255)                                                         cTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAgg                                                                                                         <  4:219348/133‑1 (MQ=255)                                                         cTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCt                                                                                         >  8:291623/1‑149 (MQ=255)                                                          tGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCtt                                                                                        >  5:289049/1‑149 (MQ=255)                                                            aTTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTTTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                         >  5:11146/1‑114 (MQ=255)                                                            aTTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTTTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                         <  6:11146/114‑1 (MQ=255)                                                            aTTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTgc                                                                                      >  5:7934/1‑149 (MQ=255)                                                                 ccGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                 <  2:6299/149‑1 (MQ=255)                                                                      gCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGGCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑Tc                                                                           <  5:316044/149‑1 (MQ=255)                                                                       cACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑Tcc                                                                          <  8:33478/149‑1 (MQ=255)                                                                       cACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑Tcc                                                                          <  4:42433/149‑1 (MQ=255)                                                                         ccGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCa                                                                        <  3:24458/149‑1 (MQ=255)                                                                            ctctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGCCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTc                                                                     >  5:321131/1‑149 (MQ=255)                                                                                   gcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCgcggcg                                                               >  3:273536/1‑148 (MQ=255)                                                                                    cagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTa                                                             >  3:19844/1‑149 (MQ=255)                                                                                             gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGa                                                         <  6:289049/144‑1 (MQ=255)                                                                                             gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                     <  4:282184/148‑1 (MQ=255)                                                                                             gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                     >  3:282184/1‑148 (MQ=255)                                                                                              aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCAt                                                   <  6:63947/149‑1 (MQ=255)                                                                                              aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCAt                                                   <  4:43568/149‑1 (MQ=255)                                                                                              aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCAt                                                   <  7:291623/149‑1 (MQ=255)                                                                                                cACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCACAGTCGCGGCGTACCGATTCCAtc                                                  >  8:14406/1‑148 (MQ=255)                                                                                                        cGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                     <  8:323709/137‑1 (MQ=255)                                                                                                        cGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                     >  7:323709/1‑137 (MQ=255)                                                                                                         gTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCa                                        >  4:334663/1‑149 (MQ=255)                                                                                                             tCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCg                                     >  2:24194/1‑148 (MQ=255)                                                                                                                 ccTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                >  3:308762/1‑149 (MQ=255)                                                                                                                  cttcAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCTCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTg                               <  5:1325/146‑1 (MQ=255)                                                                                                                  cTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTACCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                <  1:244936/149‑1 (MQ=255)                                                                                                                     cAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGGCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGgc                             <  6:174030/148‑1 (MQ=255)                                                                                                                         tGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAgg                        >  2:29064/1‑149 (MQ=255)                                                                                                                              cACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTg                                                                                       <  4:68484/82‑1 (MQ=255)                                                                                                                              cACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTg                                                                                       >  3:68484/1‑82 (MQ=255)                                                                                                                                       ggCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTg           >  8:61577/1‑148 (MQ=255)                                                                                                                                       ggCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTg           <  4:19844/148‑1 (MQ=255)                                                                                                                                           aGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt      <  4:308762/149‑1 (MQ=255)                                                                                                                                             tgTCGCTG‑CACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAg                         >  4:169131/2‑128 (MQ=255)                                                                                                                                             tgTCGCTG‑CACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAg                         <  3:169131/127‑1 (MQ=255)                                                                                                                                               tCGCTG‑CACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCa  >  3:545/1‑149 (MQ=255)                                                                                                                                                    gCCACCTGGG‑GCAGCGCCAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCTATTCCATCTCGTCACGCACCGCAGGTGGCGCTGGCGGCAGGCCAGTGGACtt      <  8:215985/140‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCA  >  NC_002947/5688367‑5688657

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