BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,484,849	G→A	R355H (CGC→CAC)	yjgP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,484,849	0	G	A	99.0%	1459.0 / NA	422	R355H (CGC→CAC)	yjgP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); major base A (204/210); minor base C (1/3); total (205/217)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 489 total aligned reads displayed.

TGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTCACCAC  >  NC_000913/4484759‑4484938

 ggACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCt                                                                                        <  2:5830967/93‑1 (MQ=255)

   aCCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgc                                                                                      >  1:3866425/1‑93 (MQ=255)

   aCCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgc                                                                                      <  1:8615584/93‑1 (MQ=255)

         cGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGccg                                                                                >  2:3278916/1‑92 (MQ=255)

         cGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCa                                                                                >  1:2826230/1‑93 (MQ=255)

          gCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAg                                                                               <  1:5345387/93‑1 (MQ=255)

            tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAggt                                                                             >  1:630629/1‑91 (MQ=255)

            tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                             >  1:7432854/1‑93 (MQ=255)

             gtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGttt                                                                            >  1:5520441/1‑93 (MQ=255)

             gtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGttt                                                                            <  1:9716909/93‑1 (MQ=255)

              tggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtttt                                                                           <  1:2613310/93‑1 (MQ=255)

                  tgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgc                                                                       >  2:562574/1‑93 (MQ=255)

                  tgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgc                                                                       <  2:1528303/93‑1 (MQ=255)

                   gtggaCGGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgcg                                                                      >  2:7202676/1‑93 (MQ=255)

                     ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                                    <  1:7165155/93‑1 (MQ=255)

                      gaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTGCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaa                                                                   <  2:5468376/93‑1 (MQ=255)

                       aCCGTTAACCTGATTTATCTGGCTTTAGCTATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaca                                                                  >  2:5126242/1‑91 (MQ=255)

                         cGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCTCCTGCGCGCCAGTTTTTCGCGTAAAgg                                                                <  2:3602295/93‑1 (MQ=255)

                         cGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAgg                                                                <  1:5612390/93‑1 (MQ=255)

                         cGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGGTTTTCGCGTGAAgg                                                                >  2:6116529/1‑93 (MQ=255)

                          gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                               >  1:251673/1‑93 (MQ=255)

                          gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                               >  1:8404023/1‑93 (MQ=255)

                           ttAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAg                                                              <  2:3915579/93‑1 (MQ=255)

                            tAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGc                                                             <  1:8143028/93‑1 (MQ=255)

                             taCCTGATTTATCTTGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCg                                                            <  1:7683546/92‑1 (MQ=255)

                             aaCCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCg                                                            <  1:9919010/93‑1 (MQ=255)

                              aCCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCAGCCACCTGCGCGCCAGCTTTTAGCGTAAAGGAGCgg                                                           >  1:3447798/1‑93 (MQ=255)

                                cTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtg                                                         <  1:709325/93‑1 (MQ=255)

                                cTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtg                                                         >  1:1029873/1‑93 (MQ=255)

                                 tGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtgt                                                        >  2:7322976/1‑93 (MQ=255)

                                 tGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtgt                                                        <  2:2027771/93‑1 (MQ=255)

                                   aTTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGa                                                      >  1:5798135/1‑93 (MQ=255)

                                   aTTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGa                                                      <  2:8497025/93‑1 (MQ=255)

                                   aTTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGa                                                      <  2:3801377/93‑1 (MQ=255)

                                    tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                                     >  1:4890497/1‑93 (MQ=255)

                                    tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                                     >  1:2266263/1‑93 (MQ=255)

                                        gcTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGTCAGCTTTTCGCGTAAATGATCGGTGTGATGCaa                                                 >  1:8938335/2‑93 (MQ=255)

                                          tGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAAcc                                               <  2:9503676/93‑1 (MQ=255)

                                           ggCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCt                                              >  1:8493499/1‑93 (MQ=255)

                                           ggCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCt                                              >  2:1405376/1‑93 (MQ=255)

                                            gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                             >  2:9713188/1‑93 (MQ=255)

                                             cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                            <  1:6510336/93‑1 (MQ=255)

                                             cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                            >  1:1399724/1‑93 (MQ=255)

                                              tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                           <  2:532126/93‑1 (MQ=255)

                                              tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                           <  1:4279630/93‑1 (MQ=255)

                                              tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCACGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                           <  1:3188539/93‑1 (MQ=255)

                                               ttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTg                                          <  2:8690031/93‑1 (MQ=255)

                                               ttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGGACAGGAGCGGTGTGATGCAACCTTTTg                                          >  1:588119/1‑93 (MQ=255)

                                                 aGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGc                                        <  1:683706/93‑1 (MQ=255)

                                                  gCGATTGTTCTCAACCTTTGGGACACCGCGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCg                                       >  2:7980069/1‑93 (MQ=255)

                                                   cGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGt                                      >  2:3526201/1‑93 (MQ=255)

                                                   cGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGt                                      <  2:6210622/93‑1 (MQ=255)

                                                   cGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGt                                      <  2:1812584/93‑1 (MQ=255)

                                                   cGATTGTTCTCAACATTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGt                                      <  2:7591200/93‑1 (MQ=255)

                                                     aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTat                                    >  1:3330903/1‑92 (MQ=255)

                                                     aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTat                                    >  1:9290507/1‑92 (MQ=255)

                                                     aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                    >  2:2575567/1‑93 (MQ=255)

                                                      ttGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGTGGTGTGATGCCACCTTTTGGCGTACt                                   >  1:1882770/1‑93 (MQ=255)

                                                         ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                                >  1:5020247/1‑93 (MQ=255)

                                                          tctcAACCTTTGGGAGACCGTGCCGGGCCGGCACCTGCGCGCCAGGTTTTCGCGTAACGGAGCGTTGTGATGCCACAGTGCGGCGTACCTGAc                               >  2:1760864/1‑93 (MQ=255)

                                                          tctcAACCTTTGGGACACCGTGACGCCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCATCCTTTTGGCGTACTTGAc                               <  2:5267864/93‑1 (MQ=255)

                                                           ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                              <  2:445457/93‑1 (MQ=255)

                                                           ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                              >  2:123816/1‑93 (MQ=255)

                                                           ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                              >  2:2452266/1‑93 (MQ=255)

                                                            tcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCg                             >  1:1659183/1‑93 (MQ=255)

                                                              aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                           >  1:2666278/1‑93 (MQ=255)

                                                               aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                          >  2:4163677/1‑93 (MQ=255)

                                                               aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                          <  2:2338892/93‑1 (MQ=255)

                                                               aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                          >  2:6789780/1‑93 (MQ=255)

                                                               aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGGGTGATGCAACCTTTTGGCGTACTTGACCGCta                          >  2:773197/1‑93 (MQ=255)

                                                                ccTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtat                         <  2:2802897/93‑1 (MQ=255)

                                                                  tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                       >  2:1135317/1‑93 (MQ=255)

                                                                  tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                       >  2:6017012/1‑93 (MQ=255)

                                                                  tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                       <  1:779333/93‑1 (MQ=255)

                                                                    tGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCg                     >  2:8844918/1‑93 (MQ=255)

                                                                     gggACACCGTGCCGGCCCGCCCCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAAGCTTTTGGCGTACTTGACCGCTCTACCgg                    >  1:5231012/1‑93 (MQ=255)

                                                                     gggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCgg                    >  1:3988929/1‑93 (MQ=255)

                                                                      ggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt                   <  1:1288876/93‑1 (MQ=255)

                                                                      ggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt                   >  2:9036732/1‑93 (MQ=255)

                                                                      ggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt                   <  2:3041639/93‑1 (MQ=255)

                                                                        gcacCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaa                 <  1:751040/92‑1 (MQ=255)

                                                                          acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa               <  2:5037092/93‑1 (MQ=255)

                                                                           ccGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAAc              >  2:6435068/1‑93 (MQ=255)

                                                                           ccGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAAc              <  1:6192431/93‑1 (MQ=255)

                                                                              tGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAt           >  2:5565099/1‑93 (MQ=255)

                                                                              tGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAt           >  2:4589652/1‑93 (MQ=255)

                                                                               gCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtt          >  1:916493/1‑93 (MQ=255)

                                                                               gCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtt          <  2:4012415/93‑1 (MQ=255)

                                                                                ccGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAttt         <  1:4428948/93‑1 (MQ=255)

                                                                                  ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc       <  2:815979/90‑1 (MQ=255)

                                                                                  ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc       >  2:4921618/4‑93 (MQ=255)

                                                                                  ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc       <  2:7856195/90‑1 (MQ=255)

                                                                                   gcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTca      >  1:10029181/3‑93 (MQ=255)

                                                                                   gcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTca      >  1:6685493/3‑93 (MQ=255)

                                                                                   gcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTca      <  1:1127125/91‑1 (MQ=255)

                                                                                   gcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTca      <  2:9839885/91‑1 (MQ=255)

                                                                                    cccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcac     >  2:6384668/2‑93 (MQ=255)

                                                                                    cccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcac     >  1:3014364/2‑93 (MQ=255)

                                                                                     ccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcacc    <  1:4634855/93‑1 (MQ=255)

                                                                                       gccaccTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcaccac  <  2:9268642/93‑1 (MQ=255)

TGGACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTCACCAC  >  NC_000913/4484759‑4484938

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