Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,819	T→C	intergenic (‑269/+16)	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,819	0	T	C	98.6%	192.7 / ‑5.5	71	intergenic (‑269/+16)	araD/araA	L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (32/38);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTC  >  NC_000913/66582‑67051                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              tagccgggcgaCAAAATCGCAGCGTGTAGGCCTGAATAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCga                                                                                                                                                                                                                                    <  2:8182/239‑2 (MQ=18)    ccgtccTATAAAATCGCGGAGTGTATGCCTGATAAGACGCGCCAGCGTCGCAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCGGCGTCGCATCAGGCGTTACACACCTGATGCGactact                                                                                                                                                                                                                               <  2:469276/247‑1 (MQ=17)       gtcTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAACTGATGCGACTACTTAg                                                                                                                                                                                                                            <  2:393852/249‑1 (MQ=18)        ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGc                                                                                                                                                                                                                           >  1:502055/1‑251 (MQ=255)           ccAAAACCGCATCGTGCAGGCCCGCTAAGACGCGCCAGCGTAGCATCAGGCGTTGAACGCCGGATGCGCTTTGCTTATCCGGCCTACCAAAGCGGAACGCGTAGGCC‑CGATAAGATGCGTTAGCGTCGCATCAGGCGTTGAATGTCGGATGCGCTTTGTTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGcga                                                                                                                                                                                                                         <  1:445710/249‑1 (MQ=21)                accGCAGCGGGTAGGCCTGAAAAGACGCGCCGGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                                                                                                                                                                   <  2:112382/249‑2 (MQ=18)                   cgagcttgtggtccggaaaaGACGCGCTATCGGCGCAACATGCGTAGAAGGCCGGGTTCGCTTTGCTTATCCGGCCTTCAAAAACGGACCGCGTAGGCC‑TGATAAGACGCGTCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCATTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCg                                                                                                                                                                                                                <  2:6456/233‑1 (MQ=16)                        gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATa                                                                                                                                                                                                           <  1:360681/251‑1 (MQ=18)                          gttggCCTGATAAGACGCGCTAGCGGCGCATCAGGCGTTGTATGCCGGTTGCGCTTTGCTTATACGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGGGTCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCTGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                                                                                                                                                          <  2:408807/247‑1 (MQ=18)                              gCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTc                                                                                                                                                                                                     <  1:145343/251‑1 (MQ=18)                                   aTAAGACGCGCCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACTCCTGATGCGACTACTGAGCGACGAAATCCGTAATACACTGCGTTcc                                                                                                                                                                                                <  2:640660/251‑1 (MQ=18)                                          gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTTTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcag                                                                                                                                                                                         <  2:429084/251‑1 (MQ=18)                                                      cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGGCTTTaa                                                                                                                                                                              >  1:555911/1‑250 (MQ=18)                                                       aTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAc                                                                                                                                                                            <  1:224452/251‑1 (MQ=18)                                                         cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAACTCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                                                                                                                                                          >  1:630660/1‑251 (MQ=18)                                                            gCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTgg                                                                                                                                                                       >  1:355268/1‑251 (MQ=18)                                                             cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGc                                                                                                                                                                      >  1:328889/1‑251 (MQ=18)                                                                   attccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                                                <  2:234386/248‑2 (MQ=18)                                                                   aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                                                <  1:467564/251‑2 (MQ=18)                                                                     tccgtattcgcTTTGCTTATCCGGCCTACAAAATCTGAACGCGTAGTCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGTCTACTTAGCGACGAAATCCGTAATACTCTTCTTTCCAGCGCAGCGCTTCTTTAAACGCTGGCAGGCggg                                                                                                                                                               <  1:609674/242‑3 (MQ=18)                                                                       cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATAAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggt                                                                                                                                                             >  2:512529/1‑246 (MQ=18)                                                                        ggATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatc                                                                                                                                                           <  1:305698/251‑7 (MQ=18)                                                                          aTGCGCGTTGCTTATCCGGCCTAAAAAATCGGAACGCGTATGCC‑GGTTAAGACGCGCATGCGTCGCTTCTGGCGTTTAATGCCGGATTCTCTTTGCTTTTCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCTTCATCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGTCAGGCGGGTA‑TCGt                                                                                                                                                        <  2:447001/251‑1 (MQ=19)                                                                            gcgcTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCG‑TCGGCCTGATAAGACGCGTCAGAGTCACATCAGGCGTTACACACCAGATGCGACTACTTAGCGACGAAATCCGAAATACACTACAGTCCAGCGCAGCGCGTCCTTAAACGCaggcaggcaggttgcgtt                                                                                                                                                        >  2:116528/1‑240 (MQ=17)                                                                               gttgtcttagcggggcTACAAAATCGGAACGCGTAGGCCGTG‑TAAGACGCGCCCGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCGTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTGAAACGCTGGCAGGCGGGTA‑TCGTTatca                                                                                                                                                   <  2:247892/236‑1 (MQ=18)                                                                                      tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGTG‑CAGGCCTGGTAAGACGCGTCCGCGTCGCATAAGGCGCTACACACCTGACGCGACGACTTCGCGTCGAAATCCGGAATACGCTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggactcgctacaacaccac                                                                                                                                             >  2:169781/1‑231 (MQ=12)                                                                                      tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTACTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CAc                                                                                                                                             <  1:234613/249‑1 (MQ=21)                                                                                       aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgggggggatagttatccatccccgg                                                                                                                                            >  2:24034/1‑228 (MQ=18)                                                                                       aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGt                                                                                                                                          <  1:24034/251‑1 (MQ=21)                                                                                          cGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTTTCAAT‑CACCGTggt                                                                                                                                       >  1:295393/1‑249 (MQ=21)                                                                                                         gaacgcgTAGGCC‑TGATAAGACGCGCCAGCGTCGCGTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAACATCGCAGCGTG‑TAGGCCTGATAAGACACGTCAGCGGCGCATCAGACGTTACACACCAGATGCGCCTACTTAGCAACGACATCCGTTATACACTTCGTTCCAGCACAGCGCGTCTTTAACAGCTGGCCGGCGGGGA‑TCGGTATCAAT‑CAccacgagtgtcatgttgtgagt                                                                                                                        >  2:432537/4‑231 (MQ=12)                                                                                                           acgcgTAGACC‑TGATAAGACGCGCCAGCGTCGCATGTGGCGTTGGATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGATAT‑CCTGCTAAGATGCGTCAGCGACGAATCATGCGTTACACACCTGATGCGAAAAATTATGGACGAACTCCGTCATACACTTCAGTCCAGCGCAGCGCGTCTTTAACCGCTTGCAGGAGGGTG‑GTGTTATACAT‑CACCGTggtgtctatgcagtgacatc                                                                                                                      >  2:73815/2‑232 (MQ=9)                                                                                                               gTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGc                                                                                                                 >  1:280445/1‑251 (MQ=25)                                                                                                                        gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGgc                                                                                                         <  1:485226/251‑1 (MQ=25)                                                                                                                                gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCg                                                                                                 <  1:203741/251‑1 (MQ=25)                                                                                                                                          cGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CAACCTGATTTTCATGTCGTA‑CATCTCGGCGacccggtgcatacacttgtggttc                                                                                        >  2:417907/1‑227 (MQ=21)                                                                                                                                           gaaTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGGCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTAGCAAT‑CACCGGGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAg                                                                                      <  2:295393/249‑1 (MQ=25)                                                                                                                                            cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgc                                                                                     <  1:512529/251‑2 (MQ=37)                                                                                                                                                aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCCTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTTGGGTTCAgcgcct                                                                                 >  1:173148/1‑246 (MQ=37)                                                                                                                                                 ggCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGGGATTTCAATGTCGGG‑CATCTCGGCGAACAGGCGCATATCGTTGAGGGTCAGCGCATg                                                                                >  1:637094/1‑251 (MQ=25)                                                                                                                                                 ggCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCTTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATg                                                                                <  1:699048/251‑1 (MQ=255)                                                                                                                                                     ttGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCG‑GAGATCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑AATCTTAGCGAACCGGCGCataacggttaggatcagcgccaggctg                                                                            >  1:380731/1‑227 (MQ=21)                                                                                                                                                             cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGGTCCGCGCATCGC‑T‑GAAGACggg                                                                  >  1:428156/1‑250 (MQ=37)                                                                                                                                                              ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACCGGCGCATATCGTTGGGGGTCCGCGCATGGC‑T‑GAAGACggg                                                                  >  1:71588/1‑249 (MQ=25)                                                                                                                                                                    gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTgc                                                           <  1:31632/251‑1 (MQ=38)                                                                                                                                                                          cTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCCATGTCGTG‑CATTTCGGCGAAATGGCGCATATCGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGccc                                                     >  1:234604/1‑250 (MQ=37)                                                                                                                                                                           ttATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGaactg                                                                                                           <  1:80130/198‑4 (MQ=37)                                                                                                                                                                             aTCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGGTTTTAAT‑CAACGTGTTTTCAATTTC‑TGTCATCTGGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCCTGGC‑TTGCAG‑CGGTCTGGTGCGCGccaccc                                                  >  2:298902/1‑250 (MQ=19)                                                                                                                                                                                cGGCCTACAAAATCGCAACGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATTAATACACGGT‑ATTTCAATGTCGTG‑CATCTTGGCGTACTCGCGCATATCGGTGGGGGTCAGCGCATGGC‑T‑GGAGACGGGATTGGGCGGGCCACCcgcg                                               >  2:347379/1‑251 (MQ=255)                                                                                                                                                                                  gCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACAAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAG‑TCGGTATCAAA‑TAACGAGATTTCATTATCGTT‑CATTTCGGCGAACTGGCCCCTATAGGTGAGGGTTACGGCATGGC‑T‑GAAGACgggagaggtggtgccacccgca                                               >  1:370894/1‑229 (MQ=255)                                                                                                                                                                                      aCAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTAGTTCCAGCGCAGCGCGTCTTTAAACGCTGGAAGGCGGGTA‑TCGTTATCCCA‑CACCGTGTTTTCCATATAGTG‑CATATCGGCGAATTGGCGCATATCGTTTGGGGTCAGCGCCTGTCAT‑GA‑GACGGTATGGGGGGGGCCACCCGCGAGGATc                                         >  2:534223/1‑251 (MQ=255)                                                                                                                                                                                         aaagcGCAGCGTG‑TAGGCCAGATAAGACGCGGCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATTGTGCGCGCCACCAGCGAGGATCCAc                                      <  2:301142/247‑1 (MQ=255)                                                                                                                                                                                         aaaacGCAGCGTG‑TAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCAc                                      <  2:328889/247‑1 (MQ=255)                                                                                                                                                                                         aaaTCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCAc                                      >  1:484964/1‑251 (MQ=255)                                                                                                                                                                                               cAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg                                <  2:355268/251‑1 (MQ=255)                                                                                                                                                                                                aGCGTG‑TAGGCCTGACAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCTGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGAACCACGCTTCgg                               >  1:432032/1‑251 (MQ=255)                                                                                                                                                                                                   gtg‑tAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GGAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGGTTCGGa                              >  1:478823/1‑249 (MQ=255)                                                                                                                                                                                                   gtg‑tAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGAATGGTGCGCGCCACCAGCCAGGATCCCCGCTTCGGa                              >  1:66449/1‑249 (MQ=255)                                                                                                                                                                                                        aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGCTTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGTGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGGGCGCCCCCAGCGGGGAGCCCCGCTTCGGAAGCAgg                        >  2:381247/1‑250 (MQ=255)                                                                                                                                                                                                         ggccttaTAAGACGCGGCAGCGTCGCATCAGTCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGGCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTGCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGtt                       <  2:238413/245‑1 (MQ=255)                                                                                                                                                                                                         ggCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcg                                                                                    <  1:586858/192‑3 (MQ=255)                                                                                                                                                                                                             tgATAAAACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAA‑TCGTTATCAAT‑CACCGTGGTTTCAATGCAGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCAGGA‑T‑TAAGACGGTAAGTTTCGCGCCCCCATCGAGGGTCCACGCTTCGGGAGCAGTTGGc                    >  1:127941/2‑250 (MQ=255)                                                                                                                                                                                                                tAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACTGGCGCATAACGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGGAGCAGTTGGCAGa                 >  1:296606/1‑250 (MQ=255)                                                                                                                                                                                                                 aaGACGCGTCAGCGTTGTATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGAt                <  2:663465/250‑1 (MQ=255)                                                                                                                                                                                                                 aaGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGAt                <  1:381247/250‑1 (MQ=255)                                                                                                                                                                                                                  aGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATcc              <  1:671835/251‑1 (MQ=255)                                                                                                                                                                                                                  aGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATcc              <  2:194814/251‑1 (MQ=255)                                                                                                                                                                                                                    aCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCgg            <  1:264377/251‑1 (MQ=255)                                                                                                                                                                                                                          cAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTgcg       <  1:290081/250‑1 (MQ=255)                                                                                                                                                                                                                             cGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑CATCTCGGCGAACTGGCGCCTATCGTTGAGGTTCAACGCATGGC‑T‑GAAGACGGGATGGGGCGCGCCACCAGGGCAGATCCACGCCTTGGAAGcgggttggagatacgggttcgcttt   >  2:454739/1‑227 (MQ=255)                                                                                                                                                                                                                              gTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTc  <  1:181078/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              TATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTC  >  NC_000913/66582‑67051

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