Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,987,499	C→A	intergenic (+71/+8)	ftnB → / ← yecJ	putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,499	0	C	A	100.0%	119.1 / NA	39	intergenic (+71/+8)	ftnB/yecJ	putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/22);  total (17/22)

GAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAAC  >  NC_000913/1987252‑1987713                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                       gAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGTGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGt                                                                                                                                                                                                                     <  2:47458/251‑1 (MQ=255)           atttttCAACCGTCCATGTCCTGCGCGATCTGTAAAAAGAACAGTTGCATGATTGTCTGGTGCTGCAAACCATTCTTGACGAAGTGCGCAGTGCGAAACTTTCGGGTATGTGCCCTGTGCAGACCGACTAACATGTTCTGAAGGTCGTGTCACACCAGCTGCATTGACCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTTAGAAACCGCTGCCTCAACTGTTTGAAGCAGAGGTTTTTTTAAtc                                                                                                                                                                                                           <  2:349756/247‑2 (MQ=255)                  aaCCGACAATTTCCTGCGCGATCTGGAACAAGAACAGCAGCATGATGGGCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCa                                                                                                                                                                                                    <  1:166607/251‑1 (MQ=255)                                    cgATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTg                                                                                                                                                                                  >  2:85026/1‑251 (MQ=255)                                          ggAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATGCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCa                                                                                                                                                                            <  1:85026/251‑1 (MQ=255)                                               aaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATa                                                                                                                                                                         >  1:175837/1‑249 (MQ=255)                                                                 gTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTAGTGTCACACCAGCTGCATTGATCATAATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCATCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGAt                                                                                                                                                     >  1:372920/1‑251 (MQ=255)                                                                 gTCTGTTGCTGAAAACCATTCTTGATGAAGAGTGCAGTGCGAAACTTGCGGGAATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAACCCGCTGCCTAATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGTGTTCAATAAGCGCGATGGCTTTTTTGAt                                                                                                                                                     >  1:129344/1‑251 (MQ=255)                                                                   cTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTg                                                                                                                                                   >  1:178736/1‑251 (MQ=255)                                                                        gctgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTCATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGc                                                                                                                                              >  2:359893/1‑251 (MQ=255)                                                                              aaCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg                                                                                                                                        >  1:94309/1‑251 (MQ=255)                                                                               aCCATTCTTGATGAAGTGCGCAGTGCGACACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGCAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAACCTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg                                                                                                                                        >  2:399174/1‑250 (MQ=255)                                                                                ccATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTAGTTGCCCTGTGCAGACGGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCAGCGGCGCTAATGCAGTGCGCCGATGAAGGTGTTGAGAAGCCGCTGCCTCATCTGTTTGAAGCAGAGGTTTGTTTAATGGTATTCACCCTGTGTGGTAAACTTGAGTTCAATAAGCGCGTTGGCTTTTTGGATTGCCCGCATGGTGAc                                                                                                                                      <  2:406897/251‑1 (MQ=255)                                                                                        gatgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTg                                                                                                                              >  1:83306/1‑251 (MQ=255)                                                                                         atgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgc                                                                                                                             <  1:359893/251‑1 (MQ=255)                                                                                                 cgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGcggcggcgg                                                                                                                     <  2:83306/251‑1 (MQ=255)                                                                                                   cAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGa                                                                                                                   <  1:172371/251‑1 (MQ=255)                                                                                                         gAAACTTGCGGGTCTGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCAGCATGGTGCCCGGGTCTGCGGCggggggcgggtt                                                                                                              >  2:213843/1‑244 (MQ=255)                                                                                                          aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAg                                                                                                            <  1:110827/251‑1 (MQ=255)                                                                                                           aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGt                                                                                                           <  2:94309/251‑1 (MQ=255)                                                                                                                  gggTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAATCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCg                                                                                                    <  2:295614/251‑1 (MQ=255)                                                                                                                      aTGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGGGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGAGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGatat                                                                                                >  1:399052/1‑251 (MQ=255)                                                                                                                            ccTGTTCAGACTGACCAACATGTTCTGAATTGCGTGTCACACCAGCTGTATTGATCATCAGCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCTCCCTGTGGTGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTTACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCa                                                                                          <  2:302387/251‑1 (MQ=255)                                                                                                                                           cAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGTTGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGCTGGTTATTAAAGTCGATATTCTTCAGCTGACTGGCCAttt                                                                           >  2:169068/1‑251 (MQ=255)                                                                                                                                                   tCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATAtt                                                                                               >  1:313177/1‑223 (MQ=255)                                                                                                                                                   tCTGAATGTCGTGTCACACCAGCTGCATTGATCAGCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATAtt                                                                                               <  2:313177/223‑1 (MQ=255)                                                                                                                                                    cTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGGTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGa                                                                   >  2:216879/1‑250 (MQ=255)                                                                                                                                                                  acacCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCAGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGACCTTCCACGGGcgcg                                                    >  2:270283/1‑251 (MQ=255)                                                                                                                                                                       aGCTGCATTGATCATCATCGGCGCTAATGCATTGCTCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGTGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTAc                                               <  1:233169/251‑1 (MQ=255)                                                                                                                                                                           gCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTGTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTGCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACAt                                             <  1:213843/249‑1 (MQ=255)                                                                                                                                                                                         cGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCTCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAAt                             <  1:277486/251‑1 (MQ=255)                                                                                                                                                                                                aaTGCATTGCGCCGATGAATGGTTTGAGAAACCGCTGCCCTATCTGTTTGAAGTAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGATGTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGGGACCGGGTCTGCGGCGGCGGGATGGGTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGAGTACATCGAGAACAACCAGAAtttgtttg                      <  2:343987/251‑1 (MQ=255)                                                                                                                                                                                                   gCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGAt                    >  1:88012/1‑250 (MQ=255)                                                                                                                                                                                                    cATTGCGCCGATGAAGGTTTTGAGAAACCTCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCATTAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATa                   <  2:69072/250‑1 (MQ=255)                                                                                                                                                                                                       tGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTGGATAAcc                <  1:380664/250‑1 (MQ=255)                                                                                                                                                                                                       tGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTGTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAAcc                <  1:399174/250‑1 (MQ=255)                                                                                                                                                                                                                    gTTTTGAGAAACCTCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTGTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAAc  <  1:355867/251‑1 (MQ=255)                                                                                                                                                                                                                                    gCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCa                                                             >  1:143818/1‑176 (MQ=255)                                                                                                                                                                                                                                    gCCTCATCTGTTTGAAGCAGAGGTTTTTGTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCa                                                             <  2:143818/176‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                       GAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAAC  >  NC_000913/1987252‑1987713

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