BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	855,672	C→G	5.5%	L174V (CTA→GTA)	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	855,672	0	C	G	5.5%	482.0 / 13.6	145	L174V (CTA→GTA)	murI	glutamate racemase
Reads supporting (aligned to +/- strand): ref base C (77/60); new base G (4/4); total (81/64)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 146 total aligned reads displayed.

GAGTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCC  >  NC_002947/855524‑855821

gagTGCCAAGTTCGCAGCCTTGCTTGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCAGCt                                                                                                                                                      <  1:3266086/150‑1 (MQ=255)

     ccAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTgg                                                                                                                                                 >  2:1163799/1‑150 (MQ=255)

     ccAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTgg                                                                                                                                                 <  1:2217742/150‑1 (MQ=255)

       aaGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCt                                                                                                                                               >  2:4820441/1‑150 (MQ=255)

               aGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTg                                                                                                                                       >  2:3908916/1‑150 (MQ=255)

               aGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTg                                                                                                                                       >  1:2432497/1‑150 (MQ=255)

                gCCGTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCacataac                                                                                                                                      <  2:5293489/150‑7 (MQ=255)

                  cTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGa                                                                                                                                    <  1:4182998/150‑1 (MQ=255)

                  cTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGa                                                                                                                                    <  2:297942/150‑1 (MQ=255)

                  cTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATGACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGa                                                                                                                                    <  2:2595377/150‑1 (MQ=255)

                    tGCTCGATCGCTTCGCAAATGACGTACAGGTTATTCCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATa                                                                                                                                  <  2:1435908/150‑1 (MQ=255)

                       tCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGc                                                                                                                               <  2:1885191/150‑1 (MQ=255)

                       tCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGc                                                                                                                               >  2:1018308/1‑150 (MQ=255)

                         gATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTg                                                                                                                             >  2:4926315/1‑150 (MQ=255)

                          aTCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGa                                                                                                                            <  1:5089558/150‑1 (MQ=255)

                          aTCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGa                                                                                                                            >  2:3979455/1‑150 (MQ=255)

                            cGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATc                                                                                                                          <  2:5051007/150‑1 (MQ=255)

                              cTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCt                                                                                                                        <  2:4049874/150‑1 (MQ=255)

                               ttCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAg                                                                                                                                           >  2:2752502/1‑130 (MQ=255)

                               ttCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTc                                                                                                                       <  2:3797832/150‑1 (MQ=255)

                                 cGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCgg                                                                                                                     <  2:2351803/150‑1 (MQ=255)

                                 cGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCgg                                                                                                                     <  1:3486546/150‑1 (MQ=255)

                                      aTGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCa                                                                                                                >  1:3647047/1‑150 (MQ=255)

                                       tGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCAc                                                                                                               >  1:2142862/1‑150 (MQ=255)

                                           gTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCAt                                                                                                           >  1:4299701/1‑150 (MQ=255)

                                           gTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCAt                                                                                                           >  1:4266667/1‑150 (MQ=255)

                                             aCAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTa                                                                                                         <  2:2363434/150‑1 (MQ=255)

                                              cAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTAc                                                                                                        <  1:525938/150‑1 (MQ=255)

                                               aGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTAATGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTAcc                                                                                                       >  2:5219943/1‑150 (MQ=255)

                                                ggTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTAccc                                                                                                      <  1:3725062/150‑1 (MQ=255)

                                                 gTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCt                                                                                                     >  2:1323075/1‑150 (MQ=255)

                                                 gTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCt                                                                                                     >  1:442836/1‑150 (MQ=255)

                                                  ttattaCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCtt                                                                                                    <  2:2077456/150‑1 (MQ=255)

                                                    attaCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTc                                                                                                  >  1:1734427/1‑150 (MQ=255)

                                                       aCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCtt                                                                                               >  1:277864/1‑150 (MQ=255)

                                                        cccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATcc                                                                                                                         <  1:1303344/123‑1 (MQ=255)

                                                        cccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATcc                                                                                                                         >  2:1303344/1‑123 (MQ=255)

                                                        cccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGACTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTc                                                                                              >  1:4602936/1‑150 (MQ=255)

                                                         ccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCg                                                                                             >  1:10346/1‑150 (MQ=255)

                                                         ccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCg                                                                                             >  1:1423969/1‑150 (MQ=255)

                                                          cAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGc                                                                                            >  2:1217938/1‑150 (MQ=255)

                                                                 gCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGc                                                                                     <  1:5027393/150‑1 (MQ=255)

                                                                  cccGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCt                                                                                    <  2:1423969/150‑1 (MQ=255)

                                                                  cccGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCt                                                                                    >  2:576534/1‑150 (MQ=255)

                                                                     ggggCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGc                                                                                 <  1:4134801/150‑1 (MQ=255)

                                                                     ggggCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGc                                                                                 <  1:1984916/150‑1 (MQ=255)

                                                                      gggCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGcc                                                                                >  2:4738015/1‑150 (MQ=255)

                                                                          tGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCa                                                                            <  1:4311394/150‑1 (MQ=255)

                                                                            gTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATg                                                                          >  2:3680471/1‑150 (MQ=255)

                                                                               gAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCac                                                                                                                                           <  1:3863667/82‑2 (MQ=255)

                                                                                 gCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTccc                                                                     >  1:3283659/1‑150 (MQ=255)

                                                                                 gCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTccc                                                                     >  2:4735846/1‑150 (MQ=255)

                                                                                    gATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGa                                                                  <  1:4735846/150‑1 (MQ=255)

                                                                                    gATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGa                                                                  >  2:744433/1‑150 (MQ=255)

                                                                                     aTCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGAt                                                                 >  2:3371423/1‑150 (MQ=255)

                                                                                     aTCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCACATTGCTGGCCGGCATGGTCCCGGAt                                                                 <  1:5435116/150‑1 (MQ=255)

                                                                                        gAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGAc                                                              >  2:458932/1‑150 (MQ=255)

                                                                                          aGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGt                                                            >  2:4891314/1‑150 (MQ=255)

                                                                                           gCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTg                                                           <  1:1266917/150‑1 (MQ=255)

                                                                                                 gACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGt                                                                        >  2:3272141/1‑131 (MQ=255)

                                                                                                  aCCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCcatca                                                    <  1:3212021/150‑1 (MQ=255)

                                                                                                    cTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCatcatc                                                  >  1:4513494/1‑150 (MQ=255)

                                                                                                    cTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCatcatc                                                  <  2:1231092/150‑1 (MQ=255)

                                                                                                    cTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCatcatc                                                  >  1:945397/1‑150 (MQ=255)

                                                                                                     gggTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCg                                                 <  2:691164/149‑1 (MQ=255)

                                                                                                       gTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCAGTGCTGGCCGGCATGGTCCCGGAGGACGGGGCCATCATCggc                                               >  2:2887423/1‑148 (MQ=255)

                                                                                                           gCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCg                                           >  2:379517/1‑150 (MQ=255)

                                                                                                           gCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCg                                           >  1:2276286/1‑150 (MQ=255)

                                                                                                             ccGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGgc                                         >  2:2426510/1‑150 (MQ=255)

                                                                                                             ccGGCGGTGCACCAGCTGCTGCCAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCAGTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGgc                                         >  1:5271023/1‑150 (MQ=255)

                                                                                                              cggcggTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGgcg                                        >  1:3526668/1‑150 (MQ=255)

                                                                                                              cggcggTGCGCCAGCTGCTGCAAGGTTATGAGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACCCCGgcg                                        <  1:5448726/150‑1 (MQ=255)

                                                                                                                     gcgcCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTGTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGt                                 <  2:3095272/150‑1 (MQ=255)

                                                                                                                     gcgcCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGACCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGt                                 <  2:2278620/150‑1 (MQ=255)

                                                                                                                      cgcCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGaca                                              >  1:2564770/1‑136 (MQ=255)

                                                                                                                      cgcCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGaca                                              <  2:2564770/136‑1 (MQ=255)

                                                                                                                        ccAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGgc                              >  2:1852043/1‑150 (MQ=255)

                                                                                                                         cAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGgcg                             >  2:5131437/1‑150 (MQ=255)

                                                                                                                          aGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCTATCATCGACACCGGCGCGGCCGTGgcgc                            >  1:1813911/1‑150 (MQ=255)

                                                                                                                              gctgcACGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGACATCATCGACACCGGCGCGGCCGTGGCGCGTCa                        >  2:5371271/1‑150 (MQ=255)

                                                                                                                                tgcAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGt                      >  2:1974682/1‑150 (MQ=255)

                                                                                                                                tgcAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATAAGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACAACGGCGAGGCCGTGGCGCGTCAGt                      >  1:2256263/1‑150 (MQ=255)

                                                                                                                                 gcAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGtt                     <  1:3845644/150‑1 (MQ=255)

                                                                                                                                  cAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTg                    >  1:3423766/1‑150 (MQ=255)

                                                                                                                                  cAAGGTTATGTGCAGACGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGTCTGCACCCATTACCCTTTCCTTCGCCCAATGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGTCGCGTCAGTTg                    >  1:3706875/1‑150 (MQ=255)

                                                                                                                                   aaGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGc                   >  1:1646728/1‑150 (MQ=255)

                                                                                                                                       ttATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACg               <  2:4266667/150‑1 (MQ=255)

                                                                                                                                        tATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACgg              <  1:3113513/150‑1 (MQ=255)

                                                                                                                                        tATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACgg              >  2:3117385/1‑150 (MQ=255)

                                                                                                                                          tgtgCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCt            >  2:662235/1‑150 (MQ=255)

                                                                                                                                           gtgCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATg                                                                          <  1:2330344/87‑1 (MQ=255)

                                                                                                                                           ctgCAGCCAGTAGCGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTg           >  1:1389663/2‑150 (MQ=255)

                                                                                                                                            tgCAGCCAGTAGCGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTg                                <  2:751434/128‑1 (MQ=255)

                                                                                                                                            tgCAGCCAGTAGCGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTg                                >  1:751434/1‑128 (MQ=255)

                                                                                                                                             gCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGtt         <  1:4730887/150‑1 (MQ=255)

                                                                                                                                              cAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTg        <  1:4926315/150‑1 (MQ=255)

                                                                                                                                               aGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTgg       >  2:4685756/1‑150 (MQ=255)

                                                                                                                                               aGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTgg       >  1:750574/1‑150 (MQ=255)

                                                                                                                                                   gCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGc   <  2:1646559/150‑1 (MQ=255)

                                                                                                                                                   gCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGc   <  2:1646728/150‑1 (MQ=255)

GAGTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCC  >  NC_002947/855524‑855821

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