Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	97.3%	114.1 / ‑3.2	37	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (21/15);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGC‑ATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCA  >  NC_000913/3485959‑3486433                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCt                                                                                                                                                                                                                                     <  1:165228/250‑1 (MQ=255)      aaTACAGTGATGTACTGCATGTATGCAAAGGAAGTCGCATTCCCGTTCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCATAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGTCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTc                                                                                                                                                                                                                                <  2:235822/250‑1 (MQ=255)        tACATTGATGTACTGCATGTATGCAAAGGCCGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCACGAGCAAGCTTATTCAcc                                                                                                                                                                                                                             >  2:252874/1‑251 (MQ=255)        tACATTGATGTACTGCATGTATGCAAAGGACGTCCCATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACACGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                             >  1:223713/1‑251 (MQ=255)             tGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAggg                                                                                                                                                                                                                         <  2:46805/250‑1 (MQ=255)                 gTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaa                                                                                                                                                                                                                     >  1:85634/1‑250 (MQ=255)                      gCATGGATGCAAAGGACGTCACATTACCGTGCAGTACAGTTTATAGCCCCTTCCCGGGTAGCGGGAAGCATATATCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                                                                                                <  1:69542/250‑1 (MQ=255)                       cATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGACTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGACAAAGCgg                                                                                                                                                                                                              >  1:239367/1‑251 (MQ=255)                       cATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                              >  1:211352/1‑251 (MQ=255)                             tGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGGATGGTTCTTGTATCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑Gc                                                                                                                                                                                                       >  1:209574/1‑251 (MQ=255)                                       gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACAt                                                                                                                                                                                             >  1:70458/1‑251 (MQ=255)                                         cacaTTACCGAGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGTAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCAGCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATc                                                                                                                                                                                            <  1:189384/250‑1 (MQ=255)                                               aCCGAGCAGTACAGTTGGTAGCACCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTGATCTGGCTCTGGAGACAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGAACAGACTCTCGAATGGTTCTTGTCGCATTGACACATTCATAAGTACCCATCCACGAGCAAGCTTATTCACCAGGGTGAACAAGCGAACAC‑ACTATACTACAtgtttaa                                                                                                                                                                                       >  2:13832/1‑243 (MQ=255)                                                    gCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGctct                                                                                                                                                                                <  2:289552/251‑1 (MQ=255)                                                     cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGTAGCATATTTCGGCAATCCAGAGAAAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCCTGGTGCTTGGCAAACAGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGACCA‑GCTGTCCTACATCGCTAACGTCTCtg                                                                                                                                                                               >  2:206912/1‑251 (MQ=255)                                                      aGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCtgt                                                                                                                                                                              >  1:192876/1‑251 (MQ=255)                                                            gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                                                                                                                                                                        <  1:252874/251‑1 (MQ=255)                                                            gTTGATAGACCCTTCCCAGGTAGCCGGCAGCATATTTCGGCAATCCAGAGACACCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGCTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCCTAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                                                                                                                                                                        >  2:51932/1‑251 (MQ=255)                                                                atgaccccGTCCCAGGTAGCGGGAAGCAGATTTCGGTAATCCAGAGAACGCGGTGTTTTCTGGCTTTGGAGAAACCTTAAAAAAAATGATCACCGCGCCTGGTTCTGGGCAACCCGCAAACAGACCCGGCTCTCGAATGGTTCTTGGCCCATTGCCCCATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGGGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGTCTCTGTGGCATTGCTg                                                                                                                                                                    <  2:70680/247‑1 (MQ=255)                                                                     cccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCa                                                                                                                                                                >  1:281356/1‑250 (MQ=255)                                                                         cccaaGGTAGCGGGAAGCATATTTCGGCAATCCAGAGAAAGCGGAGTTATCTGGCTCGGGAGAAAGCTTATAACAGAGGAGAACCGCGCATGGTGCTTGCCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATGCATAAGTACCCATCCAAGAGCAAGCGTATTCACCAGGGTGAAAAAGCGGAAAC‑GCGGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGa                                                                                                                                                            <  1:51932/246‑1 (MQ=255)                                                                                    gggAAGCATATTCCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaa                                                                                                                                                <  2:70458/251‑1 (MQ=255)                                                                                             atTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGAt                                                                                                                                       <  1:44499/251‑1 (MQ=255)                                                                                                                     cgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTAGCTGAATCAGGGTGa                                                                                                               <  1:86351/251‑1 (MQ=255)                                                                                                                       ttATCTGGCTCTGGAGAAAGCTTATAACATAGGATAACCGCGCATGGTGCTTGGCAAACCGCACACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCACAAGCAAGCTTATTCACCAGGGTGAACAAGCAGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGAAGTGCTGATCACATACGAAGATGGTACAGAACTGATCCTCTCCTATCTGAATCAGGGTGAtt                                                                                                             >  2:209250/1‑251 (MQ=255)                                                                                                                            tGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTATCCTATCTGCATCAGGGTGATTTTAtt                                                                                                        >  2:170483/1‑251 (MQ=255)                                                                                                                                     agaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTg                                                                                               >  2:113407/1‑251 (MQ=255)                                                                                                                                             tataACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGtt                                                                                       >  1:269975/1‑251 (MQ=255)                                                                                                                                                      ggATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAggg                                                                              <  1:75041/251‑1 (MQ=255)                                                                                                                                                        aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGctc                                                                                                                                                                                 <  1:63162/150‑1 (MQ=255)                                                                                                                                                        aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGctc                                                                                                                                                                                 >  2:63162/1‑150 (MQ=255)                                                                                                                                                             cgcgCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGa                                                                        >  1:66012/1‑250 (MQ=255)                                                                                                                                                                     gTGCTTGGCAAACCGCAAACAGACCAGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATCAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGA‑GAGGGCCAGGAACGTAgcgc                                                               >  2:25401/1‑250 (MQ=255)                                                                                                                                                                                                                tGTCTCATTGCCACATTCAAAAGTACCCATCCAAGAGAAAGCTTATTCACAAGGGTAAAAAAGCGGAAAC‑GCTGTACTACATCGTAAAAGGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGGTTAAAGAAATGATCCCCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGC‑ATGGGTCCGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTc                    <  1:113407/250‑1 (MQ=255)                                                                                                                                                                                                                  cctcATTGCCAAATCCATAAGTACCCCCCCAAAAGCAAGCTTAT‑CCCCCGGGGAAAAAAACGGAAACCGCTGTCCAACATCTTTAAAGCCGCGGGGGCAGGCCTAATCAAAGACGAAGAGGGAAAAGAAATGTTCCCCTCCTATATGAATCAGGGTGTTTTTTGTGGCCAACTGGGCTGGTGTGAAGAGGGCCGGGACCTTA‑CGCCAGGGCTACGGACAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTa                 <  2:209574/250‑1 (MQ=255)                                                                                                                                                                                                                   ctcATTGACACATTCATAAGTACCCATCCAAGAGCCAGCTTATTAACCAGGGTGTAAAAGCGGAAAC‑GCTGGACTACATCGTTAAAGGCTCTATGGCAGTGCTGATCAACGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGAGATTTTAATGGAGAACTGGGCCTGTTTAACGAGAGCCAGGAACATAGCGC‑ATGGGTACGTGCGAAAACCGCCTGTGAAGTGGTTGAAATTTCGTAc                >  1:196391/1‑251 (MQ=255)                                                                                                                                                                                                                                  aTAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGC‑ATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCa  >  1:184464/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGC‑ATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCA  >  NC_000913/3485959‑3486433

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