Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,688,523:1	+G	6.5%	coding (443/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,523	1	.	G	6.5%	245.6 / 11.3	62	coding (443/2061 nt)	pilJ	twitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base . (29/29);  new base G (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688376‑5688665                                                                                                                                                     |                                                                                                                                                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                                                                                                                                       >  8:167447/1‑148 (MQ=255)           gTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAg                                                                                                                                       >  7:33478/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                                                                                                  <  7:167447/149‑1 (MQ=255)                                               ccTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTc                                                                                                  <  4:412/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                                                                                  <  4:42433/149‑1 (MQ=255)                                                              cACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑Tcc                                                                                  <  8:33478/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                                                           <  4:43568/149‑1 (MQ=255)                                                                                     aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCAt                                                           <  7:291623/149‑1 (MQ=255)                                                                                     aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCAt                                                           <  6:63947/149‑1 (MQ=255)                                                                                       cACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCACAGTCGCGGCGTACCGATTCCAtc                                                          >  8:14406/1‑148 (MQ=255)                                                                                               cGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                             >  7:323709/1‑137 (MQ=255)                                                                                               cGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                             <  8:323709/137‑1 (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                   <  4:19844/148‑1 (MQ=255)                                                                                                                              ggCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTg                   >  8:61577/1‑148 (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)                                                                                                                                              aCCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGgcg  >  4:105670/1‑149 (MQ=255)                                                                                                                                                   gT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt              >  7:232105/1‑132 (MQ=255)                                                                                                                                                   gT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt              <  8:232105/132‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688376‑5688665

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