Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,610,440	C→T	E81K (GAA→AAA)	uxaB ←	tagaturonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,610,440	0	C	T	95.7%	150.8 / ‑2.1	47	E81K (GAA→AAA)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (22/23);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAATG  >  NC_000913/1610198‑1610684                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        cagagcACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGGCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGCTACCCGCTTAGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGCCCGCTTTCCCCTtct                                                                                                                                                                                                                                                 >  2:404053/5‑251 (MQ=255)    cgcACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCAt                                                                                                                                                                                                                                              <  1:391282/251‑1 (MQ=255)         cGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGCGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAgg                                                                                                                                                                                                                                         >  1:140840/1‑251 (MQ=255)               aaTGGCTAAAGCGCTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCGGATTTCACGATTAACCGAGCGAATCAGACGTGCGGCGCTGACCGCTTTCCCCTTCTTATTCAGTCCACgg                                                                                                                                                                                                                                   <  2:530151/251‑1 (MQ=255)                 tGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGAt                                                                                                                                                                                                                                 >  1:494870/1‑251 (MQ=255)                           gTTCGAACAGCAGACGGGGCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGgtggtg                                                                                                                                                                                                                       >  1:148516/1‑251 (MQ=255)                                     cagACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTAGCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCAt                                                                                                                                                                                                             <  1:379245/251‑1 (MQ=255)                                          gggTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGc                                                                                                                                                                                                        <  1:269446/251‑1 (MQ=255)                                           ggTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCg                                                                                                                                                                                                       >  1:376422/1‑251 (MQ=255)                                                gTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTc                                                                                                                                                                                                  >  2:250607/1‑251 (MQ=255)                                                     gCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTg                                                                                                                                                                                              >  1:365248/1‑250 (MQ=255)                                                                 aCCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATg                                                                                                                                                                                  >  1:9643/1‑250 (MQ=255)                                                                    gcgtgcgcGGCATAGAATTTATCGCCCGCATGGTAGCGTATACTCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGGGGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTGCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGTCCGCTTTCCCCTTCTCATACAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGt                                                                                                                                                                             <  2:302111/247‑1 (MQ=255)                                                                            gTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCg                                                                                                                                                                     <  2:376422/251‑1 (MQ=255)                                                                                        aTCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCAGCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGa                                                                                                                                                        <  2:194012/251‑1 (MQ=255)                                                                                         ccGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGAc                                                                                                                                                       <  2:473980/250‑1 (MQ=255)                                                                                           gCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGTCCATCCTGCGTGCTCAGTGACGGCGGGGATGAGAGtttaacatgaggtacaaacgac                                                                                                                                                        >  2:152900/1‑231 (MQ=255)                                                                                           gCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCa                                                                                                                                                     <  1:43539/251‑1 (MQ=255)                                                                                            cccGCATGGTAGCTGATACCCGCTTCGGTGGCGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCCCGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGttttaaacggcgacaaaacaacca                                                                                                                                                      >  2:371067/1‑229 (MQ=255)                                                                                               gCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGcc                                                                                                                                                 >  1:163208/1‑251 (MQ=255)                                                                                                   ggTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTGCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCGCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGaa                                                                                                                                             <  2:336/251‑1 (MQ=255)                                                                                                       gCTGATACCCGCTTCGGTGTTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCa                                                                                                                                         <  1:539418/251‑1 (MQ=255)                                                                                                        cTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTc                                                                                                                                          >  1:332298/1‑249 (MQ=255)                                                                                                                    tCGGGGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCCTCCTGCGTGCTCAGTGACGGCGGGGAAGA‑AGTTTCAATCGGGAGGAA‑AACGACCACGCCAGAATTCAAAACGGAGTGCTc                                                                                                                            >  2:347082/1‑251 (MQ=255)                                                                                                                        tggtgTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGG‑AGGTTCAATCGG‑ACGAACAACGACCCCGCCAGAATTCCGATCGGGGTGCTCAttt                                                                                                                        >  2:573653/1‑250 (MQ=255)                                                                                                                          gggtTAGAGAAAACAAATCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTTATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAg                                                                                                                      <  2:247630/249‑1 (MQ=255)                                                                                                                                    aaaCAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGGGTGCTCATTCAGGGGATCGAtt                                                                                                            >  1:344907/1‑251 (MQ=255)                                                                                                                                                 ttCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTcaaca                                                                                                <  2:148516/250‑1 (MQ=255)                                                                                                                                                       gTTGTGCGCCAGGTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGgcgc                                                                                         <  2:344907/251‑1 (MQ=255)                                                                                                                                                         tgtgCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGgcgcg                                                                                        <  2:365248/250‑1 (MQ=255)                                                                                                                                                            gcgtcagcttCAGGAATGCCTCGTCTGCTCTGTAGCCGCTGAGTGCACGATTAACCTAGCGAATCAGACGCGCGTCGCGGACCGCTTTCCCCTTCTCATTCAGGCCACGGATATTGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑Agg                                                                                   <  2:436746/243‑1 (MQ=255)                                                                                                                                                              gcCAGTTTCAGGAATTCATCGTATTCCCTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGCCCACGCCCGAATTCACATCGGTGTGATCCTTCAGGAGATCGCTTTGCCAGGCAACAACGGCGCGC‑CGGa                                                                                  >  2:95558/1‑250 (MQ=255)                                                                                                                                                                    ttCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGGTTCAATCGG‑ACGAACAACGACCACGCCCGAATTCAAGTCGGTGTGCTCATTCAGGAGATCGATTTTCCAGTCAACAAAGGC‑CGCCAGGAAGTTAcc                                                                           >  2:457483/1‑251 (MQ=255)                                                                                                                                                                         gcagtCATCGTAGGCACTGTGGGCGCTGATTGCACTATTAACCGAGCGAACCAGACGCGCGTCGCTGACCGCTTTCCCCTGCTCAGTCAGGCCACGGATAATGGTGGTGTACAGACCAGCCGGCGTGCTCAGTGACGGCGGGAATGA‑AGGTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGAGCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTTCCTTCAc                                                                      <  2:243420/247‑1 (MQ=255)                                                                                                                                                                                 cGTATTCACTGTATACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACTGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTg                                                               <  2:460827/250‑1 (MQ=255)                                                                                                                                                                                 cGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTg                                                               <  1:195962/250‑1 (MQ=255)                                                                                                                                                                                    aTTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATg                                                           >  1:94748/1‑251 (MQ=255)                                                                                                                                                                                        aCTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATAc                                                       >  1:473489/1‑251 (MQ=255)                                                                                                                                                                                                          aCGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATAGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTAACTTCACCAAACTGAATGATACGTTCTGGATACTGTGCAc                                     >  1:535024/1‑251 (MQ=255)                                                                                                                                                                                                                         aaTCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGcgac                      <  1:95558/251‑1 (MQ=255)                                                                                                                                                                                                                            cAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGCATGT‑AGTTTCCATCGG‑AAGAACAACGAACCCGCCCGAATTCCGATCCGTGTGCTCAATCAGGAGATCGCTTTTCCCGTCAACAAAGGCGCGC‑CGGGAGTTTCCTTCCACAACATTCATGGTACGTTATGGATACTATGCACCCGGACAAtaccgacc                     >  2:233845/1‑242 (MQ=255)                                                                                                                                                                                                                            cAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGAt                   <  1:573653/251‑1 (MQ=255)                                                                                                                                                                                                                               aCGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTa                <  1:250607/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAATg  <  1:457483/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAATg  <  2:494870/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATACTGCGTGCTCAGTGACGGCGGGAATGA‑AGGTTCAATCGG‑ACGAACAACGACCACGCCCGAATTTAGATCGGGGTGGTCCTTCAGGGGAGCGGCTTGCCCGTCAACAAAGGCGAGC‑AGGAAGTTAACTTCACCAAAATGAATGATACGTTCTGGATACTGTGCACCGGGGAAATCGCGGAGAATTAGTGTTTTCACAATg  >  1:220365/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccGCTTTCCCCTTCGCATTCAGGCCACGGATAATGGTGGTGTACAGACCAGCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAATg  <  2:290182/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        CAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGA‑AGTTTCAATCGG‑ACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGC‑AGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAATG  >  NC_000913/1610198‑1610684

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