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	97.6%	261.6 / ‑2.4	83	V34V (GTC→GTT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (46/34); total (48/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACGAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGG  >  NC_000913/67995‑68479

ggTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTGGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGGCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTTGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCg                                                                                                                                                                                                                                            <  2:435423/250‑1 (MQ=255)

       cGGCAGCGCAGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCAGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTACCAAC‑GTGCTCGGCATGtt                                                                                                                                                                                                                                     >  1:689772/1‑251 (MQ=255)

        ggCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCACGGCAAATAGCGGTGATTTCATCCGGCGCGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTCCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTg                                                                                                                                                                                                                                    >  1:357528/1‑251 (MQ=255)

         gCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTgg                                                                                                                                                                                                                                   >  1:497216/1‑251 (MQ=255)

            gcgcCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTg                                                                                                                                                                                                                                <  1:520812/251‑1 (MQ=255)

              gccgcGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGAc                                                                                                                                                                                                                              >  1:162023/1‑251 (MQ=255)

                agctgtgcaacgggggtgGAATTGCCGCGATGGTGTGTTGAGCAGGGTCAGGCCTTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCCCACCACCAGACCAGCGCAACG‑ATCGTCGTAAATCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGCTTCAACACCAGTTTGTAGGGCAGTTTCGCTTCCGTTTACAGCGCATTAACATC‑GTGCGCGGCATGTTGGGTgcccg                                                                                                                                                                                                                            <  2:482550/236‑5 (MQ=255)

                         cTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAggg                                                                                                                                                                                                                   >  1:569067/1‑251 (MQ=255)

                           gggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACTCCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGtt                                                                                                                                                                                                                 <  2:551315/251‑1 (MQ=255)

                           gggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGtt                                                                                                                                                                                                                 <  1:146569/251‑1 (MQ=255)

                             gtgtGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTc                                                                                                                                                                                                               >  1:170003/1‑251 (MQ=255)

                               gtgGAATTGGAGCCACGGTTGGTTGAGCATGTTCAGGCCGTTGATCCACATTTTGGCCGGTGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCg                                                                                                                                                                                                             <  2:625651/251‑1 (MQ=255)

                                       gcagcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATa                                                                                                                                                                                                      <  1:171739/250‑1 (MQ=255)

                                                   tgttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATgctggctg                                                                                                                                                                                         >  1:46776/1‑251 (MQ=255)

                                                       gAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCa                                                                                                                                                                                      <  1:296808/250‑1 (MQ=255)

                                                            tGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGAc                                                                                                                                                                                 >  1:241788/1‑250 (MQ=255)

                                                             ggTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGGGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGCCGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACa                                                                                                                                                                                >  2:144392/1‑250 (MQ=255)

                                                                 aGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACca                                                                                                                                                                           <  1:54400/251‑1 (MQ=255)

                                                                 aGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGAGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAAcc                                                                                                                                                                            >  1:274646/1‑250 (MQ=255)

                                                                  ggCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACcac                                                                                                                                                                          <  1:596014/251‑1 (MQ=255)

                                                                    ccGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCCGACCCGCGCAACG‑ATCGTCGGAATTCGCGTCGCGGCCAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GAGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCacac                                                                                                                                                                        >  1:286590/1‑251 (MQ=255)

                                                                    ccGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAA‑GTGCTCGGCATGTTGGGTGACCTGACGCGGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGCCAAACCacac                                                                                                                                                                        >  2:569321/1‑251 (MQ=255)

                                                                     cGTTGATCCACATTTTGGCCGGGGGGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGCCGTAAGTCGCGTCGCGGCAAATAGCAGAGGTTTCATCTGACGGGTTTCCCAGCGGTTTCAAAACCAGTGTTCAGGGCAGTTTCGCTTCCGTATTAAGCGCATTAACAAC‑GTACTCGGCATGTTGGGTAACCTGCCGCAGGGCTTCCGGGCCATACAGTAGGTGGCTTCCCATGACAAACcaaaat                                                                                                                                                                       >  2:136695/1‑247 (MQ=255)

                                                                     cGTTGATCCACAGTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGGTATTCGCGTCGCGGCAAATAGCGGTGAGTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACt                                                                                                                                                                       <  2:241788/251‑1 (MQ=255)

                                                                        tGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTc                                                                                                                                                                     >  1:36149/1‑250 (MQ=255)

                                                                         gATCCCCATTTTGGCCTGGGAGAAGGCGGGCAGCCACAACATCAGACCAGCGCAACG‑AACGTCTTAAGTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTTCCCAGCGGTTTCAACACCAGTGTGCAGGGCAGTTTCGCTTCCGTATTCATCGCATTACCGAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTGTCCGGGCCATACAGATGCTGTCTGCCAATGACAAACCACACttta                                                                                                                                                                    <  2:493274/250‑3 (MQ=255)

                                                                          aTCCCCATTTTGGACGGGGATAAGGTGTGCAGCCTCCCCACACGACCAGCGCAACGCAT‑GTCGTAATTCGGGTCGCGGCAAATAGCGGTGAGTGCAACCGTCGTGGTACCCACCGGTTTCAACCCCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGGATTAACAAC‑GTGCTCGTCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGACTGTGACAATCCACGCTTCATa                                                                                                                                                                  <  2:280154/251‑1 (MQ=255)

                                                                          aTCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGCTTTACTCAGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTACCAAC‑GTGCTCGGCATGTTGGGTGACCTGGCGCAGGGTTTCCAGGCCATACAGCTGCTGGCTGCCAATGACAAACCACAATTCAt                                                                                                                                                                   >  2:204904/1‑250 (MQ=255)

                                                                           tCCACATTTTGGCCGGGGAGAAGGAGTGCAGCCACAACACCAGACCATCGCAACG‑ATCGTCGTAATTAGCGTCGCGGCAAATAGCGGAGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGAGCAGATTCGCTTCCGTATTCAGAGCATTCACAAC‑GTGCTCGGCATGTTGAGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTAGCTGCACATGACCAACCAAACTTcgta                                                                                                                                                                  >  2:376355/1‑247 (MQ=255)

                                                                            ccACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCAACAGACCAGCGCACCG‑ATCGTCGTAATGCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAAt                                                                                                                                                                <  2:162023/251‑1 (MQ=255)

                                                                               caTTTTGGCCGGGGAGAAGGTGTGCAGCCACATCACCAGACCATCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCAGTGATTTCATCCGGCGAGGTGCCCAGCGGTTTCCACACCAGTTTGCAGTGCAGTTTCGCTTCCGTATTAAGAGCATTAACAAC‑GTGCTCGGCATGTTAGGAGACCTGACGCAGGGTTTCCGGGCCCTACCGGTGCTGGCTTCAAATGACAAACCACCCTTCATacttgt                                                                                                                                                             >  2:129313/1‑246 (MQ=255)

                                                                                        cGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCg                                                                                                                                                    <  1:342914/251‑1 (MQ=255)

                                                                                           ggAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTc                                                                                                                                                  >  2:246424/1‑250 (MQ=255)

                                                                                              gaAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACCCTTCATAATTATCAAAAATCGTCAtt                                                                                                                                               >  1:478652/1‑250 (MQ=255)

                                                                                               aaGGTGTGCAGCCACACCACCAGCACAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGGTTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTAt                                                                                                                                             <  1:215187/251‑1 (MQ=255)

                                                                                                    gtgCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCgtgt                                                                                                                                        <  1:273393/251‑1 (MQ=255)

                                                                                                    gtgCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCgtgt                                                                                                                                        >  1:98827/1‑251 (MQ=255)

                                                                                                    gtgCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCgggt                                                                                                                                        >  2:211914/1‑248 (MQ=255)

                                                                                                                caccaGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAACCCACACTTCATAATTATCAAAAATCGTAATTATCGTGTACTTATAGAGTc                                                                                                                            >  2:547767/1‑251 (MQ=255)

                                                                                                                caccaGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTc                                                                                                                            >  1:333329/1‑251 (MQ=255)

                                                                                                                    agaccagAGCAACG‑ATCGTCGTAATTCGTGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCACAGCGGTTTCAACAACAGTTTGTAGGGCAGGTTCGCGGCCGGAGTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGTCATATTGATTCTGGCTGCCAATGACAAACCACACATCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCaa                                                                                                                        <  1:315748/251‑1 (MQ=255)

                                                                                                                     gaccagCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAAc                                                                                                                       >  1:163827/1‑251 (MQ=255)

                                                                                                                      accagCGCAACG‑ACCGTCGTCATTCGCGTCGCGGCAAATAGCTGTGATTTCATCCGGCGTGCTGCCCAGCGGTTTCAACACCATTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGTATGTTGGGTGTCCTGACGCAGGGTTTCCTGTCCATACAGATGCTGGCTGCCAATGACAACCCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACg                                                                                                                      <  1:521735/251‑1 (MQ=255)

                                                                                                                         agCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGcc                                                                                                                   >  1:587491/1‑251 (MQ=255)

                                                                                                                          gcgcAACG‑TTCGTCGTAATTCGCGGCGCGGCAAATAGCTGTGATTTCAGCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCGGGCTGCCAATGACAAACCACACTTCGTAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCt                                                                                                                  <  2:66116/251‑1 (MQ=255)

                                                                                                                             cAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTggg                                                                                                               >  2:703792/1‑251 (MQ=255)

                                                                                                                              aaCG‑CTCGTAGTAATTCACGTCGCGGCAAATACCGGAGCTTTCATCCGGCGTCGTGCCCCGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCAGTCTTCAGCGCATTAAC‑ACAGTGCTCGGCATGATGGGTGACCTGACGCAGGGTTTCCGAGCCATACAGATGCTGACTGCCACTGCCCAACCACACTTAATAATTATCACAAATCGTCATTATCGTGTCCTTATCGAGTCGCAACGACCTggtc                                                                                                              >  2:331448/1‑249 (MQ=255)

                                                                                                                              aaCG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGAGTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGATTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGc                                                                                                              >  1:234939/1‑251 (MQ=255)

                                                                                                                                     cgtcgtAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCt                                                                                                                                                                                                                            <  2:469294/135‑1 (MQ=255)

                                                                                                                                     cgtcgtAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCt                                                                                                                                                                                                                            >  1:469294/1‑135 (MQ=255)

                                                                                                                                     cgtcgtAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCaca                                                                                                                                                         >  1:264070/1‑200 (MQ=255)

                                                                                                                                         gtAAATCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGcc                                                                                                    <  2:163827/251‑1 (MQ=255)

                                                                                                                                          tAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGTTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACGAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCg                                                                                                   >  1:221641/1‑251 (MQ=255)

                                                                                                                                               cgcgtcgcgGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGc                                                                                               <  1:709493/250‑1 (MQ=255)

                                                                                                                                               cgcgtcgcgGCAAATAGCGGTGAGTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCg                                                                                              <  2:46776/251‑1 (MQ=255)

                                                                                                                                                     gcgGCAAATGGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCTTACAGATGCTGGCTGCCAATGCCAAACCACACTTCATAATTATAAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGACGGGGCGGCAGtt                                                                                        >  2:304688/1‑251 (MQ=255)

                                                                                                                                                         cAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGa                                                                                     >  1:375259/1‑250 (MQ=255)

                                                                                                                                                             tAGCGGTGATTTCATCCGGAGTGGTGCCCAGCGGTTTCAATACCAGTTTGCAGAGCAGTTTCGCTTCCGTATTCAGCGCCTTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCACTGACAAACCACACTTCATAATTCTCACAAATCGTCATTATCGTGTCCTTCTAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGGAGTTGGAAgct                                                                                 >  2:43683/1‑248 (MQ=255)

                                                                                                                                                                  gTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGTCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAACAATCGTCATTATCGTGTCCTTATAGAGTCGCAACAGCCTGGGCAGCCTGTGCCGGGGCGGAAGCTGGAAGGTAGtgtt                                                                           >  2:412314/1‑251 (MQ=255)

                                                                                                                                                                            ccGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCt                                                                 <  2:497216/251‑1 (MQ=255)

                                                                                                                                                                                  tggtgCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGc                                                           >  1:26631/1‑251 (MQ=255)

                                                                                                                                                                                       cccAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGCAGATAGTGTTGTTCGGCGCTCATCGCCCAtt                                                      >  1:414361/1‑251 (MQ=255)

                                                                                                                                                                                          cggggTTTCAACACCAGTTTGCAGGTCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCt                                                   <  2:526204/248‑1 (MQ=255)

                                                                                                                                                                                            cGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTg                                                  >  2:566047/1‑250 (MQ=255)

                                                                                                                                                                                             ggTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGAt                                                <  2:569067/251‑1 (MQ=255)

                                                                                                                                                                                                  cAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAgcgg                                           <  1:703792/251‑1 (MQ=255)

                                                                                                                                                                                                               cagggtattttCGCGTCCGTAGTCAGCGCATTAACAAC‑GGGCTCGGCATGTTTGGTGACCTGACGCAGGGTTTTCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGGCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGtt                              <  2:560672/243‑1 (MQ=255)

                                                                                                                                                                                                                         tCGCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGTTAGCGGCGATAAAGCTGTTCAAAGCGtt                     <  1:219572/250‑1 (MQ=255)

                                                                                                                                                                                                                           gCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTtgtg                  >  1:395521/1‑251 (MQ=255)

                                                                                                                                                                                                                           gCTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTtgtg                  <  1:302861/251‑1 (MQ=255)

                                                                                                                                                                                                                            cTTCCGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGc                 >  1:546609/1‑251 (MQ=255)

                                                                                                                                                                                                                              tCCGTATTCAGCGCATTAACAAC‑GTGATCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAACCCACACTTCATAATTATCAAAAATCGTCATTATCGTGTACTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGGAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTAAAAGCGTTGTGccc               >  2:456171/1‑250 (MQ=255)

                                                                                                                                                                                                                                cgttatCAGCGCATTAACAAC‑GTGCCTGGCATGTTGGGGGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGAGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGTGATAAAGCTGTTCAAAGCGTTGTGCCTgt             <  2:170003/246‑2 (MQ=255)

                                                                                                                                                                                                                                cgtcttCAGCACATTAATAAC‑GTGCTAGGCTTGTTGGGTGTCCTGACGCTTGGTGTCCGTGCCATACAGATTCTGGCTGCCAATTACAAACCACTCTTCATAATTATCAAAAATCGTCATTATCGTGTCCATATAGAGTCGCAACGGCCTGGGCAGCCTGTGTCGGGGCGTAAGTTGGAAGATAGTGTTGTTCTTCGCTCAGCTCCCATTGGTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTgt             <  1:396763/247‑2 (MQ=255)

                                                                                                                                                                                                                                cGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATc                                                                                                                                            <  1:692918/124‑1 (MQ=255)

                                                                                                                                                                                                                                cGTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGACGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGGTCAAAGCGGTGTGCCttt             >  2:266284/1‑249 (MQ=255)

                                                                                                                                                                                                                                cGTATTCAGAGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATc                                                                                                                                            >  2:692918/1‑124 (MQ=255)

                                                                                                                                                                                                                                 gTATTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGACAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAGGTTGGAAGATCGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGAGATAAAGCTGTTCAAAGCGTTGTGCCTgtt            >  2:74753/1‑249 (MQ=255)

                                                                                                                                                                                                                                   aTTCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGGGTTGTTCGGCGCTCATCGCCCATTGCTGAAAGCGGCGATAAAGCTGTTCAAAGCGtttggcctgttcg          >  2:227395/1‑240 (MQ=255)

                                                                                                                                                                                                                                    ttCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGa                                        <  2:525879/220‑1 (MQ=255)

                                                                                                                                                                                                                                    ttCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGa                                        >  1:525879/1‑220 (MQ=255)

                                                                                                                                                                                                                                    ttCAGCGCATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGc         >  1:263589/1‑251 (MQ=255)

                                                                                                                                                                                                                                           cATTAACAAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGATTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCgg  >  1:481816/1‑251 (MQ=255)

TGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACG‑ATCGTCGTAATTCGCGTCGCGGCAAATAGCGGTGATTTCATCCGGCGTGGTGCCCAGCGGTTTCAACACCAGTTTGCAGGGCAGTTTCGCTTCCGTATTCAGCGCATTAACGAC‑GTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGG  >  NC_000913/67995‑68479

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