BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,224,926	+C	coding (1052/1434 nt)	ygjI →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,224,924	1	.	C	98.3%	314.5 / ‑1.4	119	coding (1050/1434 nt)	ygjI	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (59/58); total (60/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 128 total aligned reads displayed.

TCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑CC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCA  >  NC_000913/3224848‑3224994

tCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑t                                                                           <  2:1902440/76‑1 (MQ=255)

tCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑t                                                                           >  1:2838294/1‑76 (MQ=255)

 cGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑c                                                                         >  2:2334086/1‑76 (MQ=255)

 cGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑c                                                                         <  1:1045579/76‑1 (MQ=255)

  gAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cc                                                                        >  2:3969976/1‑76 (MQ=255)

    aCTGAACGAAAACGGCGTTCCGGCTC‑GCGCAGCGTGGATCCAGTTCCTGCTCGTCATCCCGCTGATGATTA‑T‑ccct                                                                      <  1:2999544/76‑3 (MQ=255)

     cTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccga                                                                     >  1:4022075/1‑73 (MQ=255)

     cTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccga                                                                     <  1:860581/76‑4 (MQ=255)

         aCGAAATCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCt                                                                 <  1:2491675/76‑1 (MQ=255)

         aCGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCt                                                                 >  1:442439/1‑76 (MQ=255)

         aCGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCt                                                                 <  2:3338059/76‑1 (MQ=255)

          cGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTc                                                                <  1:2932863/76‑1 (MQ=255)

           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCg                                                               >  1:3253521/1‑76 (MQ=255)

           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCg                                                               <  2:4134851/76‑1 (MQ=255)

           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCg                                                               <  1:1551690/76‑1 (MQ=255)

           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCg                                                               >  1:1586658/1‑76 (MQ=255)

             acaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGt                                                             <  2:628480/74‑1 (MQ=255)

             aaaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGt                                                             <  1:695032/76‑1 (MQ=255)

              aaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGtt                                                            >  1:4299930/1‑76 (MQ=255)

               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTc                                                           >  2:2401652/1‑76 (MQ=255)

               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTc                                                           <  1:4006063/76‑1 (MQ=255)

               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTc                                                           <  2:2945283/76‑1 (MQ=255)

                cGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTcc                                                          <  1:734222/76‑1 (MQ=255)

                cGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTcc                                                          >  2:685715/1‑76 (MQ=255)

                cGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTcc                                                          <  2:4306975/76‑1 (MQ=255)

                 ggCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCTTCATCCCGCTGATGATTA‑TCCC‑GCTGCTCCGTTCC‑a                                                        >  2:812070/1‑76 (MQ=255)

                 ggCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑a                                                        >  2:3363996/1‑76 (MQ=255)

                  gCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑ac                                                       >  2:564551/1‑75 (MQ=255)

                   cGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AAt                                                      <  1:847773/76‑1 (MQ=255)

                   cGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AAt                                                      >  2:3003552/1‑76 (MQ=255)

                    gTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATa                                                     >  2:2083603/1‑76 (MQ=255)

                      ttgggcgc‑gcgcaGCGTGGGTCCGGTTGCTGATCGTCAGCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATAcc                                                   <  2:291693/73‑1 (MQ=255)

                      tCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATAcc                                                   >  1:2834616/1‑76 (MQ=255)

                      gccGGCG‑GGCGCAGGGTGGATCCAGGTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATAcc                                                   <  2:284286/75‑1 (MQ=255)

                       ccGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCg                                                  >  2:4164641/1‑76 (MQ=255)

                       ccGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCg                                                  >  1:2162142/1‑76 (MQ=255)

                       ccGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGACATGCCGCTGATGATTG‑T‑CCAGATGCTCGGTTCC‑AATaaag                                                  >  1:4228896/1‑73 (MQ=255)

                        cggcgc‑gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGt                                                 >  2:4733813/1‑76 (MQ=255)

                        cggcgc‑gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGt                                                 <  1:1654389/76‑1 (MQ=255)

                        cggcgc‑gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGt                                                 <  2:4038248/76‑1 (MQ=255)

                        cggcgc‑gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGt                                                 <  1:1761597/76‑1 (MQ=255)

                         ggcgc‑gcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTg                                                >  2:2768685/1‑76 (MQ=255)

                         ggcgc‑gcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTg                                                >  1:2769007/1‑76 (MQ=255)

                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCa                                              <  1:125597/73‑1 (MQ=255)

                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCa                                              >  1:919499/4‑76 (MQ=255)

                             gcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAg                                             >  2:4601853/3‑76 (MQ=255)

                                cgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGAt                                          >  2:1573755/1‑76 (MQ=255)

                                cgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGAt                                          <  2:3320180/76‑1 (MQ=255)

                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGGTTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATc                                         <  2:658856/76‑1 (MQ=255)

                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGAt                                          <  1:1824191/75‑1 (MQ=255)

                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGAt                                          >  2:1824190/1‑75 (MQ=255)

                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATc                                         >  1:597469/1‑76 (MQ=255)

                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATc                                         >  2:839208/1‑76 (MQ=255)

                                  gaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCt                                        <  2:1391594/75‑1 (MQ=255)

                                  cAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑CATACCGTGCAGGATCt                                        >  2:1058928/1‑76 (MQ=255)

                                   cgCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtg                                       <  2:188675/75‑1 (MQ=255)

                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtg                                       >  1:3167983/1‑76 (MQ=255)

                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtg                                       <  2:3417193/76‑1 (MQ=255)

                                   aGCGTGGATCCAGTTCCTGATCGTCACCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtg                                       <  1:3293556/76‑1 (MQ=255)

                                    gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtga                                      >  2:831682/1‑76 (MQ=255)

                                      gTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtgatg                                    >  2:1604284/1‑76 (MQ=255)

                                      gTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtgatg                                    >  1:4487131/1‑76 (MQ=255)

                                       tGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCtgatga                                   >  2:1305602/1‑76 (MQ=255)

                                       tGGATCCAGTTCCTGATCGTAATCCCGCTGATGATTA‑TCCC‑GTTGCGCGGTTCC‑AATACCGTGCAAGATCtgatga                                   >  1:717257/1‑76 (MQ=255)

                                        ggATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATAcc                                                   <  2:3737942/59‑1 (MQ=255)

                                        ggATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATAcc                                                   >  1:3737942/1‑59 (MQ=255)

                                        ggATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGaa                                  <  2:3925465/76‑1 (MQ=255)

                                        ggATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGaa                                  <  2:357744/76‑1 (MQ=255)

                                         gATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAAt                                 <  1:1059435/76‑1 (MQ=255)

                                            ccAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACt                              <  1:685715/76‑1 (MQ=255)

                                            ccAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACt                              <  1:4646678/76‑1 (MQ=255)

                                             cAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACta                             <  1:1107072/76‑1 (MQ=255)

                                              cgTTCCTGATCGTCATCCTGCTGATGCTTAGC‑CC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATCCtgt                            <  2:2824425/75‑3 (MQ=255)

                                              aGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtat                            >  2:2907075/1‑76 (MQ=255)

                                              aGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtat                            <  2:4299930/76‑1 (MQ=255)

                                               gTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtatt                           <  1:3663162/76‑1 (MQ=255)

                                               gTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtatt                           >  1:1374651/1‑76 (MQ=255)

                                               gTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtatt                           >  2:2965763/1‑76 (MQ=255)

                                                ttCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtatta                          <  2:1564875/76‑1 (MQ=255)

                                                 tCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACtattat                         <  1:2293548/76‑1 (MQ=255)

                                                  ccTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTattatt                        >  1:331038/1‑76 (MQ=255)

                                                   cTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTAttatta                       <  2:3033447/76‑1 (MQ=255)

                                                   cTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTAttatta                       >  1:3436391/1‑76 (MQ=255)

                                                    tGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTaa                      >  1:326844/1‑76 (MQ=255)

                                                     gATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAat                     <  2:3680359/76‑1 (MQ=255)

                                                         gTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGa                 <  2:769882/76‑1 (MQ=255)

                                                         gTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGa                 >  2:3727214/1‑76 (MQ=255)

                                                           cATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGAcc               <  1:3003552/76‑1 (MQ=255)

                                                           cATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGAcc               >  1:4621290/1‑76 (MQ=255)

                                                           cATCCCGCTGAGGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGAcc               <  2:2672736/76‑1 (MQ=255)

                                                            aTCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccg              >  2:2752088/1‑76 (MQ=255)

                                                            aTCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccg              <  2:1677193/76‑1 (MQ=255)

                                                                 gCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCa         >  1:2965980/1‑76 (MQ=255)

                                                                 gCTGATGATTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCa         >  1:1565184/1‑76 (MQ=255)

                                                                    gatgatTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCg      >  1:3348572/1‑76 (MQ=255)

                                                                     atgatTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGCATACTATTATTAATATGACCGCCGCcgcgc     >  1:3588809/1‑75 (MQ=255)

                                                                     atgatTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGt     >  1:3553182/1‑76 (MQ=255)

                                                                       gatTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTcc   >  2:1311623/1‑76 (MQ=255)

                                                                       gatTA‑TCCC‑GATGCT‑GGTTCCAAATCCGGGGCAGGATCTGATAAATATTATTATTAATATGACCGCCGCAGCGTcc   <  1:2421947/76‑1 (MQ=255)

                                                                        atTA‑TCCC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCa  >  1:4742845/1‑76 (MQ=255)

TCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑CC‑GATGCTCGGTTCC‑AATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCA  >  NC_000913/3224848‑3224994

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