BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,113,340	G→C	6.3%	W23C (TGG→TGC)	PP_5358 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,113,340	0	G	C	6.3%	850.0 / 34.8	254	W23C (TGG→TGC)	PP_5358	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (130/108); new base C (6/10); total (136/118)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 259 total aligned reads displayed.

GCGCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGTCTA  >  NC_002947/6113191‑6113489

gcgcagcgcaGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTgg                                                                                                                                                       <  2:737334/150‑1 (MQ=255)

 cgcagcgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGc                                                                                                                                                      <  1:655566/150‑1 (MQ=255)

      cgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAAc                                                                                                                                                 >  2:4377091/1‑150 (MQ=255)

        cagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCa                                                                                                                                               <  1:3841981/150‑1 (MQ=255)

         agTATGAACTTGGTTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  2:2452176/150‑2 (MQ=255)

         agTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  2:1388854/150‑2 (MQ=255)

         agTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  2:4207345/150‑2 (MQ=255)

         agTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  2:2555049/150‑2 (MQ=255)

          gTATGATCTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATg                                                                                                                                             <  1:3272547/150‑1 (MQ=255)

           tATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATgg                                                                                                                                            <  2:3384529/150‑1 (MQ=255)

              gAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGtt                                                                                                                                         >  1:4613087/1‑150 (MQ=255)

                  ttGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGt                                                                                                                                     <  2:821296/150‑1 (MQ=255)

                  ttGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGt                                                                                                                                     <  1:5613005/150‑1 (MQ=255)

                  ttGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGt                                                                                                                                     <  2:6235294/150‑1 (MQ=255)

                  ttGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGt                                                                                                                                     <  2:785924/150‑1 (MQ=255)

                      ttGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGaaa                                                                                                                                 >  1:2643733/1‑150 (MQ=255)

                            ccTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGccc                                                                                                                           <  2:6407648/150‑1 (MQ=255)

                              tGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTg                                                                                                                         <  2:4441542/150‑1 (MQ=255)

                              tGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTg                                                                                                                         >  2:2116164/1‑150 (MQ=255)

                                   gACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACg                                                                                                                    >  2:3320951/1‑150 (MQ=255)

                                    aCAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACgg                                                                                                                   <  2:3719625/150‑1 (MQ=255)

                                     cAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACggg                                                                                                                  <  2:208505/150‑1 (MQ=255)

                                      aGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGt                                                                                                                 <  2:1017080/150‑1 (MQ=255)

                                        gCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCa                                                                                                               >  2:5520612/1‑150 (MQ=255)

                                        gCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCa                                                                                                               <  1:1682863/150‑1 (MQ=255)

                                          cAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGGAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCaaa                                                                                                             >  1:1258202/1‑150 (MQ=255)

                                              aaCTCAGTTAATCTTGATCTACCTGTTTGGAGCATATTAAAATGTTGACCCTGGATAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGc                                                                                                         <  1:6560821/150‑1 (MQ=255)

                                                   aGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgc                                                                                                    <  2:1973632/150‑1 (MQ=255)

                                                   aGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgc                                                                                                    >  1:2213629/1‑150 (MQ=255)

                                                   aGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgc                                                                                                    <  1:6453575/150‑1 (MQ=255)

                                                     ttAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgcgc                                                                                                  >  2:1114735/1‑150 (MQ=255)

                                                         gCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCaaa                                                                                              >  1:549375/1‑150 (MQ=255)

                                                            ggATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTc                                                                                           <  1:628136/150‑1 (MQ=255)

                                                                 tACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCt                                                                                      >  2:3280444/1‑150 (MQ=255)

                                                                   ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCTCCAAACTCGACCTgg                                                                                    >  2:6547219/1‑150 (MQ=255)

                                                                   ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTgg                                                                                    >  2:6194907/1‑150 (MQ=255)

                                                                   ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTgg                                                                                    >  2:5584492/1‑150 (MQ=255)

                                                                   ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTgg                                                                                    >  2:1249085/1‑150 (MQ=255)

                                                                     tGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAg                                                                                  >  2:4654738/1‑150 (MQ=255)

                                                                         tGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTCGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGc                                                                              <  1:2776775/150‑1 (MQ=255)

                                                                          ggAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCt                                                                             <  1:151848/150‑1 (MQ=255)

                                                                             gCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGa                                                                          >  2:1250887/1‑150 (MQ=255)

                                                                             gCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGa                                                                          <  1:2891900/150‑1 (MQ=255)

                                                                               atatGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACg                                                                        >  1:5199175/1‑150 (MQ=255)

                                                                                tatGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGa                                                                       >  1:6060798/1‑150 (MQ=255)

                                                                                  tGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATg                                                                     >  1:2491174/1‑150 (MQ=255)

                                                                                   gCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGc                                                                    >  1:5802686/1‑150 (MQ=255)

                                                                                     aaaTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCg                                                                  >  2:1609112/1‑150 (MQ=255)

                                                                                     aaaTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCg                                                                  >  2:697726/1‑150 (MQ=255)

                                                                                       aTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTa                                                                <  1:51158/150‑1 (MQ=255)

                                                                                       aTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTa                                                                <  1:4377091/150‑1 (MQ=255)

                                                                                            gACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCt                                                           >  2:1447289/1‑150 (MQ=255)

                                                                                             aCCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTc                                                          >  1:3798929/1‑150 (MQ=255)

                                                                                               ccTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                        >  1:6592113/1‑150 (MQ=255)

                                                                                               ccTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                        <  1:314105/150‑1 (MQ=255)

                                                                                               ccTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                        >  1:3546464/1‑150 (MQ=255)

                                                                                                 tGGATAACCTCTTCGTCCTTATGCTGGTAGCCTCGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCTTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGc                                                      <  1:1965537/150‑1 (MQ=255)

                                                                                                 tGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGc                                                      <  2:5971564/150‑1 (MQ=255)

                                                                                                          tctTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCt                                             >  2:5007198/1‑150 (MQ=255)

                                                                                                          tctTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCt                                             <  2:786195/150‑1 (MQ=255)

                                                                                                           ctTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCtt                                            >  1:256417/1‑150 (MQ=255)

                                                                                                            ttggTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           <  1:1443663/147‑1 (MQ=255)

                                                                                                            ttCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           >  1:3946707/1‑150 (MQ=255)

                                                                                                            ttCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           >  1:503323/1‑150 (MQ=255)

                                                                                                            ttCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           >  1:2297045/1‑150 (MQ=255)

                                                                                                            ttCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           >  2:347124/1‑150 (MQ=255)

                                                                                                                tCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTcc                                       >  2:1491920/1‑150 (MQ=255)

                                                                                                                  cTGATGCTGGTAGCCACGGCAGGCGCCTGGGTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCgg                                     <  1:4762514/150‑1 (MQ=255)

                                                                                                                  atGATGCTGGTAGCCACGGCAGGCACCTGGTGGTGGCACAACCCTGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATGGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCgg                                     >  2:588193/2‑150 (MQ=255)

                                                                                                                   tGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCggg                                    >  2:3737625/1‑150 (MQ=255)

                                                                                                                   tGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCggg                                    <  2:5356227/150‑1 (MQ=255)

                                                                                                                    gATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGc                                               <  1:496687/138‑1 (MQ=255)

                                                                                                                    gATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGa                                   <  1:114307/150‑1 (MQ=255)

                                                                                                                      tGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATg                                 <  2:176265/150‑1 (MQ=255)

                                                                                                                      tGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATg                                 <  2:5737648/150‑1 (MQ=255)

                                                                                                                       gCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGc                                                      >  2:6801712/1‑128 (MQ=255)

                                                                                                                       gCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATgt                                <  1:4654738/150‑2 (MQ=255)

                                                                                                                              gCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGc                         >  1:1007985/1‑150 (MQ=255)

                                                                                                                               ccACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGcc                        <  1:3277598/150‑1 (MQ=255)

                                                                                                                               cacaaccAGGCGCCT‑GCCGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCg                       >  1:2096190/7‑150 (MQ=255)

                                                                                                                                 aCGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                      >  1:3057528/1‑150 (MQ=255)

                                                                                                                                     cAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAgc                  >  2:3192091/1‑150 (MQ=255)

                                                                                                                                        gcgcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGcc               >  1:5029194/1‑150 (MQ=255)

                                                                                                                                        gcgcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGcc               <  2:4842559/150‑1 (MQ=255)

                                                                                                                                        gcgcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGcc               >  2:3771924/1‑150 (MQ=255)

                                                                                                                                        gcgcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGcc               >  2:2148338/1‑150 (MQ=255)

                                                                                                                                         cgcCTGGTTGTGGCACAACCATGGGTTGCGTGACAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt              <  2:5170386/150‑1 (MQ=255)

                                                                                                                                          gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             >  1:1535215/1‑150 (MQ=255)

                                                                                                                                          gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACACGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             >  1:4910302/1‑150 (MQ=255)

                                                                                                                                          ccatGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             >  1:3437039/4‑150 (MQ=255)

                                                                                                                                           ccTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            >  1:5886988/1‑150 (MQ=255)

                                                                                                                                           ccTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            <  2:1892593/150‑1 (MQ=255)

                                                                                                                                           catGGTTGTGCCACAACCATTGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            >  1:6390053/3‑150 (MQ=255)

                                                                                                                                           catGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            >  1:3742039/3‑150 (MQ=255)

                                                                                                                                           catGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            >  1:4937952/3‑150 (MQ=255)

                                                                                                                                           catGGTTGTGCCACAACCAGGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            >  1:572880/3‑150 (MQ=255)

                                                                                                                                             tGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCa          >  1:1949910/1‑150 (MQ=255)

                                                                                                                                                ttGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAcgc       <  1:797938/150‑1 (MQ=255)

                                                                                                                                                  gtgGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGt     >  2:2916669/1‑150 (MQ=255)

                                                                                                                                                     gCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGTCTa  <  2:450572/150‑1 (MQ=255)

GCGCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGTCTA  >  NC_002947/6113191‑6113489

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