Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,605,702	A→C	intergenic (‑39/‑102)	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,702	0	A	C	100.0%	130.5 / NA	42	intergenic (‑39/‑102)	fhuF/yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/27);  total (15/27)

ACTTAGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTG  >  NC_000913/4605463‑4605944                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   aCTTAGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGt                                                                                                                                                                                                                                         <  2:311595/251‑1 (MQ=255)       cacaTTCGGTGATGACCATTGCGCCAGTGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGAGCGGCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTAccc                                                                                                                                                                                                                                   <  2:15286/251‑1 (MQ=255)         caTTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGAGGGATCCTGCGGCTGGAGATGCGTTCGCCCGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGAAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCg                                                                                                                                                                                                                                  >  2:596902/1‑250 (MQ=255)                  atgaCCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAg                                                                                                                                                                                                                              >  1:376726/1‑245 (MQ=255)                  atgaCCAGTGCGCCAGCGTCATGGCGTTAAGGGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTGTTGCGAGCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCTCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAg                                                                                                                                                                                                                              <  2:376726/245‑1 (MQ=255)                                      aTGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATg                                                                                                                                                                                                    <  1:157380/251‑1 (MQ=255)                                                            cATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGAttt                                                                                                                                                                                                                    >  2:613867/1‑213 (MQ=255)                                                            cATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGAttt                                                                                                                                                                                                                    <  1:613867/213‑1 (MQ=255)                                                                     ggATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatcat                                                                                                                                                                     >  1:654589/1‑251 (MQ=255)                                                                       aTAAACTCCAGCAAATGTTCACGATGTTTTGCGAGCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatcatta                                                                                                                                                                   <  2:152039/251‑1 (MQ=255)                                                                         aaaCTCCATCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatcattatc                                                                                                                                                                 <  2:360748/251‑1 (MQ=255)                                                                                       gTTCACGATGTATTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAAAGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGACTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCCATGAGAGAGATg                                                                                                                                                    >  2:31899/1‑250 (MQ=255)                                                                                               tGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGata                                                                                                                                                                                <  1:592326/214‑1 (MQ=255)                                                                                               tGTTTTGCGATCGTCGACCGCCCGACATGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGAGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGCTAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCCAGGGATTTATCTGGCTTGCAAATGATAAAAATTAACCTATGata                                                                                                                                                                                >  2:592326/1‑214 (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)                                                                                                                                                                                                                           aaaTGATAATTATTGCTAATCCTAGCGATAGGTTTACCCGAGAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCATTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTg                 <  1:466813/249‑1 (MQ=255)                                                                                                                                                                                                                                 taatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCgat         >  1:586264/1‑251 (MQ=255)                                                                                                                                                                                                                                         gCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTg  <  1:502067/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   ACTTAGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTG  >  NC_000913/4605463‑4605944

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