Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	650,230	G→A	N87N (AAC→AAT)	citE ←	citrate lyase beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	650,230	0	G	A	100.0%	202.2 / NA	61	N87N (AAC→AAT)	citE	citrate lyase beta subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/26);  total (35/26)

GCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACA‑AACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGT  >  NC_000913/649984‑650470                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  gcgcACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTa                                                                                                                                                                                                                                                <  2:142099/250‑1 (MQ=255)          gTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCAtt                                                                                                                                                                                                                                      >  1:147788/1‑251 (MQ=255)                       aaggCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCCCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCTCGAACCACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAgcgc                                                                                                                                                                                                                         <  1:373750/249‑1 (MQ=255)                        gggCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGcgcg                                                                                                                                                                                                                        >  2:249559/1‑251 (MQ=255)                             aTACCGAGCAAACGCTCGGCACCCTGACGGATTTCAACTGCGCAGGTAATCCCCACCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACACCAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATGTCGAGAACATCCTGAGCGGTATCGGTTTTAGGCAGAAGCACAACGTCCGCACCACCGCGCACGACGGCTACCAGATCAATCACTCCCCATGCGGAATCCAgcccgaggcc                                                                                                                                                                                                                   >  2:218912/1‑243 (MQ=255)                                  gATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATGCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAGCGTCCGCACCACCGCGAACGACGGCTTTCAGGTCATTACCACCCCATTCGGAATCCAGCGCGATGACCCGc                                                                                                                                                                                                               >  2:359524/1‑250 (MQ=255)                                  gATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCATCATGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATTCGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTCTCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGca                                                                                                                                                                                                              <  2:137887/251‑1 (MQ=255)                                  gATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGGAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGAACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGcc                                                                                                                                                                                                              >  2:98125/1‑250 (MQ=255)                                     cAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACaa                                                                                                                                                                                                           <  2:234494/251‑1 (MQ=255)                                       aaCGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATg                                                                                                                                                                                                         <  2:304097/251‑1 (MQ=255)                                         cGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACACTGGt                                                                                                                                                                                                       >  1:197574/1‑251 (MQ=255)                                           cTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGttt                                                                                                                                                                                                     >  2:386541/1‑251 (MQ=255)                                                    gTGAGCCATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATcgc                                                                                                                                                                                            <  2:279286/251‑1 (MQ=255)                                                       aGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGa                                                                                                                                                                                          >  1:172370/1‑250 (MQ=255)                                                              tCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCCCCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGTAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCg                                                                                                                                                                                   <  1:388211/250‑1 (MQ=255)                                                                cACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACCCCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGAt                                                                                                                                                                                <  2:397814/251‑1 (MQ=255)                                                                  cTGCGCGGGTAATGCAAAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATAGGTTTTCGGCAGACGCACAACGACAGCACCACCGCGAACGAAGGATTCCAGGTCATTAACAACCAATTCGGAAACCAGCGCGTTGGCACGCACAATGGTCTCAATATcgagtaacagaggatgt                                                                                                                                                                              >  2:38976/1‑236 (MQ=255)                                                                      gcgGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGc                                                                                                                                                                           >  1:135609/1‑250 (MQ=255)                                                                            aaTCCCCAGCGGAGATTCAATAGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGcgcgg                                                                                                                                                                    >  2:176716/1‑250 (MQ=255)                                                                               ccccAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATGTCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACTCCCCATTCGGAATCCAGCGCGTTTACACTCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTgg                                                                                                                                                                  <  1:281044/250‑1 (MQ=255)                                                                                cccAGCGGAGATTCAAGCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTa                                                                                                                                                                <  2:169811/251‑1 (MQ=255)                                                                                  cAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTaa                                                                                                                                                               >  1:1109/1‑250 (MQ=255)                                                                                   aGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAAc                                                                                                                                                             <  1:304439/251‑1 (MQ=255)                                                                                       gagaTTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTAGCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATg                                                                                                                                                         <  2:242931/251‑1 (MQ=255)                                                                                         gaTTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATgcg                                                                                                                                                       <  1:249559/251‑1 (MQ=255)                                                                                             caagcgcctcCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATgcggcg                                                                                                                                                    <  2:358136/246‑1 (MQ=255)                                                                                                 cgccgcCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGgcgggcgg                                                                                                                                               <  2:188551/251‑1 (MQ=255)                                                                                                  gccgccAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATCCCGCGGCTGTTGCATCGCGTGGTAAACCATGCGgcgggcggg                                                                                                                                              >  2:330433/1‑250 (MQ=255)                                                                                                   ccgcccgcAGGACTGTGCTACCGTGTTCACGACCACAGGCTTTTTAGATCAGCAGGGTCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCAACACCGCGCACGACGGCTTCAAGGTAATTAACACCCAATTCGGAATCCAGAGAGTTGACACGCACAATGGCTTCAATATCGCAATACAGCGGCTGTTACAGCGATCGGTAAACCATGCGgcgggcgggg                                                                                                                                             >  2:293588/1‑249 (MQ=255)                                                                                                     gccAGCAGGCCGGTGCTGCCGGGTACACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCACTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTGGCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTc                                                                                                                                           <  1:183842/251‑1 (MQ=255)                                                                                                                   gctgcCGGGTTCACGACCACAGGCTGTTTCGATTCGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCTGTGTCTTTTTCACGCAAt                                                                                                                              <  1:364035/250‑1 (MQ=255)                                                                                                                      gcCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGGCTTTTTCCCGCAATGCt                                                                                                                           >  1:282720/1‑250 (MQ=255)                                                                                                                      gcCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTGCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGAACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATGAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGc                                                                                                                            <  2:118419/249‑1 (MQ=255)                                                                                                                        cGGGTTCACGACCACATGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTACCGGTTTTCGGCAGACGCACAACTTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACg                                                                                                                        <  2:5092/251‑1 (MQ=255)                                                                                                                            ttCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTTACGCAATGCTACGGAGt                                                                                                                    >  2:194680/1‑251 (MQ=255)                                                                                                                                      caGGCTTTTTCGATCCGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAgggc                                                                                                          >  1:118721/1‑249 (MQ=255)                                                                                                                                            ttttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACAt                                                                                                    >  1:305782/1‑251 (MQ=255)                                                                                                                                                cGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTGTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAg                                                                                                 <  1:274173/250‑1 (MQ=255)                                                                                                                                                               ctTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCg                                                                                      <  2:135609/246‑1 (MQ=255)                                                                                                                                                                    tCAATATCCAGAACATCCTGACCGGTATCGGTTTTCGGCAGACGCACAACGTCAGCACCACCGCGAACGACGGCTTCCAGGACATTAACAACCCATTCGGAATCCAGAGCGTTGACACGCACAATGGTTTCAATATCGAGCTACCGCGGATGTTGCAGCGCGTGGTAAACCATGTGGCGGCCTGCGTTTTTATCCTGTATTGCTTCGTGGTCTTCGTGGTCAAACATCAGGGCATCAACCgggacganga                                                                             >  2:411722/1‑243 (MQ=255)                                                                                                                                                                     cAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGGGTCTTTTTCACGCAATGCCACGGCGGCTTCGAGGTCAAACATCAGGGCATACGCCGGGTAGATgg                                                                             >  2:82968/1‑248 (MQ=255)                                                                                                                                                                         atCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGt                                                                       <  1:225113/251‑1 (MQ=255)                                                                                                                                                                           ccAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTg                                                                     <  2:282720/251‑1 (MQ=255)                                                                                                                                                                              gAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAatat                                                                                                                                                                                               <  1:75339/126‑1 (MQ=255)                                                                                                                                                                              gAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAatat                                                                                                                                                                                               >  2:75339/1‑126 (MQ=255)                                                                                                                                                                               aaCATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGCTACAGCGGATGTTGCAGCGCGTGGTAACCCATGCGGCGGGCGGTGTCTTTTTCACGCCATGCTACGGAGTCTTCGAAGTCCAACACCAGGGCATCCGCCGGGTAGATGAAGGAGTTGCTgc                                                                 >  1:67352/1‑250 (MQ=255)                                                                                                                                                                                 cATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGAcc                                                               >  1:347700/1‑251 (MQ=255)                                                                                                                                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGGCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCAt                                                             >  1:206916/1‑250 (MQ=255)                                                                                                                                                                                      tGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTAAATATCGCGATACAGCGGAGGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATcgc                                                          >  2:381045/1‑251 (MQ=255)                                                                                                                                                                                      tGAGCGGTATCGGTTTTCGGCAGACGCACAACGACCGCACCACCGCGAACGACGACTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATAGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTGCGGCGTATTCGAGGTCAAAAATCAGGGCAACAGCCGGGTAGAGGAAGGAGTTGCTGACCATcgg                                                          >  1:438254/1‑250 (MQ=255)                                                                                                                                                                                             tATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGCAGATGAAGGAGTTGCTGACCATCGCGGCATTg                                                   >  1:50149/1‑251 (MQ=255)                                                                                                                                                                                               tCGGGTTTCGGCAGACGCACAACGTCAGCACCACCGCGAACGACGGCTTCCAGGCAATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGCTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTCAACCATGCGGCGGGCGGTGCCTTTTTCACGCCATGATACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGACCGGTCGAGGAAGGAGGTCGTGACCATAGCGGCaang                                                   >  2:412116/1‑246 (MQ=255)                                                                                                                                                                                                  gTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCAcc                                              >  1:225566/1‑251 (MQ=255)                                                                                                                                                                                                            gACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGGCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCaca‑aa                                    >  2:213556/1‑250 (MQ=255)                                                                                                                                                                                                                cacaACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCCTCCGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTAGCACCAGGCcaa‑aacaaa                                >  2:44923/1‑249 (MQ=255)                                                                                                                                                                                                                 acaACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACCCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGACGGGTAGATGAAGGAGTTGCTGACCATCGCGGCCTTGGCACCAGGCACA‑AACAACAt                              >  1:91415/1‑251 (MQ=255)                                                                                                                                                                                                                   aaCGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCCTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACA‑AACAACAtgc                            >  1:316488/1‑251 (MQ=255)                                                                                                                                                                                                                    aCGTCCGCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACA‑AACAACAtgct                           >  1:442298/1‑251 (MQ=255)                                                                                                                                                                                                                          gCACCACCGCGAACGACGGCTTCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGTGTATTCGAGGTGAAACATCAGGGGATAAGCCAGGTGGAGGAAGGAGTTGCTGACCCACGCGGCATTGGCACCAGGCACACAACA‑CATGCTgcggcg                     >  2:201228/1‑251 (MQ=255)                                                                                                                                                                                                                                        gacggtgtCCAGGTCATTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTTTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACAGCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACA‑AACAACATGCTGCGGCGGGTGCGAGTTTTAc       <  1:372177/244‑1 (MQ=255)                                                                                                                                                                                                                                               gccAGGTCATTAACACCCCAGTCGGAATCCAGCGAGTGGGCACGCACAATGGTGTCAATGTCGCGATACAGCGGATGTTGCAGCGCGTGGTAAAGCATGCGGCGGGCGGTGTTGTTGTCACGCAATGCCACGGAGTCGTCGAGGTCAACCATCAGGGCATCAGCCTAGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGTCTCCAGGCACT‑AACAACATGCTGCTGCGGGTTCGCGTTTTTCgttgt  <  2:456132/248‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  GCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTAATCCCCAGCGGAGATTCAATCGCCGCCAGCAGGCCGGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGTCCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACA‑AACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGT  >  NC_000913/649984‑650470

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