Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	203.9 / NA	55	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/29);  total (26/29)

GTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652090‑1652573                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    gtgtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCCTTCATCTTATTAAACCAGCAAACGACGCTCTCCGGAGACAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaaagagg                                                                                                                                                                                                                                            >  1:811007/1‑250 (MQ=255)     gtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                                                                                                                                                                                         <  1:239751/250‑1 (MQ=255)     gtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGGATgg                                                                                                                                                                                                                                        >  1:570218/1‑251 (MQ=255)        tCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCCTTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGAACTTAAAGAAACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                                                                                                                                                                                                       >  1:475687/1‑249 (MQ=255)                    ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGaaaaa                                                                                                                                                                                                                          <  2:395935/249‑1 (MQ=255)                              cTGGCGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGAACCAAAAACACAGGCGGGAACCAATCAGGAAATTCAATTTATTAAACAAGCAAAAGACACCCTCCGGGGCCACATTAGAATAAAGAGTCTGCGAAAAGAGCAAATCATTTTGGTTCACCTCAGAGAGTATGGCAGAA‑AAGAAAAAAAAAAATacac                                                                                                                                                                                                               >  2:451485/1‑250 (MQ=255)                                 gTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCttt                                                                                                                                                                                                           <  1:708824/251‑1 (MQ=255)                                    aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGtt                                                                                                                                                                                                        <  1:245419/251‑1 (MQ=255)                                        gtgtgCACTTGGATGTGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTc                                                                                                                                                                                                    <  2:301466/251‑1 (MQ=255)                                             cACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACCAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATTCACCTTTGTTTTTCAAccc                                                                                                                                                                                               >  1:659078/1‑250 (MQ=255)                                              aCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAACCAATCGGGGCATTCATCTTATTAAGCCAGCACTCAACGCTCCCCGGAGTCAGAAGACATTAAAGAGACTGCGCAAAGAACAAAACATTTAAGTTCACCTCCGAGAGGATGGCAGAA‑ACGAAAAACAACAAGGCACCTTTGTTTTTCaaaatg                                                                                                                                                                                              >  1:781046/1‑247 (MQ=255)                                              aCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGGAAAAAAAAAATGCACCTTTGTTTTTCaaactt                                                                                                                                                                                              >  1:330091/1‑247 (MQ=255)                                                 tGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGaa                                                                                                                                                                                            <  1:513550/250‑1 (MQ=255)                                                  ggATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAAACGGGACATTCAACTTATTAAAGCAGCAATCGACGCTCTCCGGAGGCCGATGACATTAAAGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCCGAA‑CAGAAAAACAAAAAAGCACCTTTGTTTTTTAACCtta                                                                                                                                                                                             >  2:90357/1‑246 (MQ=255)                                                   gATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCCTCTTATTAAGCCAGCAATCGACGCTCTCCGGGGTCAGAAGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAACTCCGAGAGGATGGCAGAA‑AAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgt                                                                                                                                                                                          >  1:24274/1‑250 (MQ=255)                                                    aTGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTTATCTTATTAAGCCAACAATCGACGCCCTCCGGAGGCAGATGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑AAGAAAAACAAAAATTCACCTTTGTTTTTCAAACTGAAgta                                                                                                                                                                                         >  2:244539/1‑249 (MQ=255)                                                    aTGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCATCCTGAAgtgt                                                                                                                                                                                        <  1:552867/251‑1 (MQ=255)                                                    aTGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgt                                                                                                                                                                                          <  1:244539/249‑1 (MQ=255)                                                      gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa                                                                                                                                                                                      <  2:570218/251‑1 (MQ=255)                                                             cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCtt                                                                                                                                                                                                            <  2:354386/222‑1 (MQ=255)                                                             cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAA‑CAGAAAAACAAAAATGCACCtt                                                                                                                                                                                                            >  1:354386/1‑222 (MQ=255)                                                               tAGATCTGAAAAATGGAACAATGATGATAACGAGAAATTTAGCAAAAGCCCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGGCAGGTGACATTCAAGAGACTGAGCCAAGAACATAACCTTGATGATCACCTTAGAGAGTATGGCAGAA‑CAGAAAAAAAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaagga                                                                                                                                                                              >  1:610219/1‑246 (MQ=255)                                                                aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAAACAGCAAACGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATTATTGATGTTCAACCCAGAGAGTATTGCAGAAACAAAAAAACAAAA‑TGCACCTTTGTTTTTCAACCAGaaaggtccgcgtgagga                                                                                                                                                                             >  2:391648/1‑235 (MQ=255)                                                                      gAAAAATGGAACAATGATTAACAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGGGAGTaaaa                                                                                                                                                                          <  1:559823/247‑1 (MQ=255)                                                                             ggAACAATGATGAACAGGTGAAAGTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCTGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAAAAAAACAAAAATTCACCTTTGTGTTGCAACCTGAATTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                                                <  2:696669/250‑1 (MQ=255)                                                                             ggAAAAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTAGGGTg                                                                                                                                                                <  2:659078/250‑1 (MQ=255)                                                                                  aaagatgatCAGGAGAAATGTAGCAAAAGACCGTTTCACGGTACCATAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCTGAGTCAGATGACATTAACGAGACTGAGCAAAGATCATATCATTGATGTTCACCTCAGAGAGTATGGTAGAA‑CATAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTTACCAtt                                                                                                                                                          <  2:381628/248‑1 (MQ=255)                                                                                     gatgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtttt                                                                                                                                                        <  1:409366/250‑1 (MQ=255)                                                                                          tCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCttt                                                                                                                                                                                                           <  2:767471/194‑1 (MQ=255)                                                                                          tCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCttt                                                                                                                                                                                                           >  1:767471/1‑194 (MQ=255)                                                                                                          aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAGACAGCAATCGGCGCTCTCCGGAGTCAGATGACATTAACGAGACTCAGCAAAAAGCATATCATTTATGTTTACCTCAAAGAGTATGGCAGAA‑AAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTCAAATGTCAACTAGAGTAAAAAATTATTGTGACAATTTTACCATTgcatcataaaaga                                                                                                                                     >  2:483953/1‑236 (MQ=255)                                                                                                                  gTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGAc                                                                                                                          <  1:462426/251‑1 (MQ=255)                                                                                                                       aCGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGc                                                                                                                     <  1:192368/251‑1 (MQ=255)                                                                                                                             ccAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTTTTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGGAAAGAGCATATCATTTATGTTCACCTCAGCGAGTATGGCAGAA‑CAACAAAACAAAAATGCACCTTTTTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGCCTCAATAAAGGAGATGTGGGAAGCAGCGCTaa                                                                                                                >  2:359521/1‑250 (MQ=255)                                                                                                                                 aaaCACAGGCGGGAACCAATCGGGTCACTAATCTTAATAAAACAGCAATCGACGCTCTCCGGAGACAGATGCCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAACTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTAAGCGAGAGGAAAAAATTATGGTGACCATTTTAACGTTGACTCAATAAGGCAGATGTGGGACGCAGCGCTCAAACGt                                                                                                           >  2:191390/1‑251 (MQ=255)                                                                                                                                                 cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                           <  2:646120/251‑1 (MQ=255)                                                                                                                                                          cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCat                                                                                  >  2:841382/1‑251 (MQ=255)                                                                                                                                                                       aaGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGAc                                                                      <  1:396048/250‑1 (MQ=255)                                                                                                                                                                                                     ctGACATTAACGAGACTGAGCAAAGAGCATATCATTGCTGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACATTTGTTTTTCAACCTGAAGTGTAAGAGAGAGAAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCGGGTCGATGACAGCTgg                                       >  2:424198/2‑251 (MQ=255)                                                                                                                                                                                                       gACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                                 <  1:838804/191‑1 (MQ=255)                                                                                                                                                                                                       gACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                                 >  2:838804/1‑191 (MQ=255)                                                                                                                                                                                                       gACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGt                                      <  1:443683/250‑1 (MQ=255)                                                                                                                                                                                                           ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGc                                                                                                                     >  2:578053/1‑167 (MQ=255)                                                                                                                                                                                                           ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGc                                                                                                                     <  1:578053/167‑1 (MQ=255)                                                                                                                                                                                                            tAACGAGACTGAGCAAAGCGCATATCATTGATGTACAACTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGCCTCAATAAGGCAGAGGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCAATCGAAAATCATATCAGTAGAGACATAATTATGTCTGCTGGTCCCTGACAGCTGGCGCTAAc                                >  2:622715/1‑251 (MQ=255)                                                                                                                                                                                                            tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTaa                                 <  2:475687/250‑1 (MQ=255)                                                                                                                                                                                                             aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGCGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                                <  1:495784/250‑1 (MQ=255)                                                                                                                                                                                                                     tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAttt                       >  1:639457/1‑251 (MQ=255)                                                                                                                                                                                                                             gagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCaaa               <  1:90357/251‑1 (MQ=255)                                                                                                                                                                                                                                 atatCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc             >  1:457239/1‑249 (MQ=255)                                                                                                                                                                                                                                  tatCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgc                                                        <  1:284187/205‑1 (MQ=255)                                                                                                                                                                                                                                  tatCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgc                                                        >  2:284187/1‑205 (MQ=255)                                                                                                                                                                                                                                  tatCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGAGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCGATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGCGCTAACCCGGCATTTATAGCaaaac             >  2:325342/1‑246 (MQ=255)                                                                                                                                                                                                                                          aTGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  <  1:604754/251‑1 (MQ=255)                                                                                                                                                                                                                                          aTGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  <  1:68403/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    GTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAA‑CAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652090‑1652573

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: ‑