BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,236	G→A	V34V (GTC→GTT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,236	0	G	A	96.1%	157.2 / ‑6.4	52	V34V (GTC→GTT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base G (1/0); major base A (20/29); minor base C (2/0); total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACGA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGG  >  NC_000913/67942‑68528

cTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTg                                                                                                                                                                                                                                                                                                 <  1:408239/300‑1 (MQ=255)

         gCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATg                                                                                                                                                                                                                                                                                        >  1:421967/1‑300 (MQ=255)

           cATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGtt                                                                                                                                                                                                                                                                                      <  1:363259/300‑1 (MQ=255)

              gtgCAGTCTGGTTCAGGTTCATAAAGTCCATATCGCTACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTggg                                                                                                                                                                                                                                                                                   >  1:741037/1‑300 (MQ=255)

                 cAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGCCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGTTCAGGCCGTTGATCCACATTTTGGCCGGGGAGANGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGa                                                                                                                                                                                                                                                                                <  2:176853/300‑1 (MQ=255)

                   gTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCCGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTACTTCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGTCGTCTCGCTTCCGGATTTAGCGCATTAACCC‑CGTGCTCGACATGTTGGGTGaac                                                                                                                                                                                                                                                                              >  2:632088/1‑298 (MQ=255)

                    tCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCt                                                                                                                                                                                                                                                                             <  1:71145/300‑1 (MQ=255)

                    tCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCt                                                                                                                                                                                                                                                                             <  2:566623/300‑1 (MQ=255)

                                     aaGTCCATAGCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCggg                                                                                                                                                                                                                                                            <  2:478456/300‑1 (MQ=255)

                                       gTCCATATCGATACTGTCCAACGACCGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAAGGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGTTGTGCAGCCACACCAACGGACCAGCGCAACGATCGTCGTAATTCACGTCGAAGCNAACAGCGGCGGCTTCATCCAGCGTGGTG‑CCTAGCGGTTTCCAACAG‑AGTTTACAGGGCAGTTTCGATTCCGTATTTAGCGCATTC‑CGAACGTGCTCGGTATGTTGCGTGCCCTGACACCTGGTTTCCGGGcc                                                                                                                                                                                                                                                          >  2:170519/1‑300 (MQ=255)

                                            tctcaaTACTGTCCCACGGCAGCGCCGCGTTGCACTGGGTGTGGACTTGCAGCCACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGACGGGGAGAAGGTGTGCAGCCACACCAANAGACCAGCGCAACGATCGTCGTAATTCGCGTAGCGGCAAACAGCGGTGATTTCATCCGGCGNGGTG‑CCCAGCGGTTTC‑CACACCAGTTTGTAGGGCAGTTTCGATTCCGTATTCAGCGCATTAACAA‑TGTGCTCGGCATGTTGGGTGACCTGCCGCAGAGTTTACGGGCCATACa                                                                                                                                                                                                                                                     >  2:595689/6‑300 (MQ=255)

                                            tctCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCACCGGTTTGTTGAGCATGGTCAGGCCGTTGATCCCCATTTTAGAAGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAAAGATCGTCGTACTTCGCGTCGCGGCAAATAGCGGTGATTTCGTCTGGTGTGGTGTCTCA‑CGGTTTC‑AATACCAGTTTGCAGGGCAGTTTCGGTTCCGTATGCGGCGCAGTAACAA‑CGTGCTCGGAATGTTGTGTGAcccgaagcaacgcttcctggccacaca                                                                                                                                                                                                                                                     >  2:397110/3‑275 (MQ=255)

                                            tctCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACNGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCAGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTA‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTTAGCGCANTAACAA‑CGTGCTCGGCATGTTGGGTGACATGACGCAGGGTTTCCGGGCCATACa                                                                                                                                                                                                                                                     >  2:689207/3‑300 (MQ=255)

                                              tCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTACAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGa                                                                                                                                                                                                                                                   <  1:139194/300‑1 (MQ=255)

                                                gATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATg                                                                                                                                                                                                                                                 <  1:423322/300‑1 (MQ=255)

                                                 aTACTGTCCCACGGCAGCGCCGCGCTGAACTGGGTGTGGAATTGCCGCAACGGTTTNTTGAGCATGGCCCGGCCGTTGCTCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCCCACCACCAGACCAGCGCACCGATTGTCGTCATTCGTGTCGCGGTAAATCGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAAGGCAGTTCCGCTTCCGTATTCAGCGCATTAACNA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATgc                                                                                                                                                                                                                                                >  2:683049/1‑300 (MQ=255)

                                                 aTACTGTCCAACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCCGCAAAGGCTTGTTGAGCATGGTAAGGCCGTTGATCCACCTTTTGGCCGGGGAGGAGGTGTGCAGCCACACCAACAGACCAGCGCAACGATCGTCGTAATTCGTGTCGCGGGAAATCGCGGTGATTTCATCCGGCGTGGTG‑CACAGCGGTTTC‑AACACCAGCTTGAAGGACAATTTAGCTTCCGTACTCAGCGCATTAAAAA‑AGCGCTCGGAATGTCGGGTGACCTGACGCAGGGTTCACGGGCCATACAGATgc                                                                                                                                                                                                                                                >  2:139357/1‑300 (MQ=255)

                                                     tGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATgctggc                                                                                                                                                                                                                                            <  1:462269/300‑1 (MQ=255)

                                                               cAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGAc                                                                                                                                                                                                                                  <  1:229663/300‑1 (MQ=255)

                                                                      gcgTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGACAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTAGCTGCCAATGAAAAACcac                                                                                                                                                                                                                           >  1:572439/1‑300 (MQ=255)

                                                                                gggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGTAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAAt                                                                                                                                                                                                                 <  1:416727/300‑1 (MQ=255)

                                                                                gggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAAt                                                                                                                                                                                                                 <  1:14244/300‑1 (MQ=255)

                                                                                                             aGCATGGTCAGGCCGTTGATCCCCATCTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGACGTGGTG‑CCCAGCGGCTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGACGGGTTTACGGGCCATACAGATGCTGGCTGCCAATGACAAACCACAGTTCATAATTATCACAACTCGTCATTATTGTGTACTta                                                                                                                                                                                    >  1:638774/1‑300 (MQ=255)

                                                                                                                aTGGTCAGGCCGTTGATCCACATTTTGGCCGGGTAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGNTTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATag                                                                                                                                                                                 <  2:652521/300‑1 (MQ=255)

                                                                                                                aTGGTCAGGCCGTTGATCCACATTTTGGCCGGGTAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATag                                                                                                                                                                                 <  1:483176/300‑1 (MQ=255)

                                                                                                                aTGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATag                                                                                                                                                                                 <  1:57914/300‑1 (MQ=255)

                                                                                                                 tGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTACGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGAACTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATaga                                                                                                                                                                                >  1:658028/1‑300 (MQ=255)

                                                                                                                  ggTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATagag                                                                                                                                                                               <  1:621618/300‑1 (MQ=255)

                                                                                                                   gTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGt                                                                                                                                                                              <  1:222888/300‑1 (MQ=255)

                                                                                                                                                           gCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGa                                                                                                                                      <  1:105469/300‑1 (MQ=255)

                                                                                                                                                                    ccaccaGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAgtgt                                                                                                                             <  1:454217/300‑1 (MQ=255)

                                                                                                                                                                          gaccagCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTAATAATTATCAAAAATCGTCATTATCGTGTCCTTCTCGAGTCGAAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCg                                                                                                                       >  1:701192/1‑300 (MQ=255)

                                                                                                                                                                               gcgcAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACCA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCAAATGACAAACCACACTTCATAATTATCACAAATCGTCGTTATCGTGTACTTATGGAGTCGCCACGGCATGGGCAGCCTGTGCAGGGGCGGAAGTTGGAAGATAGTGTCGTTCGGCGCt                                                                                                                  >  2:315858/1‑300 (MQ=255)

                                                                                                                                                                                       aTCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGccc                                                                                                          <  1:25653/300‑1 (MQ=255)

                                                                                                                                                                                          gtcgtAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATNATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCATGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCAtt                                                                                                       >  2:659216/1‑300 (MQ=255)

                                                                                                                                                                                              tAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAACGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGGATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAAATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGAAACGGCCTGGGCAGCCTGTGCCGGGGCGGACGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTg                                                                                                   >  1:276488/1‑300 (MQ=255)

                                                                                                                                                                                                   cgcgtcgcgGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAgc                                                                                              <  1:693891/300‑1 (MQ=255)

                                                                                                                                                                                                                        gATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCNGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGc                                                                         <  1:528902/300‑1 (MQ=255)

                                                                                                                                                                                                                        gATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGc                                                                         >  1:273022/1‑300 (MQ=255)

                                                                                                                                                                                                                                       ggtg‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGc                                                          <  1:176089/300‑1 (MQ=255)

                                                                                                                                                                                                                                        gtg‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCt                                                         >  2:127191/1‑300 (MQ=255)

                                                                                                                                                                                                                                                gCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGt                                                  <  1:169692/300‑1 (MQ=255)

                                                                                                                                                                                                                                                     ttC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAg                                             >  2:126944/1‑300 (MQ=255)

                                                                                                                                                                                                                                                              cAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTtct                                     <  1:401487/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                gTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTctct                                   <  1:307386/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                  ttGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTTGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTAc                                 <  1:383426/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                  ttGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCAGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTAc                                 >  2:237040/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                       gggCAGTTTCGCTTCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCAc                            <  1:687911/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                                    tCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCATGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTg               >  2:99506/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                                    tCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGACCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTg               >  2:667412/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                                    tCCGTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCATTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCAGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCATGTTCGCTGCGCGGTTGCAGGGTTTTCTCTAACGCACTGGCAATTTTTTg               >  2:66409/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                                       gTATTCAGCGCATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTgctg            >  1:731940/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                                                 cATTAACAA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATggg  <  1:735128/300‑1 (MQ=255)

CTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTG‑CCCAGCGGTTTC‑AACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACGA‑CGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGG  >  NC_000913/67942‑68528

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