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

ACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485967‑3486451                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCa                                                                                                                                                                                                                                              <  1:311798/251‑1 (MQ=255) aCATTGATGGACTGCATGTATGCAAAGGAGGTCACATTCCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACAGTCATAAGTACCCAGCCAAGAGCAAGCTTAGTCACCa                                                                                                                                                                                                                                              <  2:116538/251‑1 (MQ=255)  aaTTGATGTAATGCATGTATGCAAAGGAAGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGCAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACAGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                               >  2:407792/2‑249 (MQ=255)   aTTGATGTACTGCATGTATGAAAAGGACGTCACATTACCGTGCAGTACAGTGGACAGACCATTCCAAGGTAGAGGGACGCATATTTCGGCAATACAACGACAGCGGCGTTATCTGGATCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCATATCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCAATCAAAGAGCAAGCTTATTCAACAgg                                                                                                                                                                                                                                            >  2:572552/1‑251 (MQ=255)          tACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaa                                                                                                                                                                                                                                      >  1:45721/1‑250 (MQ=255)           acagCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCTTTATCTGGCTCTGAAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaaa                                                                                                                                                                                                                                    <  2:151675/248‑1 (MQ=255)           aCTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaaa                                                                                                                                                                                                                                    <  1:179741/251‑1 (MQ=255)              gCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAg                                                                                                                                                                                                                                   <  1:447893/249‑1 (MQ=255)               cATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCGGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCTTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                                <  1:341994/251‑1 (MQ=255)                aTGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTGTAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                                >  2:140922/1‑250 (MQ=255)                   tATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAc                                                                                                                                                                                                                            >  2:513885/1‑251 (MQ=255)                        aaaGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGt                                                                                                                                                                                                                      >  1:257285/1‑251 (MQ=255)                         aaGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGt                                                                                                                                                                                                                      >  1:578991/1‑250 (MQ=255)                              ggTCACATTCCCGGGCTGTACAGTTGATGGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACa                                                                                                                                                                                                                <  2:481684/250‑1 (MQ=255)                                       aCCGTGCAGTACAGTTGATAGCCCCTGCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTaaa                                                                                                                                                                                                        <  2:161775/250‑1 (MQ=255)                                            gCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGctct                                                                                                                                                                                                 <  2:257285/251‑1 (MQ=255)                                               gtacagctgatagactctaaacAGGAAGCGGGTAGCATATTTCGGAAATCCAGAGACAGCAGCGTGATCTGACTATGGAGAACGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGGCTCTCGAATGGTTCTTGTCTCATTGCCACATCCATAAGAACCCATCCAAGAGCAAGCTTATTCTCCAGGGAGAAAAAGCGGAAAC‑GCTGTACTACAACGATAAA‑GGCACtgtg                                                                                                                                                                                              >  2:328108/22‑251 (MQ=255)                                                        aTAGCACCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAGTCTCCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTg                                                                                                                                                                                     <  2:398285/251‑1 (MQ=255)                                                         tAGCCCCTTCCCAGGTCGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTATGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCCCTGTGGCAgggatg                                                                                                                                                                                     >  2:525920/1‑245 (MQ=255)                                                             cccTTCCCAGGTAGCGGGAAGCATATTTCGGCAACCCAGAGACAGCGGCGGTATCTGGCTCTGGAGAATGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCaa                                                                                                                                                                                <  2:477453/251‑1 (MQ=255)                                                              cctccccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTATATCAGAGGATAACCGCGCATGGTTCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTGCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCaaa                                                                                                                                                                               <  1:221087/247‑1 (MQ=255)                                                                     aGGTAGCGGGAAGCATATTTCGGCAATCCAGACACAGCGGAGTTATCTGGCTCTGGAGAAAGCTTCTAACAGAGGCTAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCAGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGa                                                                                                                                                                          >  2:31213/1‑249 (MQ=255)                                                                         aGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAgg                                                                                                                                                                     <  1:513885/250‑1 (MQ=255)                                                                              gAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAg                                                                                                                                                               <  2:259659/251‑1 (MQ=255)                                                                                  aagatTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATATCCGTGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAt                                                                                                                                                           <  2:559748/248‑1 (MQ=255)                                                                                     atgtCGGCAATACAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGAt                                                                                                                                                        <  1:173432/248‑1 (MQ=255)                                                                                     atTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGAt                                                                                                                                                        >  2:260616/1‑251 (MQ=255)                                                                                         cGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGCGGGTAAAGAAATGATCctc                                                                                                                                                    >  2:78432/1‑251 (MQ=255)                                                                                           gCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCtctc                                                                                                                                                  <  1:405749/251‑1 (MQ=255)                                                                                              aTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTa                                                                                                                                               <  1:418842/251‑1 (MQ=255)                                                                                                catgagaatgcggtgTTACCTGGCTCTGGAGAAAGCTGAGAACAGAGGATAACCGCGAATGTTGTTTGGCAAACTGCAAACAGACCCGTCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAATTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGTAAAC‑GCTGTACTAAATCGTTAAA‑GGCTCTGTTGCAGTGCTTATCAAAGCCGAAGTGGGTAAAGAAATTATCCTCTCCTAt                                                                                                                                              <  2:574552/241‑1 (MQ=255)                                                                                                    agaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGa                                                                                                                                          >  2:6837/1‑250 (MQ=255)                                                                                                     gaCAGCGTCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGAGCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAAt                                                                                                                                        <  1:227469/251‑1 (MQ=255)                                                                                                                tATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtt                                                                                                                              <  1:489933/250‑1 (MQ=255)                                                                                                                  tCTGGCTCTGGAGAAAGCTTATAACCGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGACACATTCATAAGTACACATCCAAGAGCAAGCTTATTCCCCAGGGAGTCATAGCGGAATC‑GCTGTACTACATCGTTAAA‑GCCTCTGAAGCAGTGCTGGTCATAGACGAAGAAGCAAAAGACATGATCCTCTCCTATCTGACTCACGATGAttat                                                                                                                            >  1:118165/1‑248 (MQ=255)                                                                                                                           ggAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAAc                                                                                                                  <  2:494775/251‑1 (MQ=255)                                                                                                                           ggAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCg                                                                                                                     <  1:167063/248‑1 (MQ=255)                                                                                                                                       taACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                      >  1:206398/1‑251 (MQ=255)                                                                                                                                           agagGATAACCGCGCATGGTGCTGGTAAAACCGCAACCAGAACCGACTCTCGACTGGATCTTGTCTCCTTGACACATTCATAAGTACCCATCCAAGAGCAAGTTTATTCACCAGGAGGAAAACGCGGAAAC‑GCTGTACTACATCGTTCAA‑GGCTTTGAGGCAGTGCTGACCAAAGACGAAGGGGGAAAAGCAATGATCCTCTCCTTTATGAATCAGGGTGATTTTATTGGCGAACTGGTCCTGTTTGacgg                                                                                                  >  2:257446/1‑248 (MQ=255)                                                                                                                                                aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGCGGGc                                                                                              >  2:367872/1‑250 (MQ=255)                                                                                                                                                  aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                             >  1:505701/1‑249 (MQ=255)                                                                                                                                                  aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTATTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                             >  1:343200/1‑249 (MQ=255)                                                                                                                                                          aTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGt                                                                                     >  1:490093/1‑249 (MQ=255)                                                                                                                                                           gggTGCTTGGCAGACCGCAAACAGACCCGACTCTCGAATGTTTCGTGTTTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGTTCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAGTGGCGAACTGGGCCTGTTTGAAGATGGCCAGGAGCGTAgc                                                                                  <  2:288455/250‑1 (MQ=255)                                                                                                                                                            ggTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAgcg                                                                                 >  1:491710/1‑251 (MQ=255)                                                                                                                                                               gCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCAt                                                                              >  1:59166/1‑251 (MQ=255)                                                                                                                                                                 tcggCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGACACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATgg                                                                            >  1:349351/3‑251 (MQ=255)                                                                                                                                                                   ggCAAGCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATggg                                                                           <  1:518194/250‑1 (MQ=255)                                                                                                                                                                    gCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGt                                                                          >  1:387247/1‑250 (MQ=255)                                                                                                                                                                      aaaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACg                                                                       <  1:295347/251‑1 (MQ=255)                                                                                                                                                                         ccGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCAGTGCCTCATTCATAAGTACCCAGCCAAGAGCAAGCTCATTCACCAGGGGGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAATACGAAGAGGGTAAACAAATTATCCTCTCCTATCTGAAGCAGGGGGATTTTATTGGAGATGTGGGCCTGTTTGAAGAGGTGCAGGAACTTTGCGCATGGTTACGTga                                                                    <  2:62456/251‑2 (MQ=255)                                                                                                                                                                                    cccGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCt                                                         <  1:110735/251‑1 (MQ=255)                                                                                                                                                                                             cGAATGGTTCTTGTCTCAGTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTAAAAAAGCGGAAAC‑GCGGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAGGATCCTCTCCTATCTGAACCAGGGTGATTTTTTTGGCGAACTGGGCCTGTTGGAAGAGGCCCAGGAACGTAGCGAATGGGTACGTGCGAAAACCGCCTGTGAAGt                                                  <  1:78432/249‑1 (MQ=255)                                                                                                                                                                                                aTGGTTCTTGTCTCATTGCCACATTCAAAAGTCCCCATCCAAGAGCAAGCTTATTCACCGGGGTAAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCGGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAAGCAGGGTGATTTTATGGGCAAACTGGGCCTGTTTGAAGTGGGCCAGGAACGTAGCGTAGGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTg                                             <  1:31213/251‑1 (MQ=255)                                                                                                                                                                                                                        tCATAAGTACCCATCCAAGAGCAAGCTTATTCACCCGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTgg                                                >  2:579980/1‑224 (MQ=255)                                                                                                                                                                                                                        tCATAAGTACCCATCCAAGAGCAAGCTGATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCGCCGATCTGAATCAGGGTGATTTTATGGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTgg                                                <  1:579980/224‑1 (MQ=255)                                                                                                                                                                                                                            aaGTACCCACCAAAGAGAAAGCTTATTCCCCAGGGTAAAAAAGCGAAAAC‑GCGGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGGGTAAAGAAATGACCCCCTCCTATCTGAATCAGGGTGATTTTATGGGCGAACGGGGCCTGTTTGAAGGGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAtt                 <  1:6837/251‑1 (MQ=255)                                                                                                                                                                                                                              gTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAAATGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCCGAAATTTCGTACAAAAAATTTCGCCAATTg                >  2:326875/1‑250 (MQ=255)                                                                                                                                                                                                                                aCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt             >  1:392704/1‑251 (MQ=255)                                                                                                                                                                                                                                 cccATCCAAGAGCAAGATTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTc            >  2:28360/1‑251 (MQ=255)                                                                                                                                                                                                                                  ccATCAAAGAGCAAGTTTTTTCTCCAGGG‑GAAAAAGCGGAAACCGCGGTACTAATTCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAAAGGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTTATTTTATGGGAGAACTGGGCGTGTTGGAAGAGGCCCAGGAACGGAGCGCAGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCa           <  1:525920/251‑1 (MQ=255)                                                                                                                                                                                                                                   cATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAAATGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtta            >  2:215681/1‑248 (MQ=255)                                                                                                                                                                                                                                    aTCCAACAGCAACCTTATTCACAAAGGAGAAACAGCGGAAAC‑GCTGCAATACCTCGTTAAA‑GGCTCTGCGGCAGTGCTGATCAAAGACGAAGTGGGCAAAGAAATGATCCTCTCCTATATGACTCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGCACAAGTCCAGGCACGTAGCGCCTGGGTACGTGCGAAAACCGCATGTGACGTGGCTGAAATTTCGAACAAAAACTTTCGTCAATTGATTCAgg         >  2:365150/1‑251 (MQ=255)                                                                                                                                                                                                                                      cccagagCAAGCTTATTCACCAGGGTGAAAAAGCGAAAAC‑GCTGAACTAAATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGGAACGACATGATCCTCTCCTATCTGAAACAGGGAGATTTTACTGGCGCACTGGGCCTGTTTGGAGATGGACAGGCACGTCGCGCATGGGCACGTGCGAAAACCACCTGTGAAGAGGATGAACTTTCGAAAAAAAACTTTCGCCAATTGATTctgga        >  1:347870/4‑246 (MQ=255)                                                                                                                                                                                                                                           agCAAGCTTATTCCCAAGGGGGAAAAAGCGAAAAC‑GTTGTTCTACA‑CATTAAAAGGCTCTGGGGCAGTGCTGATCAAAGACAAAGAGGGAAAAAAATTGATCCTTTCCTACCTGAATCAGGGTGATTTTATTGGCGAACGGGGCCTGTTTGAAGAGGACCAGAAACGTAGCGCATGGTAATGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAccc  <  2:37374/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         ACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485967‑3486451

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