Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,135,883	G→C	6.1%	D229E (GAC→GAG)	PP_4521 ←	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,135,883	0	G	C	6.1%	210.1 / 3.6	66	D229E (GAC→GAG)	PP_4521	aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTG  >  NC_002947/5135738‑5136028                                                                                                                                                  |                                                                                                                                                   tctcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTc                                                                                                                                                  >  4:165117/1‑148 (MQ=255) tctcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCg                                                                                                                                                 >  3:68194/1‑149 (MQ=255) tctcGCTGCCATGCCGGTCGAGAATCTGCGCCTGCGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTATACCCCGACCCCGGCATCCACCAGGCGGCTCAGGTCAGTATGCAGACGGGCTTCTTCGCTCATGATGGCCAGGTc                                                                                                                                                  >  4:91508/1‑148 (MQ=255)    cgctgcgctgCCGGTCGAGCTTCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGc                                                                                                                                              >  6:384424/1‑149 (MQ=255)     gctgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAg                                                                                                                                               <  4:306748/147‑1 (MQ=255)          gttgcCGGTCGAGGATCTGGGCCTGGGGATCGGACAAGTCGGTTGATTGCGCCGCCCGCGACTTTACGCCGACCCCGTGATCCACCAGGCGGGTCATGGGCGTTTGCAGACGGTCTTCTTCGCTGATGTTGGCCAGCTCGACCTGGcca                                                                                                                                        <  4:68194/147‑2 (MQ=255)              ccGGTCTAGCATCTGCGCCTGGGAATCGGACAAGTCTGAAGATTGCGCCGCCCGCGACTTTACCCCTACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTTATGATGGCCAGGTCGAGCTGGCCCtt                                                                                                                                      <  8:249728/147‑1 (MQ=255)              ccGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTc                                                                                                                                     <  5:128695/148‑1 (MQ=255)               cGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAg                                                                                                                                   <  5:157917/149‑1 (MQ=255)                ggTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAgg                                                                                                                                  >  3:298799/1‑149 (MQ=255)                    gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑g                                                                                                                             <  8:59503/149‑1 (MQ=255)                     aGCATCTGCGCCTGGGCATGGGTCAACTCGCATTATTGCGCCGCCCGCTACTTTAGCCGGACCGCGGCATCCTCCAGGCGGGTCAGGGCGGTTTGCAGACGGGCTTCTTGGCTGTTGTTGGCCTGGTCGAGCTGGCCCTTCAGGCCT‑gg                                                                                                                            <  3:91508/149‑1 (MQ=255)                         tCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTg                                                                                                                        >  1:103270/1‑149 (MQ=255)                            gcgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCg                                                                                                                     >  6:148709/1‑149 (MQ=255)                             cgcCTGGGCATCGGACAGGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACTTGCTGCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGc                                                                                                                    <  5:313212/149‑1 (MQ=255)                                cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTg                                                                                                                  >  8:256264/1‑148 (MQ=255)                                   ggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTa                                                                                                                <  5:148709/147‑1 (MQ=255)                                   ggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTa                                                                                                                <  1:256221/147‑1 (MQ=255)                                           aCAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAgcgc                                                                                                       >  2:34583/1‑148 (MQ=255)                                           aCAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAgcgc                                                                                                       >  3:348085/1‑148 (MQ=255)                                             aaGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGccac                                                                                                    <  2:110293/149‑1 (MQ=255)                                               gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGccacc                                                                                                   >  2:48823/1‑148 (MQ=255)                                               gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGccacc                                                                                                   >  7:32911/1‑148 (MQ=255)                                                     aGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGAc                                                                                            <  2:103270/149‑1 (MQ=255)                                                       aTTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTg                                                                                          >  6:184743/1‑149 (MQ=255)                                                        ttGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGAGGGTCAGGGCCGTTTGCACACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGTTCAGCGCCACCAACGGACtcc                                                                                         >  2:270640/1‑147 (MQ=255)                                                          gcgcCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGt                                                                                        >  2:379162/1‑148 (MQ=255)                                                             ccgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAg                                                                                    <  1:34583/149‑1 (MQ=255)                                                              cgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAgg                                                                                   >  3:337102/1‑149 (MQ=255)                                                                     gACTTTACCCCGACCCCGGCATCCAGCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTGTTCGCTGATGATGGCCAGGTAGATCTGGCCGTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTggg                                                                            <  1:270640/149‑1 (MQ=255)                                                                     gACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTggg                                                                            <  5:384424/149‑1 (MQ=255)                                                                        tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTg                                                                              <  1:120315/144‑1 (MQ=255)                                                                        tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTg                                                                              >  2:120315/1‑144 (MQ=255)                                                                           accccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                     >  6:384273/1‑134 (MQ=255)                                                                           accccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                     <  5:384273/134‑1 (MQ=255)                                                                           accccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTc                                                                      >  6:44133/1‑149 (MQ=255)                                                                             cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAg                                                                    >  5:194238/1‑149 (MQ=255)                                                                              ccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGc                                                                   >  7:45222/1‑149 (MQ=255)                                                                                  ccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAAc                                                               >  4:243959/1‑149 (MQ=255)                                                                                  ccccgGCATCCAAGAGGCGGGTCAGGGCCGTTTGAAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GTGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCaa                                                                <  5:174451/148‑1 (MQ=255)                                                                                   cccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                               <  6:19917/116‑1 (MQ=255)                                                                                   cccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                               >  5:19917/1‑116 (MQ=255)                                                                                    ccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGt                                                             <  7:256264/149‑1 (MQ=255)                                                                                     cgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTg                                                            <  3:165117/149‑1 (MQ=255)                                                                                            caccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTg                                                                                          <  4:265984/112‑1 (MQ=255)                                                                                            caccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTg                                                                                          >  3:265984/1‑112 (MQ=255)                                                                                                         aGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTg                                        <  5:184743/149‑1 (MQ=255)                                                                                                          gggCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGc                                       <  8:32911/149‑1 (MQ=255)                                                                                                          gggCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGc                                       <  1:7060/149‑1 (MQ=255)                                                                                                                    cAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGa                             <  4:348085/149‑1 (MQ=255)                                                                                                                     aGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGAc                            >  2:220697/1‑149 (MQ=255)                                                                                                                       aCGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGcc                                      >  1:347122/1‑137 (MQ=255)                                                                                                                       aCGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGcc                                      <  2:347122/137‑1 (MQ=255)                                                                                                                       aCGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGAcc                           >  6:232499/1‑148 (MQ=255)                                                                                                                         gggCTTGTTCGCTGTTGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGGCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAg                        <  5:177604/149‑1 (MQ=255)                                                                                                                           gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAg                        >  1:84137/1‑148 (MQ=255)                                                                                                                           gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTa                      >  1:253211/1‑149 (MQ=255)                                                                                                                           gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTa                      <  4:337102/149‑1 (MQ=255)                                                                                                                                  cGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCaa                  <  1:220697/146‑1 (MQ=255)                                                                                                                                  cGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc                <  5:44133/148‑1 (MQ=255)                                                                                                                                   gCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGcc               <  3:243959/148‑1 (MQ=255)                                                                                                                                   gCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCca              <  4:298799/149‑1 (MQ=255)                                                                                                                                         gatgGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag          >  2:173403/1‑147 (MQ=255)                                                                                                                                           gggCCAGCTCGACCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCt                                <  1:336430/122‑1 (MQ=255)                                                                                                                                           gggCCAGCTCGACCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCt                                >  2:336430/2‑123 (MQ=255)                                                                                                                                               cAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTg  >  6:1216/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCT‑GGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTG  >  NC_002947/5135738‑5136028

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