Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,751	G→A	L352L (CTG→CTA)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,637,751	0	G	A	95.1%	109.5 / ‑3.4	41	L352L (CTG→CTA)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (21/18);  total (22/19)
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.84e-01

GCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGC  >  NC_000913/4637512‑4637989                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  gcgATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGcc                                                                                                                                                                                                                                        <  1:462168/251‑1 (MQ=255)   gATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCCAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATggtttgcccc                                                                                                                                                                                                                                      >  1:188853/1‑242 (MQ=255)   gATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCa                                                                                                                                                                                                                                      <  1:110990/251‑1 (MQ=255)                                aaGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCg                                                                                                                                                                                                         <  1:1657/251‑1 (MQ=255)                                          accgtaaGTGTTGGCACGTTTTACCGACAAAATTCTGACGCAAAATGCGCGTATGCAGGTATTTGTCGAAACGTTACTACGCGATGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGGTACTGGAg                                                                                                                                                                                                <  2:407765/245‑1 (MQ=255)                                                     tggCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTgggg                                                                                                                                                                                    >  1:213576/1‑250 (MQ=255)                                                       gCTCGTTTTACAGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGAGCGCACCGTGCAGTTGGCCGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGATAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTgggga                                                                                                                                                                                   >  1:105340/1‑248 (MQ=255)                                                        cTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTCGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGTAGCAGGCGCTGGGGAAt                                                                                                                                                                                 >  2:187993/1‑251 (MQ=255)                                                               taccgACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGGAAAAAACATCACTCTACATGTTATGCCCACTGGGGTTAATGTTGCTGCTACACCGGCGTTACTGGGGCAGgc                                                                                                                                                                                           >  1:250135/5‑234 (MQ=255)                                                                accgACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTTGAGCAGGCGCTGGGGGATTT‑ACTggg                                                                                                                                                                        >  2:326241/4‑250 (MQ=255)                                                                    acaacaTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATGTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGtt                                                                                                                                                                                                       <  1:436354/217‑1 (MQ=255)                                                                               aCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGAttt                                                                                                                                                         <  1:326241/251‑1 (MQ=255)                                                                                 gCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCCGCGAAGCGCGCACCGGGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGGATTT‑ACTGGATAACACCATCGAttttt                                                                                                                                                       >  2:162246/1‑250 (MQ=255)                                                                                        gcgcGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCg                                                                                                                                                 <  1:183756/250‑1 (MQ=255)                                                                                             tATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACAGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGTATAACGCCATAGATTTTACCCCCGagag                                                                                                                                             >  1:77327/1‑249 (MQ=255)                                                                                                                           ccAGGCAAGACTGGAGAATTGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGGGCAGGCGCTGGGAAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCAC‑GCTAAGCGCCGAAGTGGATc                                                                                                            >  1:80466/1‑251 (MQ=255)                                                                                                                                      tGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTc                                                                                                >  1:20378/1‑251 (MQ=255)                                                                                                                                       ggAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑Gt                                                                                                 <  1:346808/249‑1 (MQ=255)                                                                                                                                       ggAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑Gt                                                                                                 <  1:10290/249‑1 (MQ=255)                                                                                                                                       ggAGAATCGTCAGGAAGTCGTTCTGACTGTTATTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTcg                                                                                               <  1:61592/251‑2 (MQ=255)                                                                                                                                              cGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTcgccctta                                                                                         <  1:162246/250‑8 (MQ=255)                                                                                                                                              cGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTaa                                                                                        <  1:227311/251‑1 (MQ=255)                                                                                                                                              cGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGGCAGCGAAGCGCGCACCGTGCCGTTGGCAGAAAAAAACATCACTCTACCTGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCCGGCGCTGGGGAGTTTCAC‑GGATAACGCCATCGATTTTACCCCCGAGAGCGGGAGCAT‑ACAGCTAAGCGCCGAAggggcaaaaggacacgtcgccctta                                                                                          >  2:399190/1‑226 (MQ=255)                                                                                                                                                       ggcGTTCTGACTGTGGTGGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGa                                                                               <  2:304305/249‑1 (MQ=255)                                                                                                                                                       gTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTggg                                                                               >  1:142718/1‑250 (MQ=255)                                                                                                                                                          gTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGGTAACGCCATCGATTTTACCCCCGAGAGCGGGCGCATATC‑GCTAAGCGCCGAAGGGGATCAGGAACAC‑GTCGCCCTTAAGGTGcgggttac                                                                            >  2:261996/1‑244 (MQ=255)                                                                                                                                                            tCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCCCTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGGGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGGGAGCGGTCGCCTAAC‑GCTAAGCGCCCAAGTGGATCAGGAA‑ACAGTCCCCCTTAAGGTGCTGGATACCg                                                                          >  2:297947/1‑251 (MQ=255)                                                                                                                                                             cTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCgaga                                                                                                                                              >  1:376445/1‑184 (MQ=255)                                                                                                                                                             cTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCgaga                                                                                                                                              <  2:376445/184‑1 (MQ=255)                                                                                                                                                               gacgtttggtGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGATACCGGTa                                                                       <  1:292956/242‑1 (MQ=255)                                                                                                                                                                     tgTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGGGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAAAGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑CATAAGCGCCGACGTGGATCAGGAACAC‑GTCTCCCTTAAGGTGCTGGGTACCggtggtggt                                                                  >  2:472215/1‑250 (MQ=255)                                                                                                                                                                           tgtgGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGa                                                                                                                                                            <  1:49132/156‑1 (MQ=255)                                                                                                                                                                            gtgGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTg                                                           >  1:82987/1‑250 (MQ=255)                                                                                                                                                                               gcggcATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGGCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACCC‑GTCGCCCTTAAGGGGCGGGATACCGGTCGTGGTATTCCTgcct                                                        >  2:365739/1‑247 (MQ=255)                                                                                                                                                                                       atttcgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCtgt                                               <  2:369328/247‑3 (MQ=255)                                                                                                                                                                                                       gAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTgca                                 >  1:412619/1‑247 (MQ=255)                                                                                                                                                                                                        aaGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGc                              >  1:164434/1‑251 (MQ=255)                                                                                                                                                                                                           cgcgcACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGaa                                                                                                                   >  1:9114/1‑164 (MQ=255)                                                                                                                                                                                                           cgcgcACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGaa                                                                                                                   <  2:9114/164‑1 (MQ=255)                                                                                                                                                                                                                aCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActc                      <  2:82987/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaacaTCACTCTACGTGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAAC‑GCTAAGCGCCGAAGTGGATCCGGAACAC‑GTCGCCCTTAAGGTGCTGGCTACCGGCAGTGGTATTCCTGACTCCGCGCTGTCACGCATTTTTGACCGCTTTTACTCTTTGCCGCGTGCAACTGGGc  >  2:217494/5‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  GCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTT‑ACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCAC‑GCTAAGCGCCGAAGTGGATCAGGAACAC‑GTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGC  >  NC_000913/4637512‑4637989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑