Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,688,520	T→A	6.5%	Q149L (CAG→CTG)	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,520	0	T	A	6.5%	187.0 / 2.9	61	Q149L (CAG→CTG)	pilJ	twitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (2/2);  total (31/30)
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.21e-01

CGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688372‑5688665                                                                                                                                                     |                                                                                                                                                    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                                                                                                                                           <  3:112527/149‑1 (MQ=255)          gATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCg                                                                                                                                           <  4:189253/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                                                                                                  <  4:412/149‑1 (MQ=255)                                                   ccTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTc                                                                                                  <  7:167447/149‑1 (MQ=255)                                                    cTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAgg                                                                                                                 <  4:219348/133‑1 (MQ=255)                                                    cTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAgg                                                                                                                 >  3:219348/1‑133 (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                                                             >  3:282184/1‑148 (MQ=255)                                                                                        gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTcc                                                             <  4:282184/148‑1 (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                                                             <  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                   <  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)                                                                                                                                                     |                                                                                                                                                    CGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTT‑TCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688372‑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: ‑