BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,601,264	C→A	5.1%	Q171H (CAG→CAT)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,601,264	0	C	A	5.1%	558.2 / 12.6	168	Q171H (CAG→CAT)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand): ref base C (92/66); new base A (4/5); total (97/71)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGC  >  NC_002947/4601116‑4601413

cTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCt                                                                                                                                                      >  2:279480/1‑150 (MQ=255)

 ttGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTGCTGGCTGAGAAAGGCAGCGGGCGGGTCGATGAGCACCGCAGGATTGCTGCGCCGCCCCCCCTCGGGGTCGGCGAGGGCGTg                                                                                                                                                     >  2:3831842/1‑150 (MQ=255)

      ttCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTg                                                                                                                                                <  2:4553318/150‑1 (MQ=255)

        ccAGACGGTACCGAAGACTCCACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAACGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTggc                                                                                                                                              <  2:4755007/150‑1 (MQ=255)

        ccAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTggc                                                                                                                                              <  1:1207099/150‑1 (MQ=255)

         cAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTggcc                                                                                                                                             <  2:4070907/150‑1 (MQ=255)

          aGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                                            >  1:3972607/1‑150 (MQ=255)

                   cGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGc                                                                                                                                   >  1:3311166/1‑150 (MQ=255)

                    gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGcc                                                                                                                                  >  1:3321497/1‑150 (MQ=255)

                         cTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                                             >  1:10956/1‑150 (MQ=255)

                          tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                            >  2:4784402/1‑150 (MQ=255)

                          tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                                            >  1:544552/1‑150 (MQ=255)

                           cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                                           <  1:4574174/150‑1 (MQ=255)

                            aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCtt                                                                                                                          >  2:3316069/1‑150 (MQ=255)

                            aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCtt                                                                                                                          >  1:4445378/1‑150 (MQ=255)

                                      gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCgg                                                                                                                <  2:1601944/150‑1 (MQ=255)

                                      gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCgg                                                                                                                <  1:3751283/150‑1 (MQ=255)

                                      gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCgg                                                                                                                <  1:2248002/150‑1 (MQ=255)

                                           aCTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCaa                                                                                                           <  1:2557274/150‑1 (MQ=255)

                                            cTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                          <  2:1968714/150‑1 (MQ=255)

                                            cTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                          <  1:4361234/150‑1 (MQ=255)

                                            cTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                                          <  1:1545828/150‑1 (MQ=255)

                                                ccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                                      >  2:1315373/1‑150 (MQ=255)

                                                  aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                                    <  1:3280057/150‑1 (MQ=255)

                                                  aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTAGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                                    >  2:2537791/1‑150 (MQ=255)

                                                  aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTAGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                                    >  2:2537869/1‑150 (MQ=255)

                                                   gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAAc                                                                                                   >  1:2610714/1‑150 (MQ=255)

                                                   gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAAc                                                                                                   >  1:2323284/1‑150 (MQ=255)

                                                      aaTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACttt                                                                                                <  1:4777913/150‑1 (MQ=255)

                                                      aaTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACttt                                                                                                <  1:1950308/150‑1 (MQ=255)

                                                       aTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTg                                                                                               <  1:1746189/150‑1 (MQ=255)

                                                       aTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTg                                                                                               <  1:259341/150‑1 (MQ=255)

                                                        tCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGa                                                                                              <  2:4224047/150‑1 (MQ=255)

                                                         cTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                                             >  2:1863570/1‑150 (MQ=255)

                                                          tGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGacc                                                                                            <  2:1661452/150‑3 (MQ=255)

                                                           gCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCg                                                                                           >  1:2298225/1‑150 (MQ=255)

                                                             ggggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGt                                                                                         >  2:1960373/1‑150 (MQ=255)

                                                               ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGc                                                                                       <  2:2803766/150‑1 (MQ=255)

                                                               ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGc                                                                                       <  1:1640583/150‑1 (MQ=255)

                                                                     tCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCaa                                                                                 <  2:2208648/150‑1 (MQ=255)

                                                                         gCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATaaa                                                                             >  2:3473490/1‑150 (MQ=255)

                                                                           tGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgc                                                                           >  1:4227200/1‑150 (MQ=255)

                                                                           tGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgc                                                                           >  1:2890557/1‑150 (MQ=255)

                                                                           tGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgc                                                                           >  2:3555162/1‑150 (MQ=255)

                                                                            gATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgct                                                                          <  2:497118/150‑1 (MQ=255)

                                                                                 aGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGTAACTTTGATCGGTGCAGTCAATAAAGCTGgtct                                                                     >  2:749811/1‑150 (MQ=255)

                                                                                   gCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTg                                                                   >  2:3217308/1‑150 (MQ=255)

                                                                                      ccGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCg                                                                <  2:4188207/150‑1 (MQ=255)

                                                                                           cgggTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctg                                                           <  2:2127586/150‑1 (MQ=255)

                                                                                               tCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                       >  1:4571715/1‑150 (MQ=255)

                                                                                               tCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                                       <  2:1944380/150‑1 (MQ=255)

                                                                                                 gatgatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgc                                                     >  2:2111056/1‑150 (MQ=255)

                                                                                                 gatgatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgc                                                     >  2:4877002/1‑150 (MQ=255)

                                                                                                 gatgatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgc                                                     >  2:4492448/1‑150 (MQ=255)

                                                                                                  atgatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                                    >  1:1238560/1‑150 (MQ=255)

                                                                                                     atCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGccc                                                 >  1:821930/1‑150 (MQ=255)

                                                                                                      tCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCa                                                <  1:4198853/150‑1 (MQ=255)

                                                                                                      tCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCa                                                <  2:2930607/150‑1 (MQ=255)

                                                                                                        aCCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGc                                              <  1:2823744/150‑1 (MQ=255)

                                                                                                         ccGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCa                                             >  2:1631729/1‑150 (MQ=255)

                                                                                                         ccGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCa                                             >  2:4212245/1‑150 (MQ=255)

                                                                                                          cGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGGGCTGGCGGGCCCCGTGCCTGCGCCTTGTCGATGGGGGTCAACTATAAACTTTGATGGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGGGCCCAGCAc                                            <  2:928037/150‑1 (MQ=255)

                                                                                                          cGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCag                                            >  1:950667/1‑149 (MQ=255)

                                                                                                          cGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCAc                                            >  1:1838888/1‑150 (MQ=255)

                                                                                                            catcatTGCTGCTCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGa                                          <  2:4523503/150‑1 (MQ=255)

                                                                                                            catcatTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGa                                          <  2:2627311/150‑1 (MQ=255)

                                                                                                            catcatTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGa                                          <  1:1315373/150‑1 (MQ=255)

                                                                                                            catcatTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGa                                          <  2:573092/150‑1 (MQ=255)

                                                                                                             atcatTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCGACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                                         <  1:1344926/150‑1 (MQ=255)

                                                                                                                 ttGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGACACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTAc                                     >  2:3766170/1‑150 (MQ=255)

                                                                                                                   gCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCa                                   <  2:1257254/150‑1 (MQ=255)

                                                                                                                   gCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCa                                   >  1:877933/1‑150 (MQ=255)

                                                                                                                     tGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACg                                 >  2:4380810/1‑150 (MQ=255)

                                                                                                                      gATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAAt                                                                                >  2:4951767/1‑102 (MQ=255)

                                                                                                                      gATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAAt                                                                                <  1:4951767/102‑1 (MQ=255)

                                                                                                                       aTCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGc                               <  2:773700/150‑1 (MQ=255)

                                                                                                                        tcgcctaccCCGTGGGGT‑GGTGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCaa                             >  1:2414836/7‑150 (MQ=255)

                                                                                                                          caccacCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAAt                            >  1:1391981/1‑150 (MQ=255)

                                                                                                                          caccacCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAAt                            >  1:694461/1‑150 (MQ=255)

                                                                                                                           accaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCtt                                                                                                                          <  2:4167564/55‑1 (MQ=255)

                                                                                                                           accaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCtt                                                                                                                          >  1:1586962/1‑55 (MQ=255)

                                                                                                                           accaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATc                           >  2:1534862/1‑150 (MQ=255)

                                                                                                                           accaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATc                           >  1:340292/1‑150 (MQ=255)

                                                                                                                            ccaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                          <  2:3830915/150‑1 (MQ=255)

                                                                                                                            ccaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTATGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg                          >  1:3354488/1‑150 (MQ=255)

                                                                                                                             caccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGTGGTCCCACGGcgcgcgg                         >  2:3922635/1‑144 (MQ=255)

                                                                                                                             caccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCGCGGCACTCgg                         >  2:4890154/1‑150 (MQ=255)

                                                                                                                             caccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCgg                         <  2:3972607/150‑1 (MQ=255)

                                                                                                                               ccccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa                       <  1:341695/150‑1 (MQ=255)

                                                                                                                                cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTAc                      >  2:3233971/1‑150 (MQ=255)

                                                                                                                                cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTGTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTCCCACGGCAATCGGTAc                      >  2:176302/1‑150 (MQ=255)

                                                                                                                                 ccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACa                     >  2:4838552/1‑150 (MQ=255)

                                                                                                                                     ggggccggCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcgc                 >  1:2706927/7‑150 (MQ=255)

                                                                                                                                           ggCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACg           >  1:3079481/1‑150 (MQ=255)

                                                                                                                                            gccaGCGCAGGCGCTCGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGc          >  2:5111225/4‑150 (MQ=255)

                                                                                                                                            gCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGc          >  2:3527517/1‑150 (MQ=255)

                                                                                                                                              gAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGccc        >  2:1448632/1‑150 (MQ=255)

                                                                                                                                                  gcCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCg    <  1:4117054/150‑1 (MQ=255)

                                                                                                                                                    cTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCaa                                                                                 <  2:2353664/71‑1 (MQ=255)

                                                                                                                                                    cTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGc  <  1:4063337/150‑1 (MQ=255)

CTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCGCGCGACGCCCACCGGC  >  NC_002947/4601116‑4601413

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