Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,000,592	C→A	A155S (GCG→TCT)	yigE ←	DUF2233 domain‑containing protein YigE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,000,592	0	C	A	100.0%	191.3 / NA	54	A155S (GCG→TCT)	yigE	DUF2233 domain‑containing protein YigE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/22);  total (32/22)

ATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAG  >  NC_000913/4000370‑4000831                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                atGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTTCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGa                                                                                                                                                                                                                     <  1:527589/251‑1 (MQ=255)           gtgacATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGACTACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAg                                                                                                                                                                                                           <  2:508574/247‑1 (MQ=255)           gTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATcc                                                                                                                                                                                                             >  1:337119/1‑249 (MQ=255)                aTCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTGTACTGGTTTTGAAGGCATCCAGacgaa                                                                                                                                                                                                      <  2:594657/251‑1 (MQ=255)                   aaGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCTTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGacgaacga                                                                                                                                                                                                   <  1:303071/251‑1 (MQ=255)                                         ttCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGATCCCATGTTTATTAATCCCAAAACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTAAACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGAc                                                                                                                                                                             >  1:432928/1‑251 (MQ=255)                                                  gCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAacac                                                                                                                                                                    >  1:358364/1‑251 (MQ=255)                                                    tttGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGCCTTTATCTCCCGCGACATAAAACAcgc                                                                                                                                                                  >  2:108292/1‑251 (MQ=255)                                                     ttGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACTAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGTTTTTGAAGGCATCCAGCCGAACGATGCCGCCTTTATATCCCGCGACATAAAACAcgcc                                                                                                                                                                 >  1:586699/1‑251 (MQ=255)                                                      tGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACAcgcc                                                                                                                                                                 <  1:622758/250‑1 (MQ=255)                                                          aTAACAGGCAAAATCATAAAAATGTGTTGCCTGCTGGCTCAACAAAAACACGGCAGTCCCAGGTTTATTAATCCCAACTGCGGTACGGATTTTGGTTGGGGCGACGTTGGGATGAATTCGCGGATTAATTACACCTTTTTCCATCAACATTGGCTCTGGTTGCGCATGAAACTGAATCGCTTTTCAGGTTTAGAAGGCACCGTTACAAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAg                                                                                                                                                            <  1:47405/251‑1 (MQ=255)                                                                ggCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAGCACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTGTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACCCCGCCAGGACGGa                                                                                                                                                      <  2:423148/251‑1 (MQ=255)                                                                 gCAAAATCATAAAAATTTGTTTCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGAt                                                                                                                                                     <  2:525579/251‑1 (MQ=255)                                                                    aaaTCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAAGCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaaa                                                                                                                                                  <  1:151021/251‑1 (MQ=255)                                                                            aaaaTTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTcc                                                                                                                                          >  1:451698/1‑251 (MQ=255)                                                                                 ttgttgCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                                                                                      >  1:301945/1‑250 (MQ=255)                                                                                 ttgttgCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                                                                                      >  1:533274/1‑250 (MQ=255)                                                                                  tgttgCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCAc                                                                                                                                    >  2:526270/1‑251 (MQ=255)                                                                                   gttgCCTGCTGGCTCAACAAAAACACTGCGGTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCAc                                                                                                                                    <  1:479831/250‑1 (MQ=255)                                                                                     tgCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                                                                                                 >  1:71940/1‑251 (MQ=255)                                                                                      gCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGACTTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                                                                                                 >  2:224970/1‑250 (MQ=255)                                                                                       ccTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGCTAAAGAAATTCCCTTCACCTGa                                                                                                                                >  2:192643/1‑250 (MQ=255)                                                                                          gctgGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                                                                                                            >  1:631550/1‑251 (MQ=255)                                                                                          gctgGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAAGTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCATCATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                                                                                                              <  1:353796/249‑1 (MQ=255)                                                                                                     aaaaaCACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGtttt                                                                                                                                                                                                                       >  1:201049/1‑149 (MQ=255)                                                                                                     aaaaaCACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGtttt                                                                                                                                                                                                                       <  2:201049/149‑1 (MQ=255)                                                                                                       aaaCACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                                                                                                >  1:559002/1‑250 (MQ=255)                                                                                                            cGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTc                                                                                                          >  1:52571/1‑251 (MQ=255)                                                                                                             ggCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                                                                                         <  1:128377/251‑1 (MQ=255)                                                                                                                    ccATGTTTCTTAATCCCAAAACCGTTACGAATTTTGCTTGAGGCGACGTTGAGATGAATAAGCGGATTAATTCCACCGTTTTCCATCAACATTGGCCCTGTCTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCAACCAGACGAACGATGCCGTCTTTATCTCCCGCGACATAACAAACGCCGCCAGGACGGATAAAGAAATTCCCTTAACCGGAAGCGAGATTTAACGCCAACTTCTGCGGAcc                                                                                                  >  2:406747/1‑251 (MQ=255)                                                                                                                    ccATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGCCTGAACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCACGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAACGAGATTTAACGCCACCTTCTGCTGAcc                                                                                                  >  2:406734/1‑251 (MQ=255)                                                                                                                                   ccAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAACACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACa                                                                                    >  2:295105/1‑250 (MQ=255)                                                                                                                                      acacCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc                                                                                <  1:108292/251‑1 (MQ=255)                                                                                                                                         ccGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAGTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCGGCTGACCGTTTTCGATGTACAAACCGAg                                                                             <  2:337119/251‑1 (MQ=255)                                                                                                                                          cGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc                                                                            >  2:280279/1‑251 (MQ=255)                                                                                                                                               cGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc                                                                       >  1:124968/1‑251 (MQ=255)                                                                                                                                                  aTTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt                                                                     >  2:258669/1‑250 (MQ=255)                                                                                                                                                     ttGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCt                                                                 >  1:95717/1‑251 (MQ=255)                                                                                                                                                      tGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCtt                                                                >  2:15901/1‑251 (MQ=255)                                                                                                                                                               cGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTcatcat                                                         >  1:404975/1‑249 (MQ=255)                                                                                                                                                                         aTGAATACGCGGATTAATAACACCGTTTTCCATCAACATTGGCCCGGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcgg                                                                          <  1:690279/222‑1 (MQ=255)                                                                                                                                                                           gAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACACTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATCGCTTTCATCATAGATGCCGCCGTTc                                           >  2:573979/1‑251 (MQ=255)                                                                                                                                                                           gAATACACGGATTAATTACCCCGCTTTCCATCAACATTGGCCCTGACTGCACAGACAACTGACTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGCGATTTACCGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTc                                           >  2:102668/1‑251 (MQ=255)                                                                                                                                                                            aaTACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTGTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCTATGTACAAACCGAGCGGCTCATAGCTTTCATCATAGATGCCGCCGTTCa                                          <  1:630903/251‑1 (MQ=255)                                                                                                                                                                                           taCACCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAccc                           >  2:7230/1‑250 (MQ=255)                                                                                                                                                                                            acacCGTTTTCCAGCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCTcct                           <  2:432928/250‑1 (MQ=255)                                                                                                                                                                                             cacCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGacgaac                                                                                                                                                                                                     <  2:343706/79‑1 (MQ=255)                                                                                                                                                                                             cacCGTTTTCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGacgaac                                                                                                                                                                                                     >  1:343706/1‑79 (MQ=255)                                                                                                                                                                                                   tttCCATCAACATTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctg                    >  2:558202/1‑250 (MQ=255)                                                                                                                                                                                                   tttCCATCAACAGTGGCCCTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACGTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATGAAGAAATTCCCTGCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTGTTCGATGTACAAACCGAGCGTCGCATAGCTTTCATCATAGATGCCGCCGGTCATCGCCATCTGCAccttaccctg                    <  2:303406/250‑1 (MQ=255)                                                                                                                                                                                                            ataTTGGCCCTGACTGCACAGATAACTGAATTTCTTTACTGGTTTTGATGGCATCCAGACGAACGATGCCTACTTTATCTCCCGCGAAATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAATCGAGATTTAACGCCACCTTTGGCTGACCGTTTTCGATGTCTAAACCGATCGGCGCATAGCTTTCATCATATATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAata          <  1:61423/249‑1 (MQ=255)                                                                                                                                                                                                               ttGGCCCTGACTGCACAGCACACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATcc       >  2:611257/1‑251 (MQ=255)                                                                                                                                                                                                                   cccTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCAAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTAATCGCCATCTGCACCTGACCCTGCCTATTAATATCCGcc    >  2:102945/1‑250 (MQ=255)                                                                                                                                                                                                                    ccTGACTGCACAGAAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcag  <  1:462719/251‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                ATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAG  >  NC_000913/4000370‑4000831

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