Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,605,791	A→G	intergenic (‑128/‑13)	fhuF ← / → yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,605,791	0	A	G	100.0%	156.7 / NA	49	intergenic (‑128/‑13)	fhuF/yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/31);  total (18/31)

GTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATC  >  NC_000913/4605549‑4606032                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   gTTCACGATGTATTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAAAGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGACTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCCATGAGAGAGATg                                                                                                                                                                                                                                             >  2:31899/1‑250 (MQ=255)            tttGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTg                                                                                                                                                                                                                                  <  1:411996/250‑1 (MQ=255)               gCGATCGTCGCCCGCACGGCTTGCGCAAGCTTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCaa                                                                                                                                                                                                                               <  1:459120/250‑1 (MQ=255)                    cgtcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACg                                                                                                                                                                                                                         <  2:171138/251‑1 (MQ=255)                          ccGCACGGCTTGCGCAAGCGTGGGATCCTTCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCTGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATAGTGGTTATCATTATCAATGAGAGAGATGAAAACATGTTGCAACGTACGCTGGGc                                                                                                                                                                                                                   <  2:654589/251‑1 (MQ=255)                               cGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTgg                                                                                                                                                                                                              <  1:597782/251‑1 (MQ=255)                                       gcAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAg                                                                                                                                                                                                      <  1:509338/251‑1 (MQ=255)                                        cAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAgt                                                                                                                                                                                                     <  1:562344/251‑1 (MQ=255)                                            cGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTtg                                                                                                                                                                                                 <  1:226955/251‑1 (MQ=255)                                              tGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTATCGATAGGTTTACCCGATAGCAAGGGATTTCTCTGGCTTGCAAATTATAAAAATTATCATATGATCTTGGTTATCATTATCAATGCGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGGgtgtggtt                                                                                                                                                                                               >  2:382143/1‑247 (MQ=255)                                                  aTCCTGCGGCTGGCGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGAGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATAAATGAGAGAGATAAAATCATGTTg                                                                                                                                                                                                                                  >  2:510791/1‑212 (MQ=255)                                                  aTCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTg                                                                                                                                                                                                                                  <  1:510791/212‑1 (MQ=255)                                                                 aTGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGcc                                                                                                                                                                            >  1:376683/1‑250 (MQ=255)                                                                   gCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGCACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCACGGGATTTATCTGGCTTGCAAATGATACAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCGGGGCAGGGGCTGGGGGGTGTTGCTGCCGGGATTGCTGATTGCAgg                                                                                                                                                                          >  2:123983/1‑251 (MQ=255)                                                                    cGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCa                                                                                                                                                                                                                  <  1:217860/210‑1 (MQ=255)                                                                        cGCCAGATGACAGCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATGTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGa                                                                                                                                                                     <  2:120255/251‑1 (MQ=255)                                                                                    tCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGAtt                                                                                                                                                         >  1:33940/1‑251 (MQ=255)                                                                                       tCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGAttt                                                                                                                                                        <  1:662842/249‑1 (MQ=255)                                                                                                             aGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGAGGTGTTGCTGCTGGGATTGCTGATTGCAGGGCTGGTGTATGCGGATTTATCGTCAGATCAGGTGCGGATTgg                                                                                                                                >  2:579153/1‑250 (MQ=255)                                                                                                              ggCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTc                                                                                                                               >  1:268595/1‑251 (MQ=255)                                                                                                                  aTAATCGGGATAGTAATCTAAATGATAATGATTGCTAACCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAA‑TTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCt                                                                                                                          >  1:539490/1‑251 (MQ=255)                                                                                                                    aaTCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCt                                                                                                                          <  1:596902/250‑1 (MQ=255)                                                                                                                         ggATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATggg                                                                                                                    <  1:519487/251‑1 (MQ=255)                                                                                                                           aTAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTAGCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCATCGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCt                                                                                                                                                                       <  1:60519/198‑1 (MQ=255)                                                                                                                             aGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTa                                                                                                                <  1:579153/251‑1 (MQ=255)                                                                                                                             aGTAATATAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAATGGATTTAATTGGCTTGCAAATGATAAAAATTATTATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTa                                                                                                                <  2:530104/251‑1 (MQ=255)                                                                                                                                     aaaTGATAATTATTGCTAATCCTAGCGATAGGTTTACCCGAGAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCATTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑g                                                                                                         <  1:466813/249‑1 (MQ=255)                                                                                                                                           taatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCgat                                                                                                 >  1:586264/1‑251 (MQ=255)                                                                                                                                                   gCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTg                                                                                          <  1:502067/250‑1 (MQ=255)                                                                                                                                                         cctAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAACAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGGTCAGTGGCGGATTGTCATTCTGGTGGGATTAGTATT‑GACGCCGATGATGCTGTATccc                                                                                    >  2:414233/3‑248 (MQ=255)                                                                                                                                                           tAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCt                                                                                                                          >  1:262549/1‑211 (MQ=255)                                                                                                                                                           tAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCt                                                                                                                          <  2:262549/211‑1 (MQ=255)                                                                                                                                                                    gTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCTGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGctg                                                                        <  2:33940/251‑3 (MQ=255)                                                                                                                                                                      ttACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGc                                                                       <  1:414233/250‑1 (MQ=255)                                                                                                                                                                             aTAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACAtt                                                               <  1:382143/251‑1 (MQ=255)                                                                                                                                                                               aGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACAttt                                                              <  1:589070/250‑1 (MQ=255)                                                                                                                                                                                               ggcgtGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGTAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGAGGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCTGCATTACATTTTGCTATCATCGTGCCTg                                             <  2:586264/247‑1 (MQ=255)                                                                                                                                                                                                      aaaTGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTc                                                                                                                               >  1:652543/1‑163 (MQ=255)                                                                                                                                                                                                      aaaTGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTc                                                                                                                               <  2:652543/163‑1 (MQ=255)                                                                                                                                                                                                            tAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGACGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGGTTAGTATT‑GTCGACGATGCTGCTGGATTACAAACCGTTGCGGAATTCCATTTTTCTAACATCTTGCCTGGCCCGTATTggtg                                >  2:96704/1‑248 (MQ=255)                                                                                                                                                                                                             aaaaaTTATCATATGATATTGGTTATCATTATCAAGGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTgg                               <  1:440165/251‑1 (MQ=255)                                                                                                                                                                                                                    aTCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTGGCTGCCGGGAGTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGc                                                <  1:444894/227‑1 (MQ=255)                                                                                                                                                                                                                    aTCATATGATATTGGTTATCATAATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGCGTTGATGCCGGGATTGCTGATTGCCGGGCTGATGTATGCGGATTTATCGTCCGTTCAGTGGCGGGTTGTAATTCTGATGGGGTTATTATTTGCCGC‑GCTGAAGCTGTATCCCCAACCGTTTCGGCATTACATTTTGCTACCATAGTgaaggtccgaaacactgcgtgttggg                        >  2:444894/1‑226 (MQ=255)                                                                                                                                                                                                                           gATATTGGTTCTCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGATGGGCAGTGACTGGGGAGTGTTGCTGCCGGGATTGCCGCTTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGGGGCGGATTGGCATTCTGGTTGGATTAGTATT‑GACGCCGCAGAGGCTGTATCACAAAAAGATGCGGCATTACATTTTGCTACCATCGTGCCTGACCCTTATTTCTGGcccactgctcatt                  >  2:466231/1‑238 (MQ=255)                                                                                                                                                                                                                                    tatcattatcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAAt         <  1:463051/250‑1 (MQ=255)                                                                                                                                                                                                                                        attatcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGCTTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGa     >  1:504307/1‑250 (MQ=255)                                                                                                                                                                                                                                         ttatcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCg                                                   <  1:529322/203‑1 (MQ=255)                                                                                                                                                                                                                                           atcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATc  >  1:303232/1‑250 (MQ=255)                                                                                                                                                                                                                                           atcaATGAGAGAGATGAAATCATGTTGCAACGTACGCGGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATc  <  1:31899/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   GTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT‑GACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATC  >  NC_000913/4605549‑4606032

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