Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,827	C→T	D44D (GAC→GAT)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,827	0	C	T	100.0%	197.8 / NA	68	D44D (GAC→GAT)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/39);  total (29/39)

ATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCT  >  NC_000913/4636581‑4637069                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aTCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCGACGGAGGGGACGTTAATCGATACCg                                                                                                                                                                                                                                                >  1:283836/1‑251 (MQ=255) aTCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCg                                                                                                                                                                                                                                                <  1:478890/251‑1 (MQ=255)     cACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATGTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAAc                                                                                                                                                                                                                                            <  2:312264/251‑1 (MQ=255)      cccGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCGATTTTTGTCAAAGAAGTTAAACCGTGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACg                                                                                                                                                                                                                                           <  2:290044/250‑1 (MQ=255)         gttgATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGCTATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCTGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTg                                                                                                                                                                                                                                        >  1:279495/4‑251 (MQ=255)            gATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCt                                                                                                                                                                                                                                      <  1:350513/250‑1 (MQ=255)             aTACCCACCTTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTgg                                                                                                                                                                                                                                    >  1:518781/1‑251 (MQ=255)             aTACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCTGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTgg                                                                                                                                                                                                                                    <  2:455781/251‑1 (MQ=255)              tACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGc                                                                                                                                                                                                                                   >  1:158250/1‑251 (MQ=255)                 ccACACTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGGAGAGCAACGGGGGGGACGTTAATCGCTACCGCAACGTTGATGGCggg                                                                                                                                                                                                                                >  2:417533/1‑250 (MQ=255)                   acaTTAAAACGTTGCGTGCCAGGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCAGGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTAGTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACTTTGCTGGCGGAGc                                                                                                                                                                                                                              <  2:223290/251‑1 (MQ=255)                    caTTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCt                                                                                                                                                                                                                             <  1:166902/251‑1 (MQ=255)                        aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                                                                                                                                           <  1:282133/249‑1 (MQ=255)                        aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGgc                                                                                                                                                                                                                          <  1:131973/250‑1 (MQ=255)                        aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAAACGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATTGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                                                                                                                                           >  2:282133/1‑249 (MQ=255)                           acgttGCGTGCCAAGCTGCGCGCCATCAACCCCGCTCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGGGCAACGGGGGGGACGTTAATCGATCCCGCAACGTTGCTGGCGGAGCTGGcgcgg                                                                                                                                                                                                                      >  2:573140/5‑250 (MQ=255)                            cgttGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTc                                                                                                                                                                                                                     <  2:158250/248‑1 (MQ=255)                                   ggCCAAGCTGCGCGCCATCAACCCCGAGCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCGGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGAttt                                                                                                                                                                                                              <  2:328963/250‑1 (MQ=255)                                    gCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATAGATAACGCAACGTTGCTGGCGGAGCTGGCGCGTCCCCATTTg                                                                                                                                                                                                             >  1:445098/1‑251 (MQ=255)                                                    tCAACCCCGCTCTTACACCGATTAATACTCCTCCCGGCGTGGGATATAGCCTGAGGGGCCTGTAATGCGAATCGGCCTGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGATTGTACTGGCTATTTTTGTCAAAGCAGTGCAAACGGGCGAGCGAAGAGAAACGGATGGGAAGTTACTCGATACCGCAACGTTTCTGGGGGAGCTGGCGCGGCCCGATTTGttctcgaggggcaca                                                                                                                                                                                              >  2:430539/1‑240 (MQ=255)                                                       aCCCCGCTCTTTCACCGATTAATACTCATCGCGGCCTGCGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCAATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGACGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGGGCTGGCGCGTCCCGATTTGCTCTCTGGGGGCCCAACGc                                                                                                                                                                                          >  2:270469/1‑251 (MQ=255)                                                         cccGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCAt                                                                                                                                                                                        <  1:149188/251‑1 (MQ=255)                                                              tCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGc                                                                                                                                                                                    >  1:12961/1‑250 (MQ=255)                                                                   cACCGATTAATCCTCATCGCTGCATGGTATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCATCTg                                                                                                                                                                               <  2:431702/250‑1 (MQ=255)                                                                        aTTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTGTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAg                                                                                                                                                                         <  1:337207/251‑1 (MQ=255)                                                                                          aTGGGATATAGCCTGAGGGGACTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGACTGGTTCGAACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTTGGAAGAGCACCGGAGGGGGCGTTAAACGATACCGCAACGTTGCTGGCGGGGCAGGGGCGTCACGATTTGCTCTCTGGGGGCCCAAAGCAAGGGCAACGGGCGGCGGCGTTTAATCGGCTaca                                                                                                                                                        >  2:461050/1‑250 (MQ=255)                                                                                              gATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATGTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATc                                                                                                                                                   <  1:627592/251‑1 (MQ=255)                                                                                                 ataGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGcc                                                                                                                                                <  2:130873/251‑1 (MQ=255)                                                                                                        gaagggCCTGTAATGCGTATCGGCAAGCGGTTGCTGCTGGGGTATTTTTTACTGGTGGCGTTGGCGGCCTGGTGCGTGCTGGCTAGTTTTGTCAAAGAGGGTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACATCGTCGCCCGTTTcg                                                                                                                                         <  2:529579/248‑1 (MQ=255)                                                                                                        gAGGGGCCTGTAATGCTTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcg                                                                                                                                         <  2:12238/251‑1 (MQ=255)                                                                                                        gAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGttt                                                                                                                                           <  1:346002/249‑1 (MQ=255)                                                                                                             gCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGcc                                                                                                                                     >  1:43555/1‑250 (MQ=255)                                                                                                                   aaTGCGTATCGGCATGCGTGTGCTGCTGGGCTATTTTTTACGGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCGCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCATCTACAACATCGCCATTTTCGCGCCAATATc                                                                                                                               <  2:445098/250‑1 (MQ=255)                                                                                                                              gCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAAc                                                                                                                   >  1:9941/1‑251 (MQ=255)                                                                                                                               cATGCGGTTGTTGCTGGGCTATTTTTTTCTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATGTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATTCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCGGGGGACCCAACGCTTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCTATTCGCGCCAATATCTGTGTCATTAACa                                                                                                                  <  2:329279/251‑1 (MQ=255)                                                                                                                                         ctgCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTGGCTATTTTTGTCAAAGGAGTTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGGGCGTCCCGATTTGCTCTCCGGGGACCCAACGCATGGGCAAATTGGGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTAGGGCCAATATCGGTGGCATTAACAAAGAGCGCa                                                                                                         >  2:515361/2‑250 (MQ=255)                                                                                                                                              gggCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGATGAGCAACGGACGGGACCTTAATCGGTACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGAACCACCGCCTGGGCAACTGGGGCAGGCGTTTAATCAGCCACAACATCGCCCGTTTCGCCCCACAATCGGAGGCATTCACAAAGTGCGCAATggct                                                                                                   >  1:291462/1‑248 (MQ=255)                                                                                                                                                  gaTTTTTTACTGGTGGCGGTGGCGGCCGGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGGTTAATCAGCTACAACATCGCCCGGTTCGCGCCAATATCGGTGGCGTTAACAAAGTGCGCAATGAATATCa                                                                                               <  2:89430/250‑1 (MQ=255)                                                                                                                                                   aTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAATGCGCAATGAATATCAt                                                                                              >  1:113635/1‑251 (MQ=255)                                                                                                                                                        ttACTGGGTGCGGTGGCGGCCTGGGTCGGACTGGCTATTTTTGGCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATGCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCGCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCGACAACATCGCCCGTTTCTCGCCAATATCGGTGGCATTAACAAATTGCGCAATGAATATCATGTCt                                                                                          <  2:560012/250‑1 (MQ=255)                                                                                                                                                         tACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCta                                                                                         <  1:573140/250‑1 (MQ=255)                                                                                                                                                          aCTGGTGGAGGTGGCGGCCTGGGTCGTACTGGCTATTTGTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCtat                                                                                        <  2:249341/250‑1 (MQ=255)                                                                                                                                                             ggtggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATg                                                                                     <  1:283526/250‑1 (MQ=255)                                                                                                                                                                ggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAcc                                                                                  <  1:150659/250‑1 (MQ=255)                                                                                                                                                                     tggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAttcg                                                                            >  1:20482/1‑248 (MQ=255)                                                                                                                                                                        cggCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAttc                                                                             <  2:247224/245‑3 (MQ=255)                                                                                                                                                                        cggCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAttc                                                                             >  1:247224/3‑245 (MQ=255)                                                                                                                                                                                  cGTACTGGCTATTTTTGTCAAAGAAGTTAATCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCTCAGGGTAAAGTAt                                                               <  2:279495/251‑1 (MQ=255)                                                                                                                                                                                   gTACTGGCTATTTTTGTCAAAGAAGTTACCCCGGGCGTGCGCAGAGCAACGGCGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGCGCTGGCGCATACCGATTTGCTCTCTGGGGACCCAACGCATGGGCAAATGGCGAAGGCGTTTAATCAGCCACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAgggagaagaatt                                                              >  2:607748/1‑242 (MQ=255)                                                                                                                                                                                     aCTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGGCCGATTCGCAGGGTAAAGTATTGt                                                            <  2:581836/251‑1 (MQ=255)                                                                                                                                                                                         gctatgtttGTCAAAGAAGTGAAACCGGGCGTGCGAAGATCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTACAGTATTGTTCGa                                                        <  2:304237/245‑1 (MQ=255)                                                                                                                                                                                            aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTc                                                     >  1:313329/1‑251 (MQ=255)                                                                                                                                                                                                ttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGc                                                  <  2:12961/250‑1 (MQ=255)                                                                                                                                                                                                     aaaGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTaaa                                                                   >  1:542722/1‑228 (MQ=255)                                                                                                                                                                                                     aaaGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTaaa                                                                   <  2:542722/228‑1 (MQ=255)                                                                                                                                                                                                       agaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCGCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAg                                          <  2:361029/251‑1 (MQ=255)                                                                                                                                                                                                         aagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATAATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGcc                                        >  1:26926/1‑251 (MQ=255)                                                                                                                                                                                                              aaaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGAATTGTTCGATTCGGCAAATAAAGCCGTTgg                                   >  2:111103/1‑251 (MQ=255)                                                                                                                                                                                                                 ccGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgggc                                 >  1:14005/1‑248 (MQ=255)                                                                                                                                                                                                                   gggCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAgg                              <  1:502871/251‑1 (MQ=255)                                                                                                                                                                                                                       gTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGatt                           >  1:38919/1‑250 (MQ=255)                                                                                                                                                                                                                          cGAAGAGCAACGGATGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTc                       <  2:510953/251‑1 (MQ=255)                                                                                                                                                                                                                           gAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTcg                      >  1:226743/1‑251 (MQ=255)                                                                                                                                                                                                                                 cAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTgg                <  2:130886/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCCACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATg            >  1:621300/1‑251 (MQ=255)                                                                                                                                                                                                                                      gAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGa           <  1:417533/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggggACGTTAAGCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACg         <  2:218253/251‑1 (MQ=255)                                                                                                                                                                                                                                               ttaaTCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCt  >  1:571099/4‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCT  >  NC_000913/4636581‑4637069

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