BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,051,163	Δ1 bp	coding (94/231 nt)	sfa →	Sfa

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,051,162	0	A	.	99.5%	1028.7 / ‑3.8	219	coding (93/231 nt)	sfa	Sfa
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (99/119); minor base G (0/1); total (99/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

ATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAC  >  NC_000913/1051070‑1051250

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          <  2:1084651/93‑1 (MQ=255)

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          >  2:5534992/1‑93 (MQ=255)

  gcgcCGCTGGATTTCACAAAATAATATCAGGCTACCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTagc                                                                                        >  2:4934897/1‑91 (MQ=255)

   cgccgcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTagcg                                                                                       >  1:1604476/1‑90 (MQ=255)

   cgccgcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTagcg                                                                                       >  2:3312072/1‑90 (MQ=255)

    gccgcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTagcga                                                                                      <  2:6317987/93‑5 (MQ=255)

     ccgcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGAt                                                                                    <  2:900548/93‑1 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  1:5897243/1‑93 (MQ=255)

         tGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAAcc                                                                                >  1:1022342/1‑93 (MQ=255)

          ggATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCt                                                                               <  2:3345568/93‑1 (MQ=255)

          ggATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCt                                                                               >  1:5466502/1‑93 (MQ=255)

           gATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTg                                                                              <  2:1962438/93‑1 (MQ=255)

           gATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTg                                                                              <  2:706715/93‑1 (MQ=255)

            aTTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGt                                                                             <  1:4214422/93‑1 (MQ=255)

             tttCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGtt                                                                            <  2:4147144/93‑1 (MQ=255)

                cacaAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGat                                                                         <  2:2244423/93‑1 (MQ=255)

                cacaAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGat                                                                         <  1:2945492/93‑1 (MQ=255)

                cacaAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGat                                                                         <  1:2090096/93‑1 (MQ=255)

                 acaAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGatt                                                                        <  2:2822286/93‑1 (MQ=255)

                   aaaaTAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGattat                                                                      >  2:1624763/1‑93 (MQ=255)

                       taataTCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGaa                                                                  <  2:6004766/93‑1 (MQ=255)

                          tatCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCt                                                               >  1:1157531/1‑93 (MQ=255)

                          tatCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCt                                                               <  2:1833504/93‑1 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           >  1:2358019/1‑93 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           >  1:4748953/1‑93 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           <  1:2375373/93‑1 (MQ=255)

                               ggCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCgg                                                          <  1:4878849/93‑1 (MQ=255)

                                 cTCCCTCGTGGAGCCTTTTTTATATCTGCGTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCgggg                                                        <  1:6350464/93‑1 (MQ=255)

                                 cTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCgggg                                                        >  2:3206582/1‑93 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGga                                                       >  2:2730081/1‑93 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGga                                                       >  1:6208272/1‑93 (MQ=255)

                                   cccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgag                                                      >  2:2874087/1‑93 (MQ=255)

                                   cccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgag                                                      >  1:4578019/1‑93 (MQ=255)

                                    ccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgaga                                                     >  1:3060608/1‑93 (MQ=255)

                                    ccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgaga                                                     <  1:979882/93‑1 (MQ=255)

                                     cTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGatat                                                    >  1:3100152/1‑93 (MQ=255)

                                     cTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGAt                                                    <  1:3672972/93‑1 (MQ=255)

                                     cTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGAt                                                    >  1:5026912/1‑93 (MQ=255)

                                      tCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATg                                                   <  2:30631/93‑1 (MQ=255)

                                       cGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATgg                                                  >  2:1808488/1‑93 (MQ=255)

                                         tGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCt                                                <  1:6471813/93‑1 (MQ=255)

                                          ggAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCtt                                               >  2:3743997/1‑93 (MQ=255)

                                           gAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTta                                              <  1:5529516/93‑1 (MQ=255)

                                           gAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTta                                              >  1:3948269/1‑93 (MQ=255)

                                               cTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAAc                                          >  2:1295287/1‑93 (MQ=255)

                                                 tttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAt                                        <  1:3433200/93‑1 (MQ=255)

                                                  ttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       <  1:818177/93‑1 (MQ=255)

                                                    ttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTc                                     <  1:4434118/93‑1 (MQ=255)

                                                     tataTCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCt                                    <  1:5657647/93‑1 (MQ=255)

                                                     tataTCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCt                                    <  1:348820/93‑1 (MQ=255)

                                                        atCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAc                                 >  1:1558305/1‑93 (MQ=255)

                                                         tCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAcc                                >  2:4809006/1‑93 (MQ=255)

                                                         tCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAcc                                <  2:2020006/93‑1 (MQ=255)

                                                         tCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAcc                                >  2:4035128/1‑93 (MQ=255)

                                                          cTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCt                               <  2:6436317/93‑1 (MQ=255)

                                                           tGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTg                              <  1:214537/93‑1 (MQ=255)

                                                           tGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTg                              <  1:1241391/93‑1 (MQ=255)

                                                           tGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTg                              >  1:3264332/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             >  2:6305164/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             >  1:2590843/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             <  1:330097/93‑1 (MQ=255)

                                                             ccTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAc                            >  2:4607066/1‑93 (MQ=255)

                                                             ccTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAc                            <  2:2999989/93‑1 (MQ=255)

                                                             ccTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAc                            <  2:523566/93‑1 (MQ=255)

                                                              cTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAcc                           >  1:496442/1‑93 (MQ=255)

                                                              cTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAcc                           >  1:4541981/1‑93 (MQ=255)

                                                              cTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAcc                           <  2:2392702/93‑1 (MQ=255)

                                                               ttATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCa                          >  2:2315587/1‑93 (MQ=255)

                                                                tATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAg                         <  2:4905600/93‑1 (MQ=255)

                                                                 aTTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAgg                        <  2:4331539/93‑1 (MQ=255)

                                                                 aTTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAgg                        >  1:4409767/1‑93 (MQ=255)

                                                                  tttttCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAggg                       <  2:5413611/93‑1 (MQ=255)

                                                                  tttttCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAggg                       <  2:584627/93‑1 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    <  2:5783939/93‑1 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   >  1:3741781/1‑93 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   >  2:745662/1‑93 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   >  2:2047809/1‑93 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   <  1:771505/93‑1 (MQ=255)

                                                                       cttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCg                  >  2:6125373/1‑93 (MQ=255)

                                                                         tcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCggg                >  1:5226432/1‑93 (MQ=255)

                                                                         tcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCggg                >  2:2272373/1‑93 (MQ=255)

                                                                           aaCGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGaa              <  2:3569192/93‑1 (MQ=255)

                                                                           aaCGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGaa              >  2:4204210/1‑93 (MQ=255)

                                                                            aCGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAAc             >  2:4676083/1‑93 (MQ=255)

                                                                             cGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAAcc            >  1:1334039/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  1:3170274/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  1:1866069/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  1:6013711/1‑93 (MQ=255)

                                                                                  tatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAacac       >  1:2770638/1‑93 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     <  1:3502158/93‑1 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     >  1:3847832/1‑93 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     <  1:2496041/93‑1 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    <  2:2155075/93‑1 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    <  1:4315435/93‑1 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    <  2:5921365/93‑1 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    >  2:3233987/1‑93 (MQ=255)

                                                                                     gtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCtt    <  1:4931507/93‑1 (MQ=255)

                                                                                       ataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAc  <  1:2695665/93‑1 (MQ=255)

                                                                                       ataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAc  <  2:5051560/93‑1 (MQ=255)

                                                                                       ataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAc  <  2:4426987/93‑1 (MQ=255)

ATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAC  >  NC_000913/1051070‑1051250

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