BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,481,725	(C)_6→7	coding (887/1023 nt)	yheT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,481,719	1	.	C	95.2%	749.9 / 28.1	249	coding (881/1023 nt)	yheT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (5/7); new base C (114/123); total (119/130)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 276 total aligned reads displayed.

GATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑CCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCT  >  NC_000913/3481627‑3481806

gATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGtc‑t                                                                                           >  2:4336687/1‑93 (MQ=255)

 aTGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGtc‑t‑c                                                                                         >  2:2224312/1‑93 (MQ=255)

   tcTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGGTTATGGATCATCAGGTGATCCCGAAAGCGGAAAGTC‑T‑ccc                                                                                       >  1:637774/2‑93 (MQ=255)

    cTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑cccc                                                                                      <  2:2948390/93‑1 (MQ=255)

      gAACCGGATCGCCAAACCGACGCTGATTATTAACGCCAAAGACGATCCGTTTATGGGTCATCAGGTGATCCCGAAACCGGAAAGCC‑T‑ccccc‑c                                                                                   >  2:3018181/1‑93 (MQ=255)

       aaCCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑ccccc‑cc                                                                                  <  1:1026809/93‑2 (MQ=255)

       aaCCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑ccccc‑cc                                                                                  <  2:2818747/93‑2 (MQ=255)

          cGGATCGCCAACCCGACTCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑ccccc‑ccgca                                                                               <  2:4585654/93‑5 (MQ=255)

            gATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAgg                                                                              <  1:1606316/93‑1 (MQ=255)

            gATCCCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAgg                                                                              <  2:1314175/93‑1 (MQ=255)

            gATCACCANACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAgg                                                                              <  2:2895045/93‑1 (MQ=255)

                 ccAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAg                                                                         >  1:687314/1‑93 (MQ=255)

                 ccAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAg                                                                         <  1:4312140/93‑1 (MQ=255)

                  cAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAga                                                                        >  2:1343573/1‑92 (MQ=255)

                  cAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGt                                                                        <  2:3579626/93‑1 (MQ=255)

                   aaaCCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTa                                                                       >  1:1019155/1‑93 (MQ=255)

                    aaCCGACGCTGATTATTCACCCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGCAACCGGAACGTC‑TCCCCCC‑CGCAGGTGGAGCAt                                                                      >  1:1307744/1‑93 (MQ=255)

                      ccGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCa                                                                    >  2:1921700/1‑93 (MQ=255)

                      ccGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCa                                                                    >  2:3244297/1‑93 (MQ=255)

                      ccGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCa                                                                    <  1:2601826/93‑1 (MQ=255)

                         tcGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAact                                                                 >  1:4425295/2‑93 (MQ=255)

                          cGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCCGGTGGAGTATCAactg                                                                >  2:3414604/1‑93 (MQ=255)

                          cGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAactg                                                                <  2:792394/93‑1 (MQ=255)

                            cTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAactgac                                                              <  1:2700665/93‑1 (MQ=255)

                            cTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAactgac                                                              >  2:4262561/1‑93 (MQ=255)

                            cTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAactgac                                                              >  1:2375245/1‑93 (MQ=255)

                            cTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAactgac                                                              <  1:1498255/93‑1 (MQ=255)

                             tcattattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTCAC‑CCCCC‑CGCAGGTGGAGTATCAactgact                                                             <  1:4139492/91‑1 (MQ=255)

                             tGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCATGTGGAGTATCAactgact                                                             <  2:2348584/93‑1 (MQ=255)

                             tGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAactgact                                                             >  1:3378310/1‑93 (MQ=255)

                               attattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAActgactga                                                           >  2:3311036/1‑93 (MQ=255)

                                ttattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGaa                                                          >  1:4526942/1‑93 (MQ=255)

                                ttattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGaa                                                          >  2:1764370/1‑93 (MQ=255)

                                 tattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAAc                                                         >  1:1861814/1‑93 (MQ=255)

                                 tattCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAAc                                                         >  1:299696/1‑93 (MQ=255)

                                  attCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACa                                                        <  2:2563791/93‑1 (MQ=255)

                                   ttCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACAt                                                       >  2:2542033/1‑93 (MQ=255)

                                    tCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATg                                                      <  2:2061846/93‑1 (MQ=255)

                                    tCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATg                                                      >  1:3672287/1‑93 (MQ=255)

                                    tCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATg                                                      <  2:3837953/93‑1 (MQ=255)

                                     cACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATgg                                                     <  1:4365769/93‑1 (MQ=255)

                                     cACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATgg                                                     <  2:1703692/93‑1 (MQ=255)

                                     cACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATgg                                                     <  2:374975/93‑1 (MQ=255)

                                      aCGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGc                                                    <  2:1109915/93‑1 (MQ=255)

                                         ccAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGt                                                 <  1:56192/93‑1 (MQ=255)

                                          cAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTc                                                <  1:74448/93‑1 (MQ=255)

                                           aaaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGCCTGAACATGGCGGTCa                                               >  2:907311/1‑93 (MQ=255)

                                           aaaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGAATCAACTGACTGAACATGGCGGTCa                                               >  2:4412088/1‑93 (MQ=255)

                                           aaaGACGATCCGTTTATGGATCATCAGGTAATACCGAAACCGGAAAGTC‑TCCCCCC‑CGCTGGAGGAGAATCAACTGACTGAACATGGCGGTCa                                               >  2:2459515/1‑93 (MQ=255)

                                            aaGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCAt                                              >  1:873059/1‑93 (MQ=255)

                                              gtcGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGt                                            <  1:2184593/91‑1 (MQ=255)

                                              gACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGt                                            >  2:706200/1‑93 (MQ=255)

                                              gACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGt                                            <  1:1539861/93‑1 (MQ=255)

                                               aCGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGtt                                           <  2:4181847/93‑1 (MQ=255)

                                                cGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTg                                          <  2:449835/93‑1 (MQ=255)

                                                cGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTg                                          >  1:1398309/1‑93 (MQ=255)

                                                 gATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTgg                                         >  1:3229248/1‑93 (MQ=255)

                                                  aTCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGc                                        >  1:1057967/1‑93 (MQ=255)

                                                  aTCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGc                                        <  2:1777406/93‑1 (MQ=255)

                                                   tCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCt                                       <  2:3079179/93‑1 (MQ=255)

                                                   tCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCt                                       >  1:2060266/1‑93 (MQ=255)

                                                    ccGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCtt                                      <  1:654350/93‑1 (MQ=255)

                                                    ccGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCtt                                      <  1:2979251/93‑1 (MQ=255)

                                                    ccGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCtt                                      >  2:1194515/1‑93 (MQ=255)

                                                      gTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTa                                    >  2:4680632/1‑93 (MQ=255)

                                                         tATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTg                                 >  1:4215702/1‑93 (MQ=255)

                                                         tATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTg                                 <  2:2701453/93‑1 (MQ=255)

                                                         tATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTg                                 <  1:2268561/93‑1 (MQ=255)

                                                           tGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGc                               >  1:3288689/1‑93 (MQ=255)

                                                           tGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGc                               <  1:3594941/93‑1 (MQ=255)

                                                            ggATCATCAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATTAACTGACTGAACATTGCGGTCATGTTGGCTTTATTGGCg                              >  2:4300013/1‑93 (MQ=255)

                                                              atcatcAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGt                            >  1:2853452/1‑93 (MQ=255)

                                                              atcatcAGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGt                            >  2:1555269/1‑93 (MQ=255)

                                                              atcatcAGGTGATCCCGAAACCGGAAAGTC‑T‑CCCCCTCGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGt                            >  2:859294/1‑93 (MQ=255)

                                                                catcaGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTAc                          >  1:4603236/1‑93 (MQ=255)

                                                                 atcaGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGGACa                         >  2:4034682/1‑93 (MQ=255)

                                                                  tcaGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACAt                        >  1:4679008/1‑93 (MQ=255)

                                                                  tcaGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACAt                        <  1:1331780/93‑1 (MQ=255)

                                                                    aGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTa                      <  2:2021913/93‑1 (MQ=255)

                                                                    aGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTa                      <  1:964430/93‑1 (MQ=255)

                                                                    aGGTGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTa                      <  2:2274508/93‑1 (MQ=255)

                                                                       tGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACtt                   >  1:170872/1‑93 (MQ=255)

                                                                       tGATCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACtt                   <  1:1428792/93‑1 (MQ=255)

                                                                         aTCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCa                 >  1:351500/1‑93 (MQ=255)

                                                                          tCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCAt                <  1:2030823/93‑1 (MQ=255)

                                                                          tCCCGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCAt                >  1:3915926/1‑93 (MQ=255)

                                                                           cccGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTTGCGGTACGTTACTTCATc               >  1:1949302/1‑93 (MQ=255)

                                                                           cccGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATc               <  2:1355694/93‑1 (MQ=255)

                                                                           cccGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATc               <  1:1115057/93‑1 (MQ=255)

                                                                            ccGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATcc              <  1:1174825/93‑1 (MQ=255)

                                                                             cGAAACCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCg             >  2:3005611/1‑93 (MQ=255)

                                                                                aaCCGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCaa          >  2:2165825/1‑93 (MQ=255)

                                                                                  ccGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAt        <  2:754844/93‑1 (MQ=255)

                                                                                   cGGAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtg       <  1:583778/93‑1 (MQ=255)

                                                                                     gAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtgtg     <  1:2464841/93‑1 (MQ=255)

                                                                                     gAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtgtg     >  2:1861529/1‑93 (MQ=255)

                                                                                     gAAAGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAAtgtg     <  1:4071135/93‑1 (MQ=255)

                                                                                      aaaGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTgg    <  2:1748440/93‑1 (MQ=255)

                                                                                       aaGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGc   >  2:503573/1‑93 (MQ=255)

                                                                                        aGTC‑TCCCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCt  <  1:1151909/93‑1 (MQ=255)

GATGCTGAACCGGATCGCCAAACCGACGCTGATTATTCACGCCAAAGACGATCCGTTTATGGATCATCAGGTGATCCCGAAACCGGAAAGTC‑T‑CCCCC‑CGCAGGTGGAGTATCAACTGACTGAACATGGCGGTCATGTTGGCTTTATTGGCGGTACATTACTTCATCCGCAAATGTGGCT  >  NC_000913/3481627‑3481806

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: `‑`