Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,530,252	C→A	intergenic (‑9/+3)	sgcC ← / ← sgcB	putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,530,252	0	C	A	96.2%	156.0 / ‑3.1	52	intergenic (‑9/+3)	sgcC/sgcB	putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (27/23);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

GCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530008‑4530498                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gcgcAAACGTCAAAGACATGCGGGGCTAGGCCGAAACGCTAAATCATCACTTTTATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGCCCCATACCGACAAAACCAATCCCAACCTTTACCCCTCCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAtttt                                                                                                                                                                                                                                                    <  1:466754/249‑1 (MQ=255)          tCAAAGACATGCAGGGCGAGGCAGAAACGCTCAATCATCACTTGGATAGGGGGACTGAGGCGGTCGACCGCCATGATGACCACCATACCCATCCCGAAAAAACCAATCCCAACCGTTACCCCTGCCTTAACCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGAAGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACAAAATGCCTCATTTTATTTAGTTaa                                                                                                                                                                                                                                          <  1:93369/250‑1 (MQ=255)          acaaaagcaggCAGGGGTAGGCCGAAACGCCCAAACATCCTTTTGATGGGGGGACTGAGGCGGTCGACCGCCATGCCGATCCCCAGTCCCAAACCGACAAAACCAATCCCACCCTTCACCCCGGCCTTAATCGCTTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACAAAAACGTT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAAc                                                                                                                                                                                                                                         <  2:181057/241‑1 (MQ=255)               gccATGCAGGGCTAGGCCGAAACGCTCAATCATCCCTTTTATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCACCCGTTACCCCTGCCTTAAGCGCCTGAAGCCAGGGGATGCGGAAGATAAGGCCAATCACGTCCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGc                                                                                                                                                                                                                                    <  1:365695/249‑1 (MQ=255)                          cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTCTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGCTttgtt                                                                                                                                                                                                                         >  1:197475/1‑251 (MQ=255)                             ggCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGATAACAGAGCCTTGATTTgttgtt                                                                                                                                                                                                                      >  2:296429/1‑251 (MQ=255)                              cccGAAACGCCCAACCTCCACTTTGATGGGGGGACTGGGGCTGGCGATCCCCATGACGATCCCCATCCCCAAACCGACAAAACCATCCCCACCCGTTACCCCCGCCTTAACCGCCTGCAGCCAGTGGAGGCTAAAGATAAGGCCACTAACGACCATGATTATAGGGACGAAATCGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGttt                                                                                                                                                                                                                     <  1:444980/250‑1 (MQ=255)                                cGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATATAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGtttt                                                                                                                                                                                                                    >  2:176633/1‑250 (MQ=255)                                          aaccatcaCTTTGATAGGGGGACTGAGGTTGTGGACCGCCAGGACGATCACCAGCCCCATACCGACAAAACCAACCCCACCCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATGTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGcgtcgtc                                                                                                                                                                                                         <  1:163351/248‑1 (MQ=255)                                          aaTCATCACTTTGATAGGGGGACTGGGGCTTCCGATCGCCAGGACACTCACCAGCCCCATACCGCAAAAACCAATCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGcgtcgtc                                                                                                                                                                                                         <  1:120469/251‑1 (MQ=255)                                                              gACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAggg                                                                                                                                                                                     >  1:416967/1‑251 (MQ=255)                                                                   aGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGCTCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATAAAACATAATGCCTCATTTTATTGTGTTAACCGAGCCTTGATTTGTTGTTTTAACGCATCGTCGTTGATCCCTGTGAGCAgggggg                                                                                                                                                                                  >  2:383269/1‑246 (MQ=255)                                                                   aGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑AACGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGCGTTAAACGGGCCTTGATTTGTTGTTTTATCGCGTCGTCGGTGATACCTGTGCGCAGGgcggcg                                                                                                                                                                                >  2:155022/1‑251 (MQ=255)                                                                   aGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCAAACCGTTACCCCTTCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGCGACAGGATGTAATCAAACAGAATGCCTCATTTTATTGAGGTAACAGCGCCTTGATTGATGGTTTTAACGCGACGTAGTGTATCCCTGTGAACAGGgcggcg                                                                                                                                                                                >  2:383267/1‑251 (MQ=255)                                                                                   taaCGATCACCATTCCCATGCCGACAAAACCAATTCCAAGCGTTACCCCTCCCTTAATCGCCTGCAGCCCGGGGCTGGGGACAATACGGCCAAACACGCTCGTGATTTTAGGGACGAAGCCGGT‑ACCGCCGAGCGACGGGAAGTCGTCAAACATTATGCCTCATTTTATTGAGATAACAGAGCCGTGCTTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCAATTAACAGTAggg                                                                                                                                                                 >  2:423493/3‑250 (MQ=255)                                                                                                tCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACt                                                                                                                                                   >  1:294003/1‑251 (MQ=255)                                                                                                    aTACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTg                                                                                                                                               >  1:29228/1‑251 (MQ=255)                                                                                                      cccaacaaccccAATCCCAACCTTTACCCCTGACTTAATCGCCTTCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGTCCATGATGATAGGTACGAATAGGTT‑TCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCTATTAAGCGTACGTCTGCCGTACTCTCTGtttgt                                                                                                                                             <  1:383267/241‑4 (MQ=255)                                                                                                                ccAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATgg                                                                                                                                     >  1:78394/1‑249 (MQ=255)                                                                                                                 cAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGCTCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAgg                                                                                                                                  >  1:226070/1‑251 (MQ=255)                                                                                                                  aaTCCCAACCGTTACCCCTGACTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATATGGACGATAACGGT‑ACCGCCGAGAGACAGGATGTAAGAAAACCTAATGCATCATTTTATTGAGATAACAGAGCCTTGATTTGTTGTTTTAACGCGTCTGCGTTGATCCCTTTGAGCAGGGCGGCGCCATTAAGCGTAAGGAGTCCGTAATCACTGATGGGTCGCATgtgggt                                                                                                                                 >  2:97830/1‑246 (MQ=255)                                                                                                                    tCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGACTTGATTTGTTGTTTTAACGCGTCGTCGTTGGTCCCTGTGCGCAGGGCGGCGCCCTTAAGCGTAGGGATGCCGTAACCCCTGTTGGTTCGCCTGGTGGtt                                                                                                                                >  2:223131/1‑250 (MQ=255)                                                                                                                    ccccAACCGTTCCCCCGGCCTAAACCGCCTGCAGCCAGGGGATGCGGAAGAAAAGGCAAATCACGATCATGATGATAGGGACGAAAACGGT‑TCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTa                                                                                                                               <  1:322647/250‑1 (MQ=255)                                                                                                                        aaCCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAAAATAATGCCTCATTTTATTGAGTTAACAGAGCATTGATTTGTTGTTTTAACGCGTCGTCGTTGATACCTGTGAGCCTGGCGGCGCCATTACGCGTAGGGATGCCGTAATCCCTGTTGGTTCGCATGGAGGTTACAAt                                                                                                                           >  2:269946/1‑251 (MQ=255)                                                                                                                             ttACCCCGGCCTTAATCGCCTGCAGCCAGGGGCTGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGACTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCTGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGGTCGCATGGAGGTTACAATAAgg                                                                                                                       >  2:50048/1‑250 (MQ=255)                                                                                                                              tACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGt                                                                                                                      <  1:296429/250‑1 (MQ=255)                                                                                                                               cccccgcccTTAATCGCCGGCAGCCAGGGGAGGCGAAAGATAAGGCCAATCACGATCTTGATGATAGGCACCAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTcc                                                                                                                    <  2:267400/244‑1 (MQ=255)                                                                                                                                  ccTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCAt                                                                                                                  >  1:134682/1‑250 (MQ=255)                                                                                                                                        ttAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCACTCACGATCATGATGATAGGGACGACAACGGT‑ACCGCCGAGAGACAGGATGTAATCACACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGAGCAATTAAGCGTAGGGTTGCCGTAATCCCTTTGGGTTCGCATGGAGGTTACAATACGGTCCAAACCGttc                                                                                                           >  2:250541/1‑250 (MQ=255)                                                                                                                                        ttAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGtt                                                                                                            >  2:120757/1‑250 (MQ=255)                                                                                                                                          aaccGCCTGCAGCCAGGGGATGCGGAAGAAAAGGCCATACACGATCTTGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAACCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACa                                                                                                         <  1:422429/248‑1 (MQ=255)                                                                                                                                          aaTCGCCTGCAGCAAGGGGATGCGGAAGATAAGGCCAAACCCGATCATGAAGATAGGGACGAAA‑CGGTTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACa                                                                                                         <  2:170433/251‑1 (MQ=255)                                                                                                                                                        aGGGGATGCGGAAGATAAGGCCAATCACGATCATGAGGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAtct                                                                                           <  1:120757/251‑1 (MQ=255)                                                                                                                                                           ggATGCGGAAGAAAAGGCCAAACCCGATCATGATGATAGGGACAAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCAt                                                                                        <  1:376304/251‑1 (MQ=255)                                                                                                                                                               gcagAAGATAGGGCAAATCACGATCATGTTGATAGGGACAAAAACGGT‑TCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGGTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAg                                                                                    <  2:173783/248‑1 (MQ=255)                                                                                                                                                                    aGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGc                                                                               >  1:81393/1‑251 (MQ=255)                                                                                                                                                                                   cGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGAGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGGCCATACCGTTACAGTTAGGAGGGATATAATTCAGACAGCAATGGGTGGTAGTgg                                                                >  1:444317/1‑251 (MQ=255)                                                                                                                                                                                                      gAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTAGCATGGAGGTTACAATAAGGTCCGTACCGTTACCGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGGGGATGAAATGCCATGTtcagt                                             >  2:24829/1‑251 (MQ=255)                                                                                                                                                                                                                  cgccgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCtt                                  <  1:370662/251‑1 (MQ=255)                                                                                                                                                                                                                     cgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTCTTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtcag                                          >  1:321605/1‑240 (MQ=255)                                                                                                                                                                                                                     cgAGAGACAAGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTTATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTGCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtcag                                          <  2:321605/240‑1 (MQ=255)                                                                                                                                                                                                                      gagagaCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGAGTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGGAATCACTGTGGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAg                              <  2:304858/251‑1 (MQ=255)                                                                                                                                                                                                                            cAGGATGTAATCAAACATAATTCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTGTAACGCGTCGTCGTTGATCCCTGTGATCAGGGCGGCGCCATTAAGCGTAGTGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtg                        <  1:205523/251‑1 (MQ=255)                                                                                                                                                                                                                               gATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTACCGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGAGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATAACGTTACAGTTAAGAGGGATCTAATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGc                     >  1:210363/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttAATCAAACATAATGCCTCATTTTATTGAGTTAACAGATCCTTTATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGGAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCaa                   <  2:311102/249‑1 (MQ=255)                                                                                                                                                                                                                                    aaTCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCa                <  1:506282/251‑1 (MQ=255)                                                                                                                                                                                                                                      tCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATg              >  2:5961/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaaCATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGt            >  2:365272/1‑251 (MQ=255)                                                                                                                                                                                                                                         accATAATGCCTCATTTTATTGCGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCCGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCCCTGTTGGTTCGCATGGAGTTTACAATAAGGTCCATACCGTTACTGTTAAGCGGGATCTCATTCAGACAGCATTGGGCGGAAGTGGCTGAAATGCCCTGTTCAGTCAGGCACTCTTTCCGCTTGTGTGCAATcctgttt           >  2:240608/3‑245 (MQ=255)                                                                                                                                                                                                                                          aCATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGGCGTTGATCCCTGGGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGGAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGttt          <  1:423493/251‑2 (MQ=255)                                                                                                                                                                                                                                             tAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGa         >  2:160979/1‑249 (MQ=255)                                                                                                                                                                                                                                              aaTGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGt      >  2:373031/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cctcatTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAc  <  1:176633/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530008‑4530498

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