Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,784	A→T	intergenic (‑234/+51)	araD ← / ← araA	L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,784	0	A	T	100.0%	105.4 / NA	38	intergenic (‑234/+51)	araD/araA	L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/20);  total (18/20)

GTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑CAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACG  >  NC_000913/66551‑67016                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                         gTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCg                                                                                                                                                                                                                         >  1:286185/1‑251 (MQ=255)      cTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGc                                                                                                                                                                                                                     >  1:2163/1‑250 (MQ=255)        ccTTCGTGGCGGATGCGCTGTGATTAGCGGGCTAACAAAATCGCAGCGAGAAGGCCTGCGCAGACTCTCCAGCGTCGAATCAGGCGGTGAATGGCGGATGCGCTTTACGTATCCGGCCTTCAAGATCG‑GAACGGGTAGGCTTGCTAATACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGAGAAGACGCTGGAGCGTCGCATc                                                                                                                                                                                                                  <  2:313906/251‑1 (MQ=9)          ttCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATc                                                                                                                                                                                                                  >  2:43541/1‑249 (MQ=255)            cGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGc                                                                                                                                                                                                              <  1:266164/251‑1 (MQ=255)                 cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTAc                                                                                                                                                                                                         >  2:242754/1‑251 (MQ=255)                   gATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                                                                                                                       <  2:28231/251‑2 (MQ=255)                                     ggCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACtt                                                                                                                                                                                      >  2:241561/1‑250 (MQ=255)                                      gCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAg                                                                                                                                                                                    <  1:248199/251‑1 (MQ=255)                                                                 gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGt                                                                                                                                                           <  1:99404/249‑1 (MQ=18)                                                                  aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGGTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTcc                                                                                                                                                        >  1:240240/1‑251 (MQ=18)                                                                    aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcag                                                                                                                                                      <  1:138992/251‑1 (MQ=18)                                                                          cgcCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATGCACTTCGTTCcagcgcagc                                                                                                                                                <  2:181840/251‑1 (MQ=18)                                                                                        cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGGCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGCAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                                                                                                                  >  2:245600/1‑251 (MQ=18)                                                                                               tGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAACCGCTGGccgg                                                                                                                           >  2:138804/1‑248 (MQ=18)                                                                                               tGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCat                                                                                                                            <  2:170008/250‑2 (MQ=18)                                                                                                    gCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgtcgagt                                                                                                                      <  2:51700/251‑7 (MQ=18)                                                                                                                     tATCCGGCATACAAAATCGGCA‑CGCGTAGGCCTGCTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATCCGTACTACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGCTAACAATCAcc                                                                                                      >  2:116385/1‑250 (MQ=18)                                                                                                                      aTCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGt                                                                                                    <  1:241561/251‑1 (MQ=21)                                                                                                                          ggCCTACAAAATCG‑GAGCGCGTAGGCCTTATAAGACGCGTCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTGGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTTATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCTACTACTTAGCGACGAAATCCGTAATACACTTCGGTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGTGGGTATCGTTATCAATCACCTTGAtt                                                                                                <  2:282992/251‑1 (MQ=18)                                                                                                                                 caaagtcggactgcgTAGGCCTGATAAGACGTGCCTGTGTACCATCAGGCTGTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACGCTTCGTTCCAGCGCATCGTGTCTTTAAACGCTGGCATGCGGGTATCGTTATCAATCACCGTGATTTCAATg                                                                                          <  2:190500/239‑1 (MQ=21)                                                                                                                                          aacgcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCAtc                                                                                 <  2:176297/249‑1 (MQ=25)                                                                                                                                               gTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGc                                                                            >  1:281273/1‑251 (MQ=25)                                                                                                                                                       gATAAGACGCGCCTGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCGACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTAGCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGaactg                                                                      <  2:210989/249‑4 (MQ=25)                                                                                                                                                       gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGaactg                                                                      <  2:281273/249‑4 (MQ=25)                                                                                                                                                            gACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCata                                                               <  1:143835/251‑1 (MQ=25)                                                                                                                                                                cgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGGTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCgg                                                           >  1:104158/1‑250 (MQ=25)                                                                                                                                                                   cAGCTTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCCACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTACATCTCGGCGAACTGGCGCATATCGTTGa                                                        <  2:85971/251‑1 (MQ=25)                                                                                                                                                                                   gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGGGATTTTAATGTCGTGCATCCCGGCGAACCGGCGCCTATCGTTTAGGTTTAGCGCATGGcc                                        >  2:175021/1‑250 (MQ=25)                                                                                                                                                                                          gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACGGGCGCATATCGTTGAGGGTCAGCGCATGGCTGAAGAc                                  <  2:180317/250‑1 (MQ=37)                                                                                                                                                                                           ccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATCAATCACCGGGGTTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGCTGAGGTTCAGCGCCTGGCCGAAGCCgg                                >  2:284497/1‑251 (MQ=25)                                                                                                                                                                                                 gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGCTGAAGACGGTATggg                          >  1:200112/1‑250 (MQ=38)                                                                                                                                                                                                 gcgcTTTGCTTAACCGGCCTACACAATAGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAAGCTCGCAGGCGGGTATCGTTATCCATCACCGTGATTTCAATGTCGGGCATCTCGGCGAAACTGCGCATATCGTTGAGGTTCAGCGCACGGCTGAAGACGGTCTGGt                          >  2:39374/1‑251 (MQ=255)                                                                                                                                                                                                  cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCATCGCGTCGTTAAACGCTGGCAGGCGGGTATCGTTATCAATTACCGTGATTGCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTg                         <  1:299076/251‑1 (MQ=37)                                                                                                                                                                                                                  ccTACAAAATCTCATCGTGGAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCATGCGGGTATCGTTATCAATCGCTGTGATTTCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTGCAGCGCATGGCTGAATACGGTATGGGGCGGGCAACCTGCGAgg         <  2:203387/251‑1 (MQ=255)                                                                                                                                                                                                                   cTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGGTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCggggg        >  2:229452/1‑247 (MQ=255)                                                                                                                                                                                                                     aCAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATc      >  1:95748/1‑251 (MQ=255)                                                                                                                                                                                                                         aaTCGCAGCGTGTAGGCCTGATAAGACGCGTCCGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTTTCGTTATAAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCATATCGTTGAGGGTCAGCGCATGGCTGAAGAAGGGATGGTGCGCGCCACCAGAGAGGATACACg  >  2:54290/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                         GTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCG‑CAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACG  >  NC_000913/66551‑67016

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