BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	934,839	T→C	V798A (GTT→GCT)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	934,839	0	T	C	99.5%	602.3 / NA	197	V798A (GTT→GCT)	ftsK	cell division protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (99/97); total (99/98)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

CACAACAGCCGGTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAT  >  NC_000913/934749‑934927

cacaACAGCCGGTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTg                                                                                        >  2:3151009/1‑93 (MQ=255)

 ccaacaGCCGGTT‑CCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgc                                                                                       <  2:5965596/92‑1 (MQ=255)

   agcaGCCGGTT‑CCGCCGCAGCCACAATAGCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgcgc                                                                                     <  2:3319207/91‑1 (MQ=18)

      gcccGTTT‑GCGCCCCAGCCACATGATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGcgccgc                                                                                  <  1:6133662/91‑1 (MQ=18)

       ggcGGTTGGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGcgccgc                                                                                  <  2:4371739/91‑1 (MQ=18)

          ggTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagca                                                                              <  1:6506752/93‑1 (MQ=255)

          ggTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagca                                                                              >  1:2909734/1‑93 (MQ=255)

          ggTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagca                                                                              >  1:2788747/1‑93 (MQ=255)

           gTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcag                                                                             >  2:1962274/1‑93 (MQ=255)

                cgccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCAGCTGCGCCGCAGCAGCAGt                                                                         <  1:6538933/93‑1 (MQ=255)

                ccccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGt                                                                         <  1:1517158/91‑1 (MQ=255)

                ccccgcACCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGt                                                                         <  1:1168350/91‑1 (MQ=255)

                ccccgcACCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGt                                                                         <  2:4169808/91‑1 (MQ=18)

                 cccccAGCCACAATACCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTa                                                                        <  1:5677251/89‑1 (MQ=18)

                  cgccAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTAt                                                                       <  2:5505321/90‑1 (MQ=255)

                  ccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTAt                                                                       <  2:5389007/93‑1 (MQ=255)

                  ccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTAt                                                                       <  2:346749/93‑1 (MQ=255)

                  ccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTAt                                                                       <  2:1569566/93‑1 (MQ=255)

                      aGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagc                                                                   <  1:5023106/93‑1 (MQ=255)

                       gCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagca                                                                  >  1:767989/1‑93 (MQ=255)

                         cccaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCagcagc                                                                <  1:2944427/91‑1 (MQ=255)

                         cacaATATCAGCAGCCGCATCAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCagcagc                                                                <  2:6026241/93‑1 (MQ=255)

                         cacaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCagcagc                                                                >  1:6406553/1‑93 (MQ=255)

                          acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGcc                                                               >  1:1264040/1‑93 (MQ=255)

                            aaTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCcac                                                             >  2:6415719/1‑93 (MQ=255)

                              tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACCACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                           >  2:4642049/1‑93 (MQ=255)

                              tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                           >  2:4102909/1‑93 (MQ=255)

                               ataagcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaaa                                                          <  2:5889541/90‑2 (MQ=255)

                                tCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaaca                                                         >  2:3391732/1‑93 (MQ=255)

                                  agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGc                                                       >  1:4698552/1‑93 (MQ=255)

                                   gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGcc                                                      >  1:4861470/1‑93 (MQ=255)

                                   gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGcc                                                      >  2:3676936/1‑93 (MQ=255)

                                   ccagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGcc                                                      <  1:2356109/92‑1 (MQ=255)

                                       ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                  >  2:3378704/1‑93 (MQ=255)

                                        cGCAACACCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTg                                                 <  1:2008958/93‑1 (MQ=255)

                                        cGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTg                                                 <  2:2143902/93‑1 (MQ=255)

                                         gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgc                                                >  2:470284/1‑93 (MQ=255)

                                           aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                              <  2:5571934/93‑1 (MQ=255)

                                           aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                              >  2:3094854/1‑93 (MQ=255)

                                            acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGcc                                             >  1:1433032/1‑93 (MQ=255)

                                              aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcac                                           >  2:6611458/1‑93 (MQ=255)

                                              aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcac                                           >  2:3998345/1‑93 (MQ=255)

                                              aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcac                                           >  2:6710890/1‑93 (MQ=255)

                                               acCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                          >  1:3694804/1‑93 (MQ=255)

                                               acCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                          >  2:189981/1‑93 (MQ=255)

                                                ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                         >  2:2375805/1‑93 (MQ=255)

                                                ccGGTTGCCCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                         <  1:4202871/93‑1 (MQ=255)

                                                 cGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAc                                        <  2:1091834/93‑1 (MQ=255)

                                                  ggTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAcc                                       >  2:4873626/1‑93 (MQ=255)

                                                       cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAgg                                  >  1:3217606/1‑93 (MQ=255)

                                                        gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGa                                 <  2:621884/93‑1 (MQ=255)

                                                        cccACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGa                                 <  2:57200/92‑1 (MQ=255)

                                                        cccACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGa                                 <  1:4630742/92‑1 (MQ=255)

                                                         ccACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGAt                                >  1:2323443/1‑93 (MQ=255)

                                                          cacaGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATa                               <  1:242354/93‑1 (MQ=255)

                                                            cagaagCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAcc                             <  1:2581097/93‑1 (MQ=255)

                                                             agcagcATTATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                            <  1:3592307/93‑1 (MQ=255)

                                                             agcagaAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                            <  2:3986561/93‑1 (MQ=255)

                                                             agcagaAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                            <  1:1870092/93‑1 (MQ=255)

                                                             agcaccAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                            <  2:2791362/93‑1 (MQ=255)

                                                              gcagcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCt                           >  2:248965/1‑93 (MQ=255)

                                                              gcagcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCt                           <  2:3636147/93‑1 (MQ=255)

                                                              gcagcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCt                           >  2:75471/1‑93 (MQ=255)

                                                                agcaGTATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                         <  2:2696414/93‑1 (MQ=255)

                                                                agcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                         <  2:6004967/93‑1 (MQ=255)

                                                                 gcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCt                        >  2:894725/1‑93 (MQ=255)

                                                                 gcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCt                        >  1:4030887/1‑93 (MQ=255)

                                                                 gcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCt                        >  1:4949803/1‑93 (MQ=255)

                                                                  caATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCtt                       >  1:2179328/1‑93 (MQ=255)

                                                                   aaTATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                      >  2:456543/1‑93 (MQ=255)

                                                                     tatCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCAt                    >  1:5815877/1‑93 (MQ=255)

                                                                      atCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                   >  1:2259061/1‑93 (MQ=255)

                                                                      atCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                   <  2:1866820/93‑1 (MQ=255)

                                                                      accagcagCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCaac                   <  2:3009438/91‑3 (MQ=255)

                                                                       tCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                  <  2:5781883/93‑1 (MQ=255)

                                                                       tCAGCACCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                  <  1:5359100/93‑1 (MQ=255)

                                                                         agcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGc                <  2:697806/93‑1 (MQ=255)

                                                                          gcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCt               >  1:4769737/1‑93 (MQ=255)

                                                                            agcCGCACAAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt             <  2:1324869/93‑1 (MQ=255)

                                                                            agcCGCAACAGCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt             <  2:6291209/93‑1 (MQ=255)

                                                                            agcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt             <  2:2909734/93‑1 (MQ=255)

                                                                            agcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt             <  2:2863291/93‑1 (MQ=255)

                                                                            acccGCATCAACCGGCTGCGCCGAAGCAGCAGTATCAGCACCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt             <  1:1597998/91‑1 (MQ=255)

                                                                             gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt            >  1:5510160/1‑93 (MQ=255)

                                                                             gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt            >  2:394759/1‑93 (MQ=255)

                                                                             gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt            >  2:2536344/1‑93 (MQ=255)

                                                                             gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt            <  1:882326/93‑1 (MQ=255)

                                                                              ccGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg           >  1:5881742/1‑93 (MQ=255)

                                                                               cGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa          <  2:502967/93‑1 (MQ=255)

                                                                                gCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt         >  1:3491877/1‑93 (MQ=255)

                                                                                gCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt         >  2:4036910/1‑93 (MQ=255)

                                                                                 caacaaCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg        >  1:4595182/1‑93 (MQ=255)

                                                                                  aacaacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGc       >  2:5559991/1‑93 (MQ=255)

                                                                                  aacaacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGc       >  1:377130/1‑93 (MQ=255)

                                                                                   acaacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCg      >  1:1110265/1‑93 (MQ=255)

                                                                                   acaacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCg      >  1:3349980/1‑93 (MQ=255)

                                                                                     aacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa    <  2:2563261/93‑1 (MQ=255)

                                                                                      acCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTaa   >  1:566853/1‑93 (MQ=255)

                                                                                       ccGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt  <  2:5079272/93‑1 (MQ=255)

                                                                                       ccGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt  >  2:1692313/1‑93 (MQ=255)

CACAACAGCCGGTT‑GCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAT  >  NC_000913/934749‑934927

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`