Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,694,030	C→A	6.9%	R713L (CGA→CTA) ‡	ligD ←	DNA ligase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,694,030	0	C	A	6.9%	226.5 / 6.5	72	R713L (CGA→CTA) ‡	ligD	DNA ligase D
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base A (2/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAGGT  >  NC_002947/3693885‑3694177                                                                                                                                                   |                                                                                                                                                    gcCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTGATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGcc                                                                                                                                                  >  4:125661/1‑149 (MQ=255) gcCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGcc                                                                                                                                                  >  2:59239/1‑149 (MQ=255)   ccTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                                <  2:102205/149‑1 (MQ=255)    cTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCggg                                                                                                                                               <  3:195565/149‑1 (MQ=255)    cTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCggg                                                                                                                                               <  3:191069/149‑1 (MQ=255)       gCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCATCGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGt                                                                                                                                              >  1:227646/1‑148 (MQ=255)        ccGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTc                                                                                                                                             >  7:251566/1‑147 (MQ=255)        ccGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtgg                                                                                                                                            <  7:309827/148‑3 (MQ=255)                   gTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCa                                                                                                                                >  7:83301/1‑149 (MQ=255)                    tAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCAc                                                                                                                               <  8:83301/149‑1 (MQ=255)                      gTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTGATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  3:125661/149‑1 (MQ=255)                      gTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  8:251566/149‑1 (MQ=255)                      gTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  4:106257/149‑1 (MQ=255)                      gTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  1:59239/149‑1 (MQ=255)                      gTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGACACAACCACGt                                                                                                                             <  7:17017/149‑1 (MQ=255)                        cGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTg                                                                                                                            <  8:175615/148‑1 (MQ=255)                         gATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             >  3:315354/1‑146 (MQ=255)                         gATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  4:315354/146‑1 (MQ=255)                           tGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtgg                                                                                                                                            <  3:285740/129‑3 (MQ=255)                           tGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtgg                                                                                                                                            >  4:285740/1‑127 (MQ=255)                               aaTGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCt                                                                                                                    <  8:192686/149‑1 (MQ=255)                                aTGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                                   <  1:328699/149‑1 (MQ=255)                                    gTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             >  2:12176/1‑135 (MQ=255)                                    gTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGcca                                                                                                               >  6:193202/1‑149 (MQ=255)                                     tCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccac                                                                                                              >  7:273188/1‑149 (MQ=255)                                             cgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  8:39132/127‑1 (MQ=255)                                             cgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             >  7:39132/1‑127 (MQ=255)                                               aTTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGttt                                                                                                     >  4:150495/1‑149 (MQ=255)                                                tttttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtttt                                                                                                    <  8:273188/149‑1 (MQ=255)                                                 ttttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtttt                                                                                                    >  7:93442/1‑148 (MQ=255)                                                   ttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCa                                                                                                  <  5:341159/148‑1 (MQ=255)                                                   ttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAg                                                                                                 >  4:243077/1‑149 (MQ=255)                                                   ttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGAGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCa                                                                                                  <  1:347648/148‑1 (MQ=255)                                                      ggCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGaa                                                                                              <  3:150495/149‑1 (MQ=255)                                                                 cgcgCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTATTTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGc                                                                                   >  3:74508/1‑149 (MQ=255)                                                                 cgcgCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTAAGGAACACCTTCAGGc                                                                                   >  5:196211/1‑149 (MQ=255)                                                                    gcTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCa                                                                                <  8:93442/149‑1 (MQ=255)                                                                      tCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGt                                                                              >  8:113172/1‑149 (MQ=255)                                                                            gcTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTc                                                                         >  8:327105/1‑148 (MQ=255)                                                                                cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCcag                                                                      >  7:104313/1‑147 (MQ=255)                                                                                cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTGCAGGAACACCTTCAGGCCCAGTTCGGCcag                                                                      >  6:254339/1‑147 (MQ=255)                                                                                        gCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAgg                                                            <  3:243077/149‑1 (MQ=255)                                                                                               cAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGtt                                                     >  1:337924/1‑149 (MQ=255)                                                                                                     aTTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACTACCACGTGAATACCCTTGCCACCACTGGGTTTCAGGCACCCCTTCAGGCCCAGGTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCg                                               >  1:274976/1‑149 (MQ=255)                                                                                                        gACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGcc                                             <  7:327105/148‑1 (MQ=255)                                                                                                        gACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCt                                            >  2:99641/1‑149 (MQ=255)                                                                                                            aTGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTccac                                        <  2:227646/149‑1 (MQ=255)                                                                                                                 gTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATgg                                   <  4:74508/149‑1 (MQ=255)                                                                                                                 gTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATgg                                   <  7:113172/149‑1 (MQ=255)                                                                                                                  tGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGGTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGc                                  <  1:99641/149‑1 (MQ=255)                                                                                                                        aaGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTc                              >  8:11233/1‑147 (MQ=255)                                                                                                                           tCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAggg                         <  5:193202/149‑1 (MQ=255)                                                                                                                             cTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCa                       <  2:337924/149‑1 (MQ=255)                                                                                                                                gACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAgcg                    >  3:29294/1‑149 (MQ=255)                                                                                                                                   cTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTTAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcag                 <  6:196211/149‑1 (MQ=255)                                                                                                                                      gTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATCCCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGTCCATGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggt              <  5:254339/149‑1 (MQ=255)                                                                                                                                      gTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggt              >  2:17421/1‑149 (MQ=255)                                                                                                                                        ttgaCCCGGCGAGCCGGGTCAAGGGCACTACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggt              >  1:198577/4‑147 (MQ=255)                                                                                                                                         ccACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCAGTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGa                                                              <  1:246234/98‑1 (MQ=255)                                                                                                                                         ccACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGa                                                              >  2:246234/1‑98 (MQ=255)                                                                                                                                         ccACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcag           >  1:158271/1‑149 (MQ=255)                                                                                                                                          cACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGt                                                          >  6:21020/1‑101 (MQ=255)                                                                                                                                          cACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGt                                                          <  5:21020/101‑1 (MQ=255)                                                                                                                                           aCCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcaggg         >  8:196249/1‑149 (MQ=255)                                                                                                                                           aCCCGGCGAGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcagg          <  2:198577/148‑1 (MQ=255)                                                                                                                                             ccGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcag           >  4:275911/1‑145 (MQ=255)                                                                                                                                             ccGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcag           <  3:275911/145‑1 (MQ=255)                                                                                                                                               ggCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCa     <  8:104313/149‑1 (MQ=255)                                                                                                                                               ggCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCa     <  1:17421/149‑1 (MQ=255)                                                                                                                                                gCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcaggg         >  2:13840/1‑144 (MQ=255)                                                                                                                                                gCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcaggg         <  1:13840/144‑1 (MQ=255)                                                                                                                                                  tCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAGGt  >  6:27020/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAGGT  >  NC_002947/3693885‑3694177

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