BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	855,671	G→A	5.6%	P173P (CCG→CCA)	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	855,671	0	G	A	5.6%	447.2 / 10.6	143	P173P (CCG→CCA)	murI	glutamate racemase
Reads supporting (aligned to +/- strand): ref base G (77/58); new base A (4/4); total (81/62)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 145 total aligned reads displayed.

GAGTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGC  >  NC_002947/855524‑855820

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)

           tCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCac                                                                                                                                          <  1:3466299/150‑2 (MQ=255)

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

               aGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTg                                                                                                                                      <  2:2875054/150‑1 (MQ=255)

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

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

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

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

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

                          aTCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTCCTGGCTGCCGGTTGCGATACGCTGa                                                                                                                           >  2:3979463/1‑150 (MQ=255)

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

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

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

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

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

                                   cAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCt                                                                                                                  <  1:681104/150‑1 (MQ=255)

                                   cAAATGACGTACAGGTTATGACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCt                                                                                                                  <  1:4239792/150‑1 (MQ=255)

                                      aTGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACACTGATCCTCGGCTGCa                                                                                                               >  1:2237922/1‑150 (MQ=255)

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

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

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

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

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

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

                                                  ttattaCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCtt                                                                                                   <  1:839950/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)

                                                         ccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCAccc                                                                                                            >  2:2388896/1‑134 (MQ=255)

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

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

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

                                                           aGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGcc                                                                                          >  1:4650772/1‑150 (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:3680528/1‑150 (MQ=255)

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

                                                                            gTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATg                                                                         <  2:1297034/150‑1 (MQ=255)

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

                                                                               gAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCac                                                                                                                                          >  2:3863667/1‑81 (MQ=255)

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

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

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

                                                                                     aTCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGAt                                                                <  1:5435080/150‑1 (MQ=255)

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

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

                                                                                       cGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGa                                                              >  1:5247454/1‑150 (MQ=255)

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

                                                                                        gAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGAc                                                             >  2:1484415/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                                                 >  1:945397/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                  ggTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGc                                   >  2:1742490/1‑150 (MQ=255)

                                                                                                                    tGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCg                                 >  2:4945790/1‑150 (MQ=255)

                                                                                                                     gcgcCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGt                                <  1:4489101/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                           gctgctGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGcgcg                          >  2:1044481/1‑150 (MQ=255)

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

                                                                                                                                tgcAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGt                     >  2:1974682/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)

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

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

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

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

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

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

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

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

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

GAGTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGC  >  NC_002947/855524‑855820

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