BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,619	C→A	6.4%	A97A (GCC→GCA)	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,619	0	C	A	6.4%	488.4 / 17.1	157	A97A (GCC→GCA)	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base C (67/80); new base A (6/4); total (73/84)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 170 total aligned reads displayed.

GACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGC  >  NC_002947/549470‑549767

gACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGcc                                                                                                                                                      >  2:1075321/1‑150 (MQ=32)

gACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGcc                                                                                                                                                      >  1:3266669/1‑150 (MQ=32)

 aCAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGccc                                                                                                                                                     >  1:4707364/1‑150 (MQ=32)

    aGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGa                                                                                                                                                  <  2:316074/150‑1 (MQ=32)

         gACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGAc                                                                                                                                             <  2:3683585/150‑1 (MQ=32)

            tCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgc                                                                                                                                          >  2:4926937/1‑150 (MQ=32)

             cGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgcg                                                                                                                                         <  2:4141295/150‑1 (MQ=32)

                 tCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGAt                                                                                                                                     >  2:1750213/1‑150 (MQ=32)

                   cGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATcc                                                                                                                                   >  2:2250634/1‑150 (MQ=32)

                    ggAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCt                                                                                                                                  >  2:2608604/1‑150 (MQ=32)

                    ggAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCt                                                                                                                                  >  1:2814316/1‑150 (MQ=32)

                     gaagaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTg                                                                                                                                 <  1:2608604/150‑1 (MQ=32)

                      aagaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTgg                                                                                                                                >  2:4116283/1‑150 (MQ=32)

                       agaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGt                                                                                                                               <  2:4237662/150‑1 (MQ=32)

                       agaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGt                                                                                                                               <  2:4012890/150‑1 (MQ=32)

                           aaaGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGc                                                                                                                           >  2:3641710/1‑150 (MQ=32)

                           aaaGCGCGCGGTATCACCATCAACACCGCTCACCTATAGAACAACACGAACATTCTAAACAACGCTCACGATGACTGACCAGGTCACGCTGACTACGAGACGAACATTATCACCGGAGCAGCCCAGATGGACGGCGCGATCCAGGttgtc                                                                                                                           >  1:4574220/1‑147 (MQ=11)

                              gcgcgcGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCg                                                                                                                        >  1:1785132/1‑150 (MQ=32)

                                  gcgGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGgccg                                                                                                                    >  1:2042817/1‑150 (MQ=32)

                                   cgGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGgccgc                                                                                                                   <  1:3101198/150‑1 (MQ=17)

                                        tCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATg                                                                                                              <  2:932364/150‑1 (MQ=32)

                                         cACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATgg                                                                                                             <  1:3351160/150‑1 (MQ=17)

                                           ccATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTc                                                                                                           >  2:369224/1‑150 (MQ=32)

                                           ccATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTc                                                                                                           <  1:3745122/150‑1 (MQ=32)

                                           ccATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTc                                                                                                           >  1:3948116/1‑150 (MQ=32)

                                            cATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                                                                                                          >  1:580654/1‑150 (MQ=32)

                                            cATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                                                                                                          <  1:796781/150‑1 (MQ=32)

                                            cATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                                                                                                          <  2:2635039/150‑1 (MQ=32)

                                            cATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                                                                                                          <  2:1647533/150‑1 (MQ=32)

                                             aTCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCg                                                                                                         >  2:1234137/1‑150 (MQ=32)

                                                aaCACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATg                                                                                                      <  2:1785132/150‑1 (MQ=32)

                                                       cTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCaaa                                                                                               <  1:3700396/150‑1 (MQ=21)

                                                       cTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCaaa                                                                                               <  2:1999348/150‑1 (MQ=34)

                                                        tCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAAc                                                                                              >  1:4211650/1‑150 (MQ=34)

                                                        tCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAAc                                                                                              <  2:4885759/150‑1 (MQ=34)

                                                        tCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGGGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAAc                                                                                              >  2:1013498/1‑150 (MQ=21)

                                                          aCGTCGAGTACAACTCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAAccc                                                                                            <  1:2507653/150‑1 (MQ=21)

                                                          aCGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAAccc                                                                                            <  2:4392046/150‑1 (MQ=34)

                                                            gTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTCCGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGt                                                                                          <  2:4430541/150‑1 (MQ=21)

                                                            gTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGt                                                                                          <  2:2880934/150‑1 (MQ=34)

                                                              cGAGTACAACTCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCCCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                        <  1:1614810/150‑1 (MQ=21)

                                                              cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                        <  1:779861/150‑1 (MQ=34)

                                                              cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                        <  2:282916/150‑1 (MQ=34)

                                                              cGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                        <  1:4465545/150‑1 (MQ=34)

                                                                     aaCTCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATc                                                                                 <  1:1626968/150‑1 (MQ=21)

                                                                     aaCTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATc                                                                                 >  1:3270179/1‑150 (MQ=34)

                                                                       cTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCct                                                                               <  1:4756543/150‑1 (MQ=34)

                                                                        tCGAACATTCGTCACTACGCTCCCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctg                                                                              <  1:2377245/150‑1 (MQ=21)

                                                                        tCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctg                                                                              >  1:2233047/1‑150 (MQ=34)

                                                                          gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                            <  1:1381118/150‑1 (MQ=35)

                                                                          gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                            >  1:4192415/1‑150 (MQ=35)

                                                                          gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                            <  2:1101623/150‑1 (MQ=35)

                                                                          gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                            >  1:4623512/1‑150 (MQ=35)

                                                                          gAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCctgct                                                                            <  1:3985444/150‑1 (MQ=35)

                                                                                  gTCACTACGCTCACGTTGACTGCCCAGGTCACACTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTc                                                                    <  2:166022/150‑1 (MQ=21)

                                                                                    cACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAg                                                                  >  1:1218076/1‑150 (MQ=34)

                                                                                    cACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAg                                                                  >  1:1218019/1‑150 (MQ=34)

                                                                                       tACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTa                                                               >  1:1317540/1‑150 (MQ=34)

                                                                                       tACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTa                                                               <  1:4065666/150‑1 (MQ=34)

                                                                                          gCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGc                                                            <  2:4298417/150‑1 (MQ=34)

                                                                                          gCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGc                                                            <  2:1302538/150‑1 (MQ=34)

                                                                                            tCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGt                                                          >  1:2816638/1‑150 (MQ=34)

                                                                                                gTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCg                                                      >  1:5012051/1‑150 (MQ=34)

                                                                                                      tGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATc                                                >  2:3816489/1‑150 (MQ=34)

                                                                                                       gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                               >  2:4107793/1‑150 (MQ=34)

                                                                                                         ccAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTg                                             <  2:1871716/150‑1 (MQ=34)

                                                                                                         ccAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTg                                             <  1:1935070/150‑1 (MQ=34)

                                                                                                            ggTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTc                                          >  2:4655137/1‑150 (MQ=34)

                                                                                                              tCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                        >  2:4018906/1‑102 (MQ=32)

                                                                                                                aCGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTcc                                      <  2:2233047/150‑1 (MQ=34)

                                                                                                                aCGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTcc                                      >  1:1551799/1‑150 (MQ=34)

                                                                                                                aCGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTcc                                      <  2:2192481/150‑1 (MQ=34)

                                                                                                                 cGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCt                                     <  1:3240564/150‑1 (MQ=34)

                                                                                                                  gCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTg                                    >  2:2435026/1‑150 (MQ=34)

                                                                                                                        tACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAg                              >  1:1712462/1‑150 (MQ=34)

                                                                                                                        tACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAg                              >  1:3607245/1‑150 (MQ=34)

                                                                                                                          cGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGc                            <  2:580654/150‑1 (MQ=34)

                                                                                                                          cGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGc                            <  2:162135/150‑1 (MQ=34)

                                                                                                                                gaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCt                      >  2:1550158/1‑150 (MQ=37)

                                                                                                                                gaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCt                      >  1:1295845/1‑150 (MQ=37)

                                                                                                                                gaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCt                      >  2:4733490/1‑150 (MQ=37)

                                                                                                                                 aaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGa                                                                                        <  2:4325766/83‑1 (MQ=33)

                                                                                                                                  aCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTgg                    >  1:2297090/1‑150 (MQ=37)

                                                                                                                                  aCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTgg                    <  2:4707364/150‑1 (MQ=37)

                                                                                                                                  aCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTgg                    <  2:3891770/150‑1 (MQ=37)

                                                                                                                                    aTGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTa                  <  2:3813335/150‑1 (MQ=37)

                                                                                                                                    aTGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTa                  <  2:3813486/150‑1 (MQ=37)

                                                                                                                                      gATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAga                <  2:656219/150‑1 (MQ=37)

                                                                                                                                      gATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAga                <  2:2177703/150‑1 (MQ=37)

                                                                                                                                       aTCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgac               <  2:2829902/150‑1 (MQ=37)

                                                                                                                                          aCCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgacgac            >  2:4857366/1‑150 (MQ=37)

                                                                                                                                           ccGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGacgacg           >  1:4526336/1‑150 (MQ=37)

                                                                                                                                                gctgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGaa      <  2:3102439/150‑2 (MQ=2)

                                                                                                                                                 ctgcACAGATGGACGGCGCGATCCTGGTTGTTTCGTCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTAAGGTTTGCGTTACTGCGATCTTCGTGTTCCTGTACAAGTGCGATCAGGTTGACGATGCTGaac     <  1:3455175/150‑3 (MQ=2)

                                                                                                                                                 ctgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGaac     <  1:2530517/150‑3 (MQ=2)

                                                                                                                                                  tgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAact    >  2:499446/1‑150 (MQ=2)

                                                                                                                                                  tgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAact    >  2:1037575/1‑150 (MQ=2)

                                                                                                                                                  tgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAact    >  1:2949122/1‑150 (MQ=2)

                                                                                                                                                    cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc  >  2:2259582/1‑150 (MQ=37)

                                                                                                                                                    cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc  >  2:133389/1‑150 (MQ=37)

GACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGAACATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGC  >  NC_002947/549470‑549767

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