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	100.0%	192.5 / NA	58	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/31);  total (27/31)

AATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGT  >  NC_000913/3485964‑3486445                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       aaTACATTGATGTACTGCATGTATCCAAAGGACGTCACATTACCGTGCAGTACAGTTCATAGCCCCTTCCCAGGTAGCGGGCAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCCTGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTc                                                                                                                                                                                                                                                >  2:292273/1‑250 (MQ=255)      aTTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                                          >  1:653958/1‑251 (MQ=255)          aTGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGACCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCAAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACACATTCAAGAGCAAGCTTATTCACCAGGGTg                                                                                                                                                                                                                                       >  2:356591/1‑250 (MQ=255)            gTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaa                                                                                                                                                                                                                                    <  1:174994/251‑1 (MQ=255)            gTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaa                                                                                                                                                                                                                                    <  1:205427/251‑1 (MQ=255)              aCTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAAGCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAGCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaa                                                                                                                                                                                                                                   <  1:622186/250‑1 (MQ=255)                  cATGTATGCAAAGGACGGCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTCCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATACCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCgg                                                                                                                                                                                                                             <  2:446808/251‑1 (MQ=255)                       aTGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAAc                                                                                                                                                                                                                         <  2:514425/250‑1 (MQ=255)                        tGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGCCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACg                                                                                                                                                                                                                        <  1:328951/250‑1 (MQ=255)                                    cacaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCg                                                                                                                                                                                                           >  2:99204/1‑251 (MQ=255)                                    cacaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCg                                                                                                                                                                                                           >  1:355824/1‑251 (MQ=255)                                     acaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGt                                                                                                                                                                                                          <  1:99204/251‑1 (MQ=255)                                                    ccAGTTGATAGCCCCTTCCCATGTGGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCAGGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACTCTGGACTTCAACGTTAAAGGCTCTGTgg                                                                                                                                                                                            <  2:371061/249‑1 (MQ=255)                                                    aCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                           <  1:282349/251‑1 (MQ=255)                                                     cAGTTGATAGCCCTTTCCCTGGGAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGACGTAATCTGGCTCTGGAGAAAGCTTATAACAGCGGATAACCGCGCATGGTGCTTGGAAAACCGCAAACAGAACAGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGCACCCATCCAAGAGCAAGCTTATCCACCCGGGTGACAAC‑GCGGAAACGGTGTGCAACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                           >  2:2243/1‑250 (MQ=255)                                                      aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                           >  1:89342/1‑249 (MQ=255)                                                         tGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                                                                                                                                                                                      >  2:559440/1‑251 (MQ=255)                                                          gATAGCACCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCt                                                                                                                                                                                     >  2:294208/1‑251 (MQ=255)                                                                 ccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaa                                                                                                                                                                                                                                     <  2:363959/197‑1 (MQ=255)                                                                 ccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaa                                                                                                                                                                                                                                     >  1:363959/1‑197 (MQ=255)                                                                    ccccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGAc                                                                                                                                                                          <  2:8144/250‑1 (MQ=255)                                                                     cccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGTTTATAACAGAGGAAAACCGCGCATGGTGCTTGGCAAACCGCAAACAGTCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCGGTGCTGA‑TCAAAGACg                                                                                                                                                                         <  1:244174/251‑1 (MQ=255)                                                                      ccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCCGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGCGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCCTCCAAGAGCAAGCTTATTCACCAGGGTGAACAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACg                                                                                                                                                                         >  2:214489/1‑250 (MQ=255)                                                                       cAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGa                                                                                                                                                                        >  1:171661/1‑250 (MQ=255)                                                                        aGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGAATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACgca                                                                                                                                                                       >  1:115285/1‑248 (MQ=255)                                                                             gCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGt                                                                                                                                                                 <  1:138650/251‑1 (MQ=255)                                                                             gCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGt                                                                                                                                                                 <  1:129399/251‑1 (MQ=255)                                                                              ggggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTa                                                                                                                                                                <  2:554759/250‑1 (MQ=255)                                                                                ggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTaaa                                                                                                                                                              <  2:129907/251‑1 (MQ=255)                                                                                        atTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAACTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGAt                                                                                                                                                     >  1:412363/1‑251 (MQ=255)                                                                                              gCAATCCAGAGTCATCGGCGTTATCTGGCTTAGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACATACCCGACTCGCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑tctc                                                                                                                                              <  2:470294/251‑1 (MQ=255)                                                                                                           aGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAg                                                                                                                                 <  2:264108/251‑1 (MQ=255)                                                                                                                cgTTATCGTGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGAATGGTGTTTGGCAACCCGCAACCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGa                                                                                                                            <  2:223959/251‑1 (MQ=255)                                                                                                                    atcgggCTCTGGACAAAGCTTATAACAGAGGATAACCACGAATGGTTCTTGGCAAACCACAAACAGTCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATGAGGGTGAttt                                                                                                                         <  2:441769/246‑1 (MQ=255)                                                                                                                        ggCTCTGCAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTTAtttt                                                                                                                        <  1:634510/247‑1 (MQ=255)                                                                                                                             ggaagaAAGTTTATAACAGAGAATAACCCCGCATGTTGCTTGCCAAACCCCAAACGAACCCGACTCTCGAATGGTTCTTGTCTCATTGCCCCAGTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTTATTTTATTGGCGaa                                                                                                               <  2:256804/248‑1 (MQ=255)                                                                                                                                 aaaaGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                                                                                                           <  1:404182/250‑1 (MQ=255)                                                                                                                                       ttATAACAGAGGATAACCGCGCAAGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATCGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGGGAAAAA‑GCGGAAACGCTGGACTACATCGTTAAAGGCTATGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TAGTCC‑TCTCCTATCTGAATCAGGGTGACTTGATTGGCGAACGGGGCCTG‑t                                                                                                    >  2:356084/1‑251 (MQ=255)                                                                                                                                            aCAGAGGATAAACGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGa                                                                                                >  2:186820/1‑250 (MQ=255)                                                                                                                                                     aaCCGGGCATGTTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCa                                                                                       <  2:412363/250‑1 (MQ=255)                                                                                                                                                             aTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAg                                                                              <  1:436010/251‑1 (MQ=255)                                                                                                                                                                             cGCAAACAGACCCGACTCTCGAAGGGTTCTTGTCTCATTGCCACATTCATAAGTACGCATCCAAGAGCAAGCTTGTTCACCAGGGTGAGAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGAGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑GCTCCTATATGAATCAGGGTGATTTTATTGGCGAACTGGGGCGG‑TTTGAAGAGAACCAGAAACGTAGCGAATGGG‑TACGtacg                                                             >  2:280415/1‑248 (MQ=255)                                                                                                                                                                                   cAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAAcc                                                       >  2:583628/1‑251 (MQ=255)                                                                                                                                                                                     gACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGc                                                     >  2:538442/1‑251 (MQ=255)                                                                                                                                                                                                 gattGGTTCTTGTCTCATTGCCACATTCAAAAGTCCCCATCCAAGAGCAAGCTTATCCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCTTAAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGAAAAAAAA‑TGACCC‑TCTCCTATTTGAACCAGGGTGATTTTTTTGGCGACCGGGGCTTG‑TTTGAAGGGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGc                                         <  2:171661/248‑1 (MQ=255)                                                                                                                                                                                                 gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGACTGTGAAGTGGc                                         >  2:475208/1‑251 (MQ=255)                                                                                                                                                                                                        tCTTGTCTCATTGCCACCTTCATAAGTACCCATCCAAGAGCAAGCTTATTCTCCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCa                                                                                                                                  >  1:308608/1‑157 (MQ=255)                                                                                                                                                                                                        tCTTGTCTCATTGCCACATGCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTTGCATTGCTGA‑TCAAAGACGAAGAGGGTAAATAAA‑TGATCC‑TCTCCTATCTTAATCa                                                                                                                                  <  2:308608/157‑1 (MQ=255)                                                                                                                                                                                                                 aTTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaa                          >  1:124897/1‑250 (MQ=255)                                                                                                                                                                                                                   gcccacagtcaaaagcccccaccaaaGAGCAACCTTATTCACCGGGGACAAAAG‑GCGGAAACGGTGTATAACATCGTAAAAGGTTCTGGGGCAGTGCAGA‑TCAAAGA‑GGAGGGGGTAAAAAAAATGACCT‑TCTCCTATCAGAATCGGGGTGATTTTATTGGCGA‑CGGGGCCTGGTTAGAAGAGGCCCAGGAAGTTAGCCTATGGG‑TCCGTCCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                       <  1:214489/227‑1 (MQ=255)                                                                                                                                                                                                                        aagtCATAAGTACCCATCAAAGAGAAATCTTATCCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGTCTCTGGGGCAGTGCTGA‑TCAAAGACGAAGAGGGAAAAAAAA‑TGACCC‑CCCCCTATCTGAATCAGGGTGTTTTTATTGGCGAACTGGGCCGG‑TTGGAAGAGGGCCAGAACCGTAGCGCATGGG‑TCCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTc                  <  1:408517/248‑1 (MQ=255)                                                                                                                                                                                                                             aTAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGAACAAAAAATTTCGCCa              >  1:168265/1‑250 (MQ=255)                                                                                                                                                                                                                               aagaacccaagcAAGAGCCAGCTTATTCACCAGGGAGAAAAC‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGAAGTGCTAC‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGCGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCACGGACCGTAGCGCATGGG‑AACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTCCAAAAAATTTCGCCAAtt           >  2:416596/12‑251 (MQ=255)                                                                                                                                                                                                                                  tACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCACTTGAt        >  2:392861/1‑251 (MQ=255)                                                                                                                                                                                                                                   ccccATCCAAGAGAAAGCTTATTCACCAGGGTAAAAAA‑GCGGAAACGCTTTACTACATCTTAAAAGGCTCTGTGGCAGTGCTGA‑TAAAAGACGAAGAGGTTAAAAAAA‑TTATCC‑TCTCCTATCTGAACCAGGGTGATTTTATTGGGGAATTGGTCTTG‑TTTGAAGAGGGCCAGGAACGTGGCGCATGGG‑TACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt       <  1:559440/250‑1 (MQ=255)                                                                                                                                                                                                                                        ccaaaGAGCAAGCTTATTCACCAGG‑TGAAAAAAGCGGAAACGTTGTACTACATCTTAAAAGGCTCTGTGGCATTG‑TGAATAAAAGACGAAGATGTTAAAGAAA‑T‑ATCCCTCTCCTATCTGAATCAGGGTGTTTTTATTGGTGAGCTGGGCCTG‑TTTAAAGGGGGCAAGGAACGTAGCG‑ATGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGt  <  2:513781/248‑1 (MQ=255)                                                                                                                                                                                                                                         ccAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGt  >  1:262253/1‑250 (MQ=255)                                                                                                                                                                                                                                         ccAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGGGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TCCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGt  >  2:496603/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       AATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGA‑TCAAAGACGAAGAGGGTAAAGAAA‑TGATCC‑TCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTG‑TTTGAAGAGGGCCAGGAACGTAGCGCATGGG‑TACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGT  >  NC_000913/3485964‑3486445

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