BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,781	G→A	9.9%	E151E (GAG→GAA)	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,781	0	G	A	9.9%	740.7 / 53.4	261	E151E (GAG→GAA)	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base G (111/124); new base A (15/11); total (126/135)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

GCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGC  >  NC_002947/549632‑549930

gcgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgag                                                                                                                                                       <  1:6738674/150‑1 (MQ=39)

 cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                                                                                                      <  2:8984650/150‑1 (MQ=39)

 cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                                                                                                      <  1:5101217/150‑1 (MQ=39)

 cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                                                                                                      <  2:456249/150‑1 (MQ=39)

 cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                                                                                                      <  2:186056/150‑1 (MQ=39)

  gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                                                                                                                                     <  2:6506456/150‑1 (MQ=39)

   aTCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATg                                                                                                                                                    <  2:7888932/150‑1 (MQ=39)

       tGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAg                                                                                                                                                <  2:2336545/150‑1 (MQ=39)

       tGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAg                                                                                                                                                >  1:2204668/1‑150 (MQ=39)

       tGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAg                                                                                                                                                >  2:2228269/1‑150 (MQ=39)

       tGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAg                                                                                                                                                >  1:479581/1‑150 (MQ=39)

        ggttgtttCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGt                                                                                                                                               <  1:8259573/143‑1 (MQ=255)

            tGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                                                           <  2:2081389/150‑1 (MQ=39)

            gtttCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTcgc                                                                                                                                           <  1:552661/147‑1 (MQ=255)

              cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                                                                                                                                         >  2:7488576/1‑150 (MQ=39)

               tCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAc                                                                                                                                        >  2:72628/1‑150 (MQ=39)

               tCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAc                                                                                                                                        >  2:1299436/1‑150 (MQ=39)

               tCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAc                                                                                                                                        >  2:4609549/1‑150 (MQ=39)

               tCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAc                                                                                                                                        >  1:1637602/1‑150 (MQ=39)

                  gccgccGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctg                                                                                                                                     >  2:5035433/1‑150 (MQ=39)

                  gccgccGATGCTCCGAAGCCGCAAACCCGTGAGCACATCCAGCTGTCCCGACAGGTAGCCGTTCCGTACAACGTGGACTTCCACAACAAGGCCGACCAGGTAGACGACGCTCAGCAGCAGGAACAGGTCGAGAACGAAGTACGCGACctc                                                                                                                                     >  2:1983624/1‑149 (MQ=255)

                   ccgcagACGGCCCGATGCCGCAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTCGCGTTCCGCCAATCGACGTGTTCCTGAACAATTGCGACCAGGTTGACGATGCAGAACTGCTCGAACTGGTTGAACTGGAAGTACGCGAActtc                                                                                                                                    >  2:1514564/1‑148 (MQ=255)

                    cgccgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgct                                                                                                                                   <  2:2595851/150‑1 (MQ=39)

                     gccgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgctg                                                                                                                                  <  1:6713723/150‑1 (MQ=39)

                       tgacggcccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTc                                                                                                                                >  2:629323/8‑150 (MQ=255)

                        gacggcccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTcg                                                                                                                               <  1:2356179/144‑2 (MQ=255)

                            gcccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTcgaagt                                                                                                                           <  2:6855597/148‑6 (MQ=255)

                              ccGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                                                         <  2:2210104/150‑1 (MQ=39)

                              ccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTcgaagtac                                                                                                                         >  2:3621520/1‑143 (MQ=255)

                                gATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTACCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGa                                                                                                                       >  2:2801383/1‑150 (MQ=255)

                                    ccacagaCCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTc                                                                                                                   <  2:3813019/144‑1 (MQ=255)

                                    ccGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTc                                                                                                                   <  1:4030976/150‑1 (MQ=39)

                                    ccGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTc                                                                                                                   <  2:7536580/150‑1 (MQ=39)

                                      gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                 >  2:507426/1‑150 (MQ=38)

                                      acagaCCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTccc                                                                                                                 >  1:8759981/5‑150 (MQ=255)

                                       cAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCa                                                                                                                >  1:7175956/1‑150 (MQ=37)

                                       cAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCa                                                                                                                <  2:758940/150‑1 (MQ=37)

                                          aCCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCTTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGc                                                                                                             <  2:5722081/150‑1 (MQ=37)

                                          aCCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGc                                                                                                             <  2:9157025/150‑1 (MQ=37)

                                               tGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAGGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGcgacga                                                                                                        >  1:6169459/1‑150 (MQ=37)

                                               tGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGcgacga                                                                                                        <  2:5407390/150‑1 (MQ=37)

                                                   cacaTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACt                                                                                                    >  1:6625894/1‑150 (MQ=37)

                                                    acaTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTc                                                                                                   <  1:2110848/150‑1 (MQ=37)

                                                     caTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcc                                                                                                  >  1:8670112/1‑144 (MQ=2)

                                                      aTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCg                                                                                                 >  2:2018302/1‑150 (MQ=37)

                                                      aTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCg                                                                                                 >  2:2772387/1‑150 (MQ=37)

                                                      aTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGTGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcca                                                                                                 >  1:1402155/1‑143 (MQ=2)

                                                      aTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcca                                                                                                 >  1:3073296/1‑143 (MQ=2)

                                                         ctgctgTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatc                                                                                              <  1:9083041/150‑1 (MQ=37)

                                                         ctgctgTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatc                                                                                              <  2:2494018/150‑1 (MQ=37)

                                                          tgctgTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGaca                                                                                                      >  2:6151655/1‑141 (MQ=37)

                                                          tgctgTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatca                                                                                             >  1:4999398/1‑150 (MQ=37)

                                                          tgctgTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatca                                                                                             >  1:706768/1‑150 (MQ=37)

                                                             tcgctcGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGCGACGAAGTTCCAatcatca                                                                                          >  2:7350555/6‑150 (MQ=2)

                                                               tCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatcatc                                                                                        <  1:2440496/150‑1 (MQ=34)

                                                                  cGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGt                                                                                     >  1:8721388/1‑150 (MQ=34)

                                                                    tCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTc                                                                                   <  1:8346558/150‑1 (MQ=34)

                                                                         tAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTc                                                                              <  1:6260992/150‑1 (MQ=34)

                                                                         tAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTc                                                                              >  2:827239/1‑150 (MQ=34)

                                                                             cGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTAt                                                                          >  1:6458716/1‑150 (MQ=34)

                                                                                  cGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGccc                                                                     <  1:7832265/150‑1 (MQ=35)

                                                                                   gTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                                                                    >  1:7424683/1‑150 (MQ=37)

                                                                                       aTCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGaa                                                                <  1:8506456/150‑1 (MQ=37)

                                                                                                ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgac                                                       >  2:5936212/1‑150 (MQ=37)

                                                                                                  ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     <  2:4215360/150‑1 (MQ=37)

                                                                                                  ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     <  2:4823798/150‑1 (MQ=37)

                                                                                                  ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     >  2:2904226/1‑150 (MQ=37)

                                                                                                  ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     >  1:2992491/1‑150 (MQ=37)

                                                                                                  ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     <  1:3574197/150‑1 (MQ=37)

                                                                                                    tGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACa                                                   >  1:376227/1‑150 (MQ=37)

                                                                                                      aaCAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAAc                                                 <  2:7714569/150‑1 (MQ=37)

                                                                                                      aaCAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAAc                                                 >  1:2296345/1‑150 (MQ=37)

                                                                                                        cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGGCCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa                                               >  2:703846/1‑150 (MQ=37)

                                                                                                        cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa                                               >  1:5060196/1‑150 (MQ=37)

                                                                                                        cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa                                               >  1:5354078/1‑150 (MQ=37)

                                                                                                            gCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATg                                           <  2:4081652/150‑1 (MQ=37)

                                                                                                            gCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATg                                           >  1:4328804/1‑150 (MQ=37)

                                                                                                             ccGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                          <  2:1044317/150‑1 (MQ=35)

                                                                                                                  cTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACt                                     <  1:3410639/150‑1 (MQ=34)

                                                                                                                     gTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTAcc                                  <  2:739271/150‑1 (MQ=34)

                                                                                                                        gacgacGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCt                               <  2:948713/150‑1 (MQ=34)

                                                                                                                        gacgacGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCt                               <  1:4057228/150‑1 (MQ=34)

                                                                                                                          cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                             <  2:3076881/150‑1 (MQ=34)

                                                                                                                          cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                             <  2:5260331/150‑1 (MQ=34)

                                                                                                                            acgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCa                           <  2:3310642/150‑1 (MQ=35)

                                                                                                                             cgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaa                          <  2:2708699/150‑1 (MQ=37)

                                                                                                                              gCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaag                         >  1:2731229/1‑150 (MQ=37)

                                                                                                                                    ctgctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTg                   >  2:7565001/1‑150 (MQ=37)

                                                                                                                                    ctgctgGAACTGGTCGAGAGGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTg                   <  1:7264257/150‑1 (MQ=37)

                                                                                                                                      gctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGt                 <  1:4609549/150‑1 (MQ=37)

                                                                                                                                          gAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGaa             >  1:6055843/1‑150 (MQ=37)

                                                                                                                                                gTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTg       <  2:1308842/150‑1 (MQ=34)

                                                                                                                                                gTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTg       >  1:5714961/1‑150 (MQ=34)

                                                                                                                                                gTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTg       <  2:1684167/150‑1 (MQ=34)

                                                                                                                                                 tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg      <  2:2642609/150‑1 (MQ=34)

                                                                                                                                                 tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg      >  2:3424747/1‑150 (MQ=34)

                                                                                                                                                   gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt    >  1:8616984/1‑150 (MQ=35)

                                                                                                                                                   gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt    >  1:1331855/1‑150 (MQ=35)

                                                                                                                                                   gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt    >  2:3700928/1‑150 (MQ=35)

                                                                                                                                                     gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc  <  1:5035433/150‑1 (MQ=37)

GCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGC  >  NC_002947/549632‑549930

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