Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,957,001	T→A	6.4%	H135Q (CAT→CAA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,957,001	0	T	A	6.4%	187.4 / 3.1	62	H135Q (CAT→CAA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/2);  total (32/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 = 6.87e-01

AGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCAGCGCGATGGCCAGCT  >  NC_002947/5956857‑5957147                                                                                                                                                 |                                                                                                                                                   aGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCAGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCAc                                                                                                                                                >  8:2748/1‑149 (MQ=255)      cccAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCgg                                                                                                                                            <  8:16689/148‑1 (MQ=255)       ccAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCt                                                                                                                                          >  6:110315/1‑149 (MQ=255)       ccAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCt                                                                                                                                          >  3:70770/1‑149 (MQ=255)        cAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTATTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCtt                                                                                                                                         <  4:188273/149‑1 (MQ=255)           cGCCTGGCACGCGCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGcc                                                                                                                                      <  2:31540/149‑1 (MQ=255)                 gCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGaa                                                                                                                                  <  7:2748/147‑1 (MQ=255)                      cTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAgctgc                                                                                                                           <  4:70770/149‑1 (MQ=255)                         gTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                                                                        <  1:266733/149‑1 (MQ=255)                         gTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                                                                        >  8:283814/1‑149 (MQ=255)                         gTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                                                                        >  6:300749/1‑149 (MQ=255)                         gTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGAAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTgg                                                                                                                        >  8:34860/1‑149 (MQ=255)                                  cAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTgc                                                                                                               >  5:333682/1‑149 (MQ=255)                                  cAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTgc                                                                                                               >  1:151674/1‑149 (MQ=255)                                    aGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTgcgc                                                                                                             <  8:172692/149‑1 (MQ=255)                                     ggCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgcc                                                                                                            >  8:269016/1‑149 (MQ=255)                                     ggCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgcc                                                                                                            <  5:110315/149‑1 (MQ=255)                                     ggCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGcgcc                                                                                                            >  2:290467/1‑149 (MQ=255)                                           gCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGt                                                                                                      <  5:175882/149‑1 (MQ=255)                                                  aaCCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGaa                                                                                                                                  >  6:327029/1‑114 (MQ=255)                                                  aaCCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGaa                                                                                                                                  <  5:327029/114‑1 (MQ=255)                                                              ccGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTg                                                                                   <  4:158928/149‑1 (MQ=255)                                                                    cACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGggcgg                                                                              >  2:326406/1‑148 (MQ=255)                                                                    cACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGggcgg                                                                              <  7:34860/148‑1 (MQ=255)                                                                      cGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCaa                                                                           <  5:300749/149‑1 (MQ=255)                                                                             ttCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTggg                                                                    <  7:283814/149‑1 (MQ=255)                                                                               ccTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGgcg                                                                  >  1:130688/1‑149 (MQ=255)                                                                                cTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGgcgc                                                                 <  3:30510/149‑1 (MQ=255)                                                                                cTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCAAGCACCCGGCTTGCCCGaa                                                                                                                                  <  4:202677/84‑1 (MQ=255)                                                                                cTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCAAGCACCCGGCTTGCCCGaa                                                                                                                                  <  4:198574/84‑1 (MQ=255)                                                                                cTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCAAGCACCCGGCTTGCCCGaa                                                                                                                                  >  3:202677/1‑84 (MQ=255)                                                                                cTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCAAGCACCCGGCTTGCCCGaa                                                                                                                                  >  3:198574/1‑84 (MQ=255)                                                                                      ggCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCTGAAg                                                            <  7:6389/148‑1 (MQ=255)                                                                                             tGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGATGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACttt                                                    >  5:311260/1‑149 (MQ=255)                                                                                                gCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTc                                                 <  7:269016/149‑1 (MQ=255)                                                                                                 cTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTcc                                                <  6:333682/149‑1 (MQ=255)                                                                                                 cTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTcc                                                <  2:151674/149‑1 (MQ=255)                                                                                                 cTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGGGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTcc                                                >  1:229734/1‑149 (MQ=255)                                                                                                  tCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCt                                               >  3:258411/1‑149 (MQ=255)                                                                                                   cGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCt                                               >  1:277000/1‑148 (MQ=255)                                                                                                    gCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCt                                               <  4:258411/147‑1 (MQ=255)                                                                                                     cTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGGTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGACTGGGCGCCGAAGCACACTTTTTCCTGAt                                            <  8:200226/149‑1 (MQ=255)                                                                                                                cGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTg                                                                                   >  5:57908/1‑99 (MQ=255)                                                                                                                cGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTg                                                                                   <  6:57908/99‑1 (MQ=255)                                                                                                                cGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAgg                                 >  3:156584/1‑149 (MQ=255)                                                                                                                cGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAgg                                 >  2:57678/1‑149 (MQ=255)                                                                                                                    cacaGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCtt                             >  7:12031/1‑149 (MQ=255)                                                                                                                       aGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTACGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGc                          >  4:36858/1‑149 (MQ=255)                                                                                                                         cGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGcc                         >  4:259747/1‑148 (MQ=255)                                                                                                                            cATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAg                     <  2:277000/149‑1 (MQ=255)                                                                                                                              tGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAg                     <  8:60986/147‑1 (MQ=255)                                                                                                                              tGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAg                     >  7:60986/1‑147 (MQ=255)                                                                                                                                gACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCa                 <  1:326406/149‑1 (MQ=255)                                                                                                                                     gtgGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGTGCTTCGCCCAAGGTCAccgc             >  8:105602/1‑148 (MQ=255)                                                                                                                                       gggTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCAGCGCGAt          <  7:93310/149‑1 (MQ=255)                                                                                                                                         gtgtGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCAGCGCGATgg        >  5:123553/1‑149 (MQ=255)                                                                                                                                         gtgtGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCAGCGCGATgg        >  3:310757/1‑149 (MQ=255)                                                                                                                                           gtgCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCaa                      >  2:147883/1‑133 (MQ=255)                                                                                                                                           gtgCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCaa                      <  1:147883/133‑1 (MQ=255)                                                                                                                                             gCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCg                                                               <  7:17898/90‑1 (MQ=255)                                                                                                                                             gCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCg                                                               >  8:17898/1‑90 (MQ=255)                                                                                                                                               cATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCAGCGCGATGGCCAGCt  >  6:265241/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AGATGCCCAACGCCTGGCACGCTCGTTCACCTACAAGGCCCAGCATGGCAACCCGCCACGCCCGATCCACGGCCTGTTCCTGATGGGCAGCCTGGGCTCGCTGGCGCAGGCCGAGCACAGCGACATGGACCTGTGGGTGTGCCATGCACCCGGCTTGCCCGAACAGCTGCTGGGTGAGCTGCGCCGCAAGTGCCAGCTGCTGGAGGCCTGGGCGGCAAGCCTGGGCGCCGAAGCACACTTTTTCCTGATCGACACACAGGGCTTCGCCCAAGGTCAGCGCGATGGCCAGCT  >  NC_002947/5956857‑5957147

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