Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,694,029	T→G	7.0%	R713R (CGA→CGC) ‡	ligD ←	DNA ligase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,694,029	0	T	G	7.0%	217.2 / 5.6	71	R713R (CGA→CGC) ‡	ligD	DNA ligase D
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base G (2/3);  total (33/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 = 9.66e-01

TAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAGGT  >  NC_002947/3693881‑3694177                                                                                                                                                      |                                                                                                                                                     tAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCt                                                                                                                                                      >  3:106257/1‑149 (MQ=255)     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                                                                                                                             <  4:106257/149‑1 (MQ=255)                          gTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                             <  8:251566/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                                                                                                                                            >  4:285740/1‑127 (MQ=255)                               tGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtgg                                                                                                                                            <  3:285740/129‑3 (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           <  3:275911/145‑1 (MQ=255)                                                                                                                                                 ccGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcag           >  4:275911/1‑145 (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)                                                                                                                                                      |                                                                                                                                                     TAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCA‑CGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAGGT  >  NC_002947/3693881‑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: ‑