BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,939	G→C	8.5%	G310A (GGC→GCC) ‡	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,939	0	G	C	8.5%	476.8 / 26.0	149	G310A (GGC→GCC) ‡	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base G (71/64); new base C (11/2); total (82/67)
Fisher's exact test for biased strand distribution p-value = 3.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

GTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCAC  >  NC_002947/5754790‑5755085

gTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACgg                                                                                                                                                    <  1:2676822/150‑1 (MQ=255)

gTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACgg                                                                                                                                                    <  1:1376877/150‑1 (MQ=255)

 tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                                                                                                   >  2:1806377/1‑150 (MQ=255)

    ggCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTc                                                                                                                                                <  1:3371683/150‑1 (MQ=255)

     gCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCa                                                                                                                                               <  2:721754/150‑1 (MQ=255)

      cGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                                                                                                                              <  2:3254073/150‑1 (MQ=255)

      cGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                                                                                                                              <  2:4659799/150‑1 (MQ=255)

             ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                       >  1:197947/1‑150 (MQ=255)

             ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                       <  1:3039082/150‑1 (MQ=255)

             ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                       >  1:2474429/1‑150 (MQ=255)

             ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                       >  1:1032759/1‑150 (MQ=255)

                  gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCtt                                                                                                                                                 >  1:2739447/1‑135 (MQ=255)

                  gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGAc                                                                                                                                  <  1:3108748/150‑1 (MQ=255)

                  gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGAc                                                                                                                                  >  1:5124812/1‑150 (MQ=255)

                        gCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGTTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                                                                                                                            <  2:1786810/150‑1 (MQ=255)

                          aCGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGt                                                                                                                          <  2:3586627/150‑1 (MQ=255)

                            gTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATTGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTAc                                                                                                                        >  2:2668411/1‑150 (MQ=255)

                              gCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACaa                                                                                                                      <  2:1943905/150‑1 (MQ=255)

                                  ttctCGATCTGGAAAGAAGGCAAGCCAACGTTAGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCAATAAGCGt                                                                                                                  >  2:4021242/4‑150 (MQ=255)

                                  aTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGt                                                                                                                  <  1:4732319/150‑1 (MQ=255)

                                   tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                                                                                                 <  2:3907007/150‑1 (MQ=255)

                                     tCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTg                                                                                                               >  1:2025948/1‑150 (MQ=255)

                                      cGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTgg                                                                                                              <  1:1849958/150‑1 (MQ=255)

                                      cGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTgg                                                                                                              <  1:1113884/150‑1 (MQ=255)

                                       gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                                                                                                             >  1:4865/1‑150 (MQ=255)

                                       gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                                                                                                             <  2:1353469/150‑1 (MQ=255)

                                       gATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGt                                                                                                             >  2:2501139/1‑150 (MQ=255)

                                        aTCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTATGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTc                                                                                                            >  1:4803874/1‑150 (MQ=255)

                                          cTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTccc                                                                                                          <  1:2630438/150‑1 (MQ=255)

                                           tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                         >  1:3932204/1‑150 (MQ=255)

                                           tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                         <  1:5240303/150‑1 (MQ=255)

                                           tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                         >  2:230523/1‑150 (MQ=255)

                                           tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                         >  2:230665/1‑150 (MQ=255)

                                            ggAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                                    <  2:4350140/138‑1 (MQ=255)

                                            ggAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                                    >  1:4350140/1‑138 (MQ=255)

                                                aGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCtt                                                                                                    <  1:760135/150‑1 (MQ=255)

                                                    ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGATCGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                                <  2:3429364/150‑1 (MQ=255)

                                                    ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                                <  2:4225981/150‑1 (MQ=255)

                                                    ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                                >  1:1548501/1‑150 (MQ=255)

                                                     gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                                                                                               >  2:3308831/1‑150 (MQ=255)

                                                     gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                                                                                               <  2:3766908/150‑1 (MQ=255)

                                                     gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                                                                                               <  1:4322215/150‑1 (MQ=255)

                                                         gAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcccc                                                                                           <  1:3644694/150‑1 (MQ=255)

                                                               cTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAAt                                                                                     >  1:3946779/1‑150 (MQ=255)

                                                                ttCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTaggc                                                                                    >  2:457558/1‑147 (MQ=255)

                                                                   tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                                                                                                                            <  2:3356326/107‑1 (MQ=255)

                                                                   tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                 <  2:925893/150‑1 (MQ=255)

                                                                      ggTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGcc                                                                              <  1:4676378/150‑1 (MQ=255)

                                                                          aaGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACt                                                                          <  1:230523/150‑1 (MQ=255)

                                                                          aaGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACt                                                                          <  1:230665/150‑1 (MQ=255)

                                                                            ggCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCg                                                                        >  1:3658575/1‑150 (MQ=255)

                                                                             gCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCgg                                                                       >  1:1574506/1‑150 (MQ=255)

                                                                               tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGAcc                                                                                                                                 >  1:5334272/1‑90 (MQ=255)

                                                                               tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCt                                                                     >  1:4205371/1‑150 (MQ=255)

                                                                               tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCt                                                                     >  2:504822/1‑150 (MQ=255)

                                                                                 tGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCg                                                                   >  1:1375584/1‑150 (MQ=255)

                                                                                  gCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCgt                                                                  >  1:3088564/1‑149 (MQ=255)

                                                                                  gCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCgt                                                                  >  1:3088834/1‑149 (MQ=255)

                                                                                     ggCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAAc                                                               >  1:1861110/1‑150 (MQ=255)

                                                                                      gCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAAcc                                                              >  2:2154271/1‑150 (MQ=255)

                                                                                        cTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGt                                                            >  1:1600652/1‑150 (MQ=255)

                                                                                         tGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGtt                                                           >  1:908393/1‑150 (MQ=255)

                                                                                           tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGc                                                                               >  1:2150189/1‑128 (MQ=255)

                                                                                           tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGc                                                                               <  2:2150189/128‑1 (MQ=255)

                                                                                               acacCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTc                                                                         <  1:2227319/130‑1 (MQ=255)

                                                                                                  ccGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGa                                                  >  1:3814812/1‑150 (MQ=255)

                                                                                                    tgccTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgc                                                >  2:4940629/3‑146 (MQ=255)

                                                                                                       cTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGt                                             >  1:4033511/1‑150 (MQ=255)

                                                                                                        tGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTa                                            <  2:523877/150‑1 (MQ=255)

                                                                                                          ttcTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTAtt                                          <  1:1873268/148‑1 (MQ=255)

                                                                                                            cTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTcc                                        >  2:4751188/1‑150 (MQ=255)

                                                                                                               cATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTc                                                                                                                                     <  2:1554871/54‑1 (MQ=255)

                                                                                                               cATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTa                                     >  1:1851759/1‑150 (MQ=255)

                                                                                                                aTCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTAc                                    >  2:1563053/1‑150 (MQ=255)

                                                                                                                 tCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACg                                   >  2:4037974/1‑150 (MQ=255)

                                                                                                                  cggcggTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGt                                  >  1:626464/1‑150 (MQ=255)

                                                                                                                       gTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCTACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCa                             <  1:727067/150‑1 (MQ=255)

                                                                                                                       gTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCa                             <  2:2844833/150‑1 (MQ=255)

                                                                                                                          tcatcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGcc                          <  2:3689691/150‑1 (MQ=255)

                                                                                                                           catcaAGCACGGTAAGGCCCTGAACGTCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGccc                         <  2:2851692/150‑1 (MQ=255)

                                                                                                                           catcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGccc                         <  2:3932204/150‑1 (MQ=255)

                                                                                                                            atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                 >  2:3722822/1‑93 (MQ=255)

                                                                                                                            atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                 <  1:3722822/93‑1 (MQ=255)

                                                                                                                            atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTcc                                                         >  2:4771212/1‑117 (MQ=255)

                                                                                                                             tcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGa                       >  2:4081964/1‑150 (MQ=255)

                                                                                                                             tcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGa                       >  1:901279/1‑150 (MQ=255)

                                                                                                                               aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                     <  2:4203126/150‑1 (MQ=255)

                                                                                                                               aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                     <  2:1138155/150‑1 (MQ=255)

                                                                                                                                aGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                     >  2:4168868/1‑149 (MQ=255)

                                                                                                                                aGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                     <  1:4169106/149‑1 (MQ=255)

                                                                                                                                     ggTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGa               >  2:1891539/1‑150 (MQ=255)

                                                                                                                                     ggTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGa               >  1:1400450/1‑150 (MQ=255)

                                                                                                                                      gTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGAc              <  1:2441365/150‑1 (MQ=255)

                                                                                                                                           gCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATc         >  1:3526737/1‑150 (MQ=255)

                                                                                                                                             ccTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGa       >  2:4312410/1‑150 (MQ=255)

                                                                                                                                              cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa      <  1:2740215/150‑1 (MQ=255)

                                                                                                                                              cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa      <  2:5325937/150‑1 (MQ=255)

                                                                                                                                              cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa      >  2:550145/1‑150 (MQ=255)

                                                                                                                                                gAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAACGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGc    >  1:396115/1‑150 (MQ=255)

                                                                                                                                                 aaCGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCa   <  1:5260865/150‑1 (MQ=255)

GTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCAC  >  NC_002947/5754790‑5755085

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