BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,880,756	G→T	5.3%	R156R (CGA→AGA)	recF ←	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,880,756	0	G	T	5.3%	337.3 / 4.5	107	R156R (CGA→AGA)	recF	gap repair protein
Reads supporting (aligned to +/- strand): ref base G (54/47); new base T (2/4); total (56/51)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 115 total aligned reads displayed.

ATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAGTTCCG  >  NC_000913/3880607‑3880905

aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCg                                                                                                                                                       >  1:604146/1‑150 (MQ=255)

aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCg                                                                                                                                                       <  1:127955/150‑1 (MQ=255)

  cAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCt                                                                                                                                                     <  1:958542/150‑1 (MQ=255)

  cAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCt                                                                                                                                                     <  1:1703491/150‑1 (MQ=255)

   aGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCATGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCtt                                                                                                                                                    <  2:1342213/150‑1 (MQ=255)

    gCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTg                                                                                                                                                   <  2:734194/150‑1 (MQ=255)

          tCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTg                                                                                                                                             <  2:1717993/150‑1 (MQ=255)

          tCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTg                                                                                                                                             <  2:610306/150‑1 (MQ=255)

            ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCTCAAGCGATTGCt                                                                                                                                           <  1:662900/150‑1 (MQ=255)

            ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCTCAAGCGATTGCt                                                                                                                                           <  1:611217/150‑1 (MQ=255)

            ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCt                                                                                                                                           <  1:655874/150‑1 (MQ=255)

              ccGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCTCAAGCGATTGCTca                                                                                                                                         <  1:1242573/150‑2 (MQ=255)

              ccGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCTCAAGCGATTGCTca                                                                                                                                         <  1:856273/150‑2 (MQ=255)

              ccGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCGTGAGATTGCTcc                                                                                                                                         >  1:923484/1‑150 (MQ=255)

               cgcgATACCGGCGCTATACTCCGCGCACCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCa                                                                                                                                        <  2:631949/150‑1 (MQ=255)

                gcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAg                                                                                                                                       <  2:429112/150‑1 (MQ=255)

                gcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCACATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAg                                                                                                                                       <  2:1292311/150‑1 (MQ=255)

                      ccGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTg                                                                                                                                 <  2:1554699/150‑1 (MQ=255)

                      ccGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGAGCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTg                                                                                                                                 <  2:306803/150‑1 (MQ=255)

                        ggCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCTGCGGGATCAGCGCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGaa                                                                                                                               <  1:1074904/150‑1 (MQ=255)

                        ggCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCATGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGaa                                                                                                                               <  2:342268/150‑1 (MQ=255)

                              ataCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATcc                                                                                                                         >  2:275661/1‑150 (MQ=255)

                               taCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCg                                                                                                                        <  2:587952/150‑1 (MQ=255)

                                aCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCTGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGATAATTCCgg                                                                                                                       >  1:1101360/1‑150 (MQ=255)

                                 cTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCggg                                                                                                                      >  2:791064/1‑150 (MQ=255)

                                 cTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCGCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCggg                                                                                                                      <  1:78394/150‑1 (MQ=255)

                                         gcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTtgt                                                                                                              <  1:1186891/150‑1 (MQ=255)

                                          ccAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGACACCTGGCGCAGCGCCGCATAGCGCTGCTTGCGCAATCGCTTGAGATAGCTCCAGGCGGTGAAAAATCCGGGTTCGTtgtg                                                                                                             >  2:194027/1‑150 (MQ=255)

                                            aGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGaa                                                                                                           >  2:1519504/1‑150 (MQ=255)

                                              gTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAg                                                                                                         >  1:369726/1‑150 (MQ=255)

                                                 cTGATTTGCTCCGCCAGCGGGATCCGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCAt                                                                                                      >  1:1539519/1‑150 (MQ=255)

                                                   gATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATcc                                                                                                    >  2:1225861/1‑150 (MQ=255)

                                                   gATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATcc                                                                                                    >  2:1356705/1‑150 (MQ=255)

                                                   gATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATcc                                                                                                    <  2:914429/150‑1 (MQ=255)

                                                    aTTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATccc                                                                                                   <  1:1582913/150‑1 (MQ=255)

                                                     tttGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATcccc                                                                                                  <  2:975723/150‑1 (MQ=255)

                                                     tttGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATcccc                                                                                                  <  1:592399/150‑1 (MQ=255)

                                                     tttGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATcccc                                                                                                  <  1:1711174/150‑1 (MQ=255)

                                                         cTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTc                                                                                              >  1:3774/1‑150 (MQ=255)

                                                          tCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCg                                                                                             <  1:1220045/150‑1 (MQ=255)

                                                          tCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCg                                                                                             <  1:1290302/150‑1 (MQ=255)

                                                           ccGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCgc                                                                                            >  1:1287370/1‑149 (MQ=255)

                                                           ccGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGa                                                                                            <  1:466056/150‑1 (MQ=255)

                                                             gCCAGCTGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACTTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAgg                                                                                          <  2:496969/150‑1 (MQ=255)

                                                             gCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAgg                                                                                          >  2:828952/1‑150 (MQ=255)

                                                                     gATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGctct                                                                                  >  2:494178/1‑150 (MQ=255)

                                                                       tCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCtcttc                                                                                <  2:40693/150‑1 (MQ=255)

                                                                       tCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCtcttc                                                                                <  1:197669/150‑1 (MQ=255)

                                                                       tCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCtcttc                                                                                >  2:794882/1‑150 (MQ=255)

                                                                        cAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTACAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCtcttct                                                                               >  1:1579852/1‑150 (MQ=255)

                                                                          gCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGt                                                                             >  2:807418/1‑150 (MQ=255)

                                                                            tctTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTAt                                                                           >  1:1703369/1‑150 (MQ=255)

                                                                             ctTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTAtt                                                                          >  1:167647/1‑150 (MQ=255)

                                                                                tATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAACATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTgg                                                                       >  1:151445/1‑150 (MQ=255)

                                                                                aaTCCCACGGGCGTAGCTGTTCGTAACGTATCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGATTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTgg                                                                       >  1:467527/2‑150 (MQ=255)

                                                                                  tCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTgggg                                                                     >  1:12251/1‑150 (MQ=255)

                                                                                   cccgctggCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGgc                                                                    <  2:200396/144‑1 (MQ=255)

                                                                                    ccACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGgcc                                                                   >  1:464521/1‑150 (MQ=255)

                                                                                    ccACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGgcc                                                                   >  1:382058/1‑150 (MQ=255)

                                                                                       cGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGgccgcc                                                                >  1:434107/1‑150 (MQ=255)

                                                                                        gggCGTAGCTGTTCGTAACGTGTCACGTGGCGCAGCGCCGCATTGCGCTGCTTGATCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGATGAATGCTCTTCTGTATTTGGGGccgccg                                                               <  1:1089457/150‑1 (MQ=255)

                                                                                           cGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTg                                                            <  2:166478/150‑1 (MQ=255)

                                                                                                  gTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGagtaaag                                                     <  1:910696/150‑1 (MQ=255)

                                                                                                  gTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGagtaaag                                                     >  1:590008/1‑150 (MQ=255)

                                                                                                      gTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCTTCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAgtaaagtaaa                                                 >  1:490082/1‑150 (MQ=255)

                                                                                                      gTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGAATTTGGGGCCGCCGTGGAgtaaagtaaa                                                 >  1:496714/1‑150 (MQ=255)

                                                                                                       gcaCGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTGCGTTGTGAAAGCATCCCCAGGCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAc                                                <  1:831253/148‑1 (MQ=255)

                                                                                                         aCGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAccc                                              >  2:956584/1‑150 (MQ=255)

                                                                                                         aCGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCGTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAccc                                              <  2:293744/150‑1 (MQ=255)

                                                                                                          cGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCt                                             >  2:874439/1‑150 (MQ=255)

                                                                                                            tgtCACCGGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTc                                           <  2:1533001/150‑1 (MQ=255)

                                                                                                             gtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGCATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAccggtct                                          >  1:86191/1‑145 (MQ=255)

                                                                                                             gtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAccagtct                                          >  1:369714/1‑145 (MQ=255)

                                                                                                             gtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCt                                          >  2:1553570/1‑150 (MQ=255)

                                                                                                               cACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTgg                                        <  1:791064/150‑1 (MQ=255)

                                                                                                                aCCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGc                                       <  2:115829/150‑1 (MQ=255)

                                                                                                                aCCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGc                                       <  1:1711546/150‑1 (MQ=255)

                                                                                                                 ccTGGCGCTGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAACCCCTTCTGGCg                                      >  1:1620038/1‑150 (MQ=255)

                                                                                                                     gcgcAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTAt                                  >  2:572170/1‑150 (MQ=255)

                                                                                                                          gcgcCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACt                             >  2:1531623/1‑150 (MQ=255)

                                                                                                                             ccgcATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCa                          >  1:105576/1‑150 (MQ=255)

                                                                                                                             caatgcggcgcTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCa                          >  2:1287374/8‑150 (MQ=255)

                                                                                                                                   tGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCa                    <  2:1324249/150‑1 (MQ=255)

                                                                                                                                   tGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGGGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGTCGTTATCAACTGCATTGGCa                    <  1:1538159/150‑1 (MQ=255)

                                                                                                                                    gcgcTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCAt                   >  2:123102/1‑150 (MQ=255)

                                                                                                                                         gctTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGt              >  1:1457307/1‑150 (MQ=255)

                                                                                                                                           tggAGCAATCTCAAGCGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTgc            >  2:1484956/3‑150 (MQ=255)

                                                                                                                                            tGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAGCTGCATTGGCATCAGGTgcg           >  1:1065487/1‑150 (MQ=255)

                                                                                                                                            tGAGCAATCGCTTGAGATTGCGCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGTAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTgcg           <  1:1055622/150‑1 (MQ=255)

                                                                                                                                               gCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGGTCTTCTGTATTTGGGGCCGCCGTTCAGTAAAGTAAACCCATCTGCCTTTATCATCTACATTGGCATCACGTGCGCCa        >  2:849024/1‑150 (MQ=255)

                                                                                                                                               gCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCa        >  1:397517/1‑150 (MQ=255)

                                                                                                                                                cAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAACCCCATCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAg       >  1:855467/1‑150 (MQ=255)

                                                                                                                                                 aaTCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAGt      >  2:1029142/1‑150 (MQ=255)

                                                                                                                                                   tCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGGGCGTTGTGAAAGCATCGCCAGTCGAGGAATGCTCTTCTGTCTTTGGGGCCGCCGTGGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAGTTc    <  1:194027/150‑1 (MQ=255)

                                                                                                                                                    cGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAGTTcc   >  2:706095/1‑150 (MQ=255)

                                                                                                                                                     gCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGATTGCTCTTCTGTATTTGGTGCCGCGGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGACAGTTCcg  >  1:917268/1‑150 (MQ=255)

                                                                                                                                                     gCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAATCCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCCGTTCcg  >  1:1149797/1‑150 (MQ=255)

                                                                                                                                                     gCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTATGTATTTGGGGCCGCCGTTTAGTAAAGTAAACCCTACTGTCGTTATCAACTGCATTGTCATCAGGTGCGCAAGTTCcg  >  1:1634275/1‑150 (MQ=255)

                                                                                                                                                     gCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTACGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCAGTATTAGGGGCCGCCGATGAGAAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCAAGTTCcg  >  2:1089339/1‑150 (MQ=255)

                                                                                                                                                     gCTTGAGATTGCTCCAGGCGGTGAAAAATCAGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAGTTCcg  >  1:691096/1‑150 (MQ=255)

ATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTTATCAACTGCATTGGCATCAGGTGCGCCAGTTCCG  >  NC_000913/3880607‑3880905

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`