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	96.2%	172.1 / NA	52	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (28/22);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAC‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGG  >  NC_000913/3485960‑3486444                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  cAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCtt                                                                                                                                                                                                                                              >  1:372009/1‑250 (MQ=255)              aTGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTGGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGa                                                                                                                                                                                                                                <  2:507345/251‑1 (MQ=255)                     gCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCg                                                                                                                                                                                                                         >  1:340865/1‑251 (MQ=255)                      cATGTATGCGAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCCTAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                        >  2:261953/1‑251 (MQ=255)                      cATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                        >  1:542831/1‑251 (MQ=255)                           aTGCAAAGGACGTAACATTACCGAGCAGTACAGTTGATAGACCCTTACCAGGTAGCGGGAAGCATATGTCGGCAATCCAGAGACAGCGGAGTTATCTGGCTCTGGCGAACGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCCGGGGGAACAAGCGGAAACg                                                                                                                                                                                                                   >  2:477641/1‑251 (MQ=255)                             gCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCt                                                                                                                                                                                                                 >  2:371773/1‑251 (MQ=255)                             gCAAAGGACGTCACATTACCGTGCAGTACAGTGGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCTGAAACGCt                                                                                                                                                                                                                 <  1:551326/251‑1 (MQ=255)                                          aaTTACCGTGCAGGACAGTTGATAGCCCCTTCCCAGGGAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGtt                                                                                                                                                                                                    <  2:209624/250‑1 (MQ=255)                                              aCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTCCGGCAATCCAGAGACAGCGGCGTTATCTGGCGCTGGAAAAAGCTTATAACAGAAGATAACCGCGCATGGTGCTTGTCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTaaa                                                                                                                                                                                                 <  2:372009/250‑1 (MQ=255)                                                    cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCtg                                                                                                                                                                                         <  1:684327/251‑1 (MQ=255)                                                        aCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGc                                                                                                                                                                                     >  1:657545/1‑251 (MQ=255)                                                          aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATCTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCCTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGc                                                                                                                                                                                     >  1:412968/1‑249 (MQ=255)                                                                                                  gCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAACTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGAc                                                                                                                                                                     <  1:165767/225‑1 (MQ=255)                                                                                                  gCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAACTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGAc                                                                                                                                                                     >  2:165767/1‑225 (MQ=255)                                                                                                    aaTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCt                                                                                                                                         >  2:368914/1‑251 (MQ=255)                                                                                                    aaTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCt                                                                                                                                         >  1:163271/1‑251 (MQ=255)                                                                                                    aaTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGAAAACCGCGCATGGTGCTTGGCAAACCTCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGTAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCGGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCt                                                                                                                                         <  1:656953/251‑1 (MQ=255)                                                                                                        cAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCAGCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTGCTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCt                                                                                                                                     <  1:261953/251‑1 (MQ=255)                                                                                                             aCAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                                                                                                <  2:340865/251‑1 (MQ=255)                                                                                                                gcggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAgg                                                                                                                             <  1:499611/251‑1 (MQ=255)                                                                                                                   gcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTg                                                                                                                          <  1:371773/251‑1 (MQ=255)                                                                                                                        aTCTGGCTCTGGCGAAAGCTTACACCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATACACGAGCAC‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTCAAGAAATGCTCCTCTCCTATCTGAATCAGGGTGAtttt                                                                                                                     >  2:319877/1‑251 (MQ=255)                                                                                                                           tGGCTCTGGAGAAAGCTTATAACAGTGGATAACCGCGCATGGTGCTTGGCAAAACGCAAACAGACTCGACTCTAGAATGGTTCTTGTCTCATTGACACATTCATAAGTACACATCCAAGAG‑ACAGCTTATTCAACAGGGTGAAAAAGCGGAACCGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTAtt                                                                                                                  >  1:325358/1‑251 (MQ=255)                                                                                                                            ggCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAt                                                                                                                   <  1:106778/249‑1 (MQ=255)                                                                                                                                ctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGc                                                                                                               <  1:111316/249‑1 (MQ=255)                                                                                                                                   gagaAAGCTTATAACAGCGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACt                                                                                                          >  2:50994/1‑251 (MQ=255)                                                                                                                                     gaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAgg                                                                                                                                                                                                                                    >  1:439582/1‑128 (MQ=255)                                                                                                                                     gaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCGTATTCACCAgg                                                                                                                                                                                                                                    <  2:439582/128‑1 (MQ=255)                                                                                                                                      aaaGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGAGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTggg                                                                                                       >  2:496059/1‑251 (MQ=255)                                                                                                                                          cTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTg                                                                                                   >  1:295893/1‑251 (MQ=255)                                                                                                                                           ttATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGt                                                                                                  >  1:366040/1‑251 (MQ=255)                                                                                                                                            tataACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGt                                                                                                  <  1:318783/250‑1 (MQ=255)                                                                                                                                            taaaccAGAGGATAACCGCCCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTGCGTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCGTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CGCTGTGGCAGTGCTGATCAAAGGCGAAGAGGGTAAAGAAATGATCCACTCCTATCTGAAGCAGGCTGATTTGATTTGCGAACTGGGCCtt                                                                                                   <  2:476787/244‑2 (MQ=255)                                                                                                                                             ataACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGtt                                                                                                 <  2:568576/250‑1 (MQ=255)                                                                                                                                                    agGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAgg                                                                                         >  2:8043/1‑251 (MQ=255)                                                                                                                                                       aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAGGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                      >  1:606438/1‑251 (MQ=255)                                                                                                                                                         aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                      >  1:528696/1‑249 (MQ=255)                                                                                                                                                            cgcgCATGGTGCTTGGCAAAACGCAAACCGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGCAGAGGGCCAGGa                                                                                  >  2:63216/1‑250 (MQ=255)                                                                                                                                                                          ggCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTTAAGAGGGCCAGTAACGTAGCGCATGGGt                                                                   <  2:559701/251‑1 (MQ=255)                                                                                                                                                                            cAAACCGAAAAAAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCACCATTCATAAGTACCAATCCAAGAGCAA‑GCTTATTCCCCAGAGTGAAAAAGCGGAAACGCTGTACTACATCGTTACAGG‑CTCTGTGGCAGTGAAGATCAAAGAAGAAGCGGGTCAAGAAATGATCCTATCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCCGTTTGAAGAGGGCCAGGAACGTCGCGCCTGGGTAc                                                                 >  2:251810/1‑251 (MQ=255)                                                                                                                                                                             aaaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACg                                                                <  2:542831/251‑1 (MQ=255)                                                                                                                                                                                  gCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCAACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCg                                                            >  2:492786/1‑250 (MQ=255)                                                                                                                                                                                                   tcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGAAA‑GATTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTg                                          >  1:624257/1‑251 (MQ=255)                                                                                                                                                                                                       aTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTg                                      >  2:571768/1‑251 (MQ=255)                                                                                                                                                                                                        tGGTTCTTGCCTCTTTGCCACATTCATAAGTACCCACCCAAGAGCAA‑GCTTATTCACAAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGTTTTTATTGGCGAACTGGCCCTGTTTGAAGAGGGCCAGGAACGTAGCGCAGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTg                                      <  1:571768/250‑1 (MQ=255)                                                                                                                                                                                                            tCTTGTCTCATTGCCACATTCATACCAACCCATCCAAGAGCAA‑GCTTATTCACAAGGGTGAAAACGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGACGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAAtt                                 >  1:559929/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaGTCCCCACCCAAGACCAA‑GCTTATTCCCCAGGGTAAAAAAGGGGAAACGCTTTATTACATCGTTAAAGG‑TTCGGGGGCAGTGCTAAACAAAGAGGAAGAGGGTAAAGAAAGGACCCCCTCTTATCTAAATAAGGGTTTTTTAATTGGCGAACTGGGCTTGTTTAAAGAGGGCAAGAAACGTAGCGAAGGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAtt          <  2:528696/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCCCATCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGAGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt      >  1:480604/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCCAAGAGCAA‑GCTTATTCACCAGGGTGAAAAAGCGGCAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAgg  >  1:462277/1‑251 (MQ=255)                                                                                                                                                                                                                                            tCCAAGAGCAA‑GTTTATTCACGAGGGAAAAAAAGGGGAAACTCTGTACTACAT‑GTTTAAGGACTTTGTGGCAGTGCTGATAAAAGACGAAGAGGGAAAAGAAAAGTTCCTCCCCTATCGGAATCAGGTTGATTTTATTGGGAAACGGGGCTTGTTGTAGGAGGGCAAGAAACGTAGCGCAGGGTTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAg   <  2:485754/249‑1 (MQ=255)                                                                                                                                                                                                                                                agagCAA‑GCTTATTAACCAGGGCGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaa                         >  2:84437/1‑223 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  CAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAC‑GCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGG‑CTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGG  >  NC_000913/3485960‑3486444

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