BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,971	T→C	G259G (GGT→GGC)	gntT →	gluconate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,544,971	0	T	C	96.7%	536.4 / 4.4	180	G259G (GGT→GGC)	gntT	gluconate permease
Reads supporting (aligned to +/- strand): ref base T (3/3); new base C (87/86); total (91/89)
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.47e-01

Only 100 of 190 total aligned reads displayed.

GGTGCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGA  >  NC_000913/3544881‑3545060

ggTGCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGa                                                                                         <  1:4653893/93‑1 (MQ=255)

   gCCGGTAGTCGGGATGGCGATGCTTGCGCTGGCCGAATGGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGAccc                                                                                      <  2:2967624/93‑1 (MQ=255)

   gCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGAccc                                                                                      <  1:645430/93‑1 (MQ=255)

    ccGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCg                                                                                     >  1:185759/1‑93 (MQ=255)

      ggTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGt                                                                                   >  1:4732030/1‑93 (MQ=255)

          gtaCTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATg                                                                               >  2:3519186/1‑93 (MQ=255)

           taCTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATgg                                                                              >  2:4410958/1‑93 (MQ=255)

            aCTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGc                                                                             <  1:2245399/93‑1 (MQ=255)

            aCTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGc                                                                             >  2:3602685/1‑93 (MQ=255)

             cTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCa                                                                            >  2:797722/1‑93 (MQ=255)

              tGATGGCGATGCGTGCGATTGCCTAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCaa                                                                           <  2:3392558/93‑1 (MQ=255)

              tGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCaa                                                                           >  1:3144163/1‑93 (MQ=255)

               gATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAAc                                                                          >  2:1739597/1‑93 (MQ=255)

                 tGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGc                                                                        >  2:2879575/1‑93 (MQ=255)

                 gggCGAGGCGTGCGATTGCCAAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGc                                                                        <  1:1061095/92‑1 (MQ=255)

                  ggCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCt                                                                       >  2:338682/1‑93 (MQ=255)

                   gCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTg                                                                      >  1:1621789/1‑93 (MQ=255)

                    ggATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGa                                                                     <  1:1660743/92‑1 (MQ=255)

                     gATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGAt                                                                    >  2:4017499/1‑93 (MQ=255)

                     gATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGAt                                                                    >  2:3369768/1‑93 (MQ=255)

                     gATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTCAt                                                                    >  2:341254/1‑93 (MQ=255)

                      aTGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGAtt                                                                   >  2:1099643/1‑93 (MQ=255)

                      aTGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGAtt                                                                   <  1:2592629/93‑1 (MQ=255)

                       tgcgtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTg                                                                  >  2:1414205/1‑93 (MQ=255)

                        gcgtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGc                                                                 >  2:412713/1‑93 (MQ=255)

                         cgtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGcc                                                                >  2:1973394/1‑93 (MQ=255)

                          gtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGACACGCTGATTGCCg                                                               >  2:2122174/1‑93 (MQ=255)

                           tgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGt                                                              >  1:412200/1‑93 (MQ=255)

                            gcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTg                                                             >  2:323687/1‑93 (MQ=255)

                            gcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTg                                                             >  2:3876570/1‑93 (MQ=255)

                               aTTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTg                                                          >  1:3662541/1‑93 (MQ=255)

                                ttGCCGAAATGATCCTGCCGACAGGTCACGCTTTACTGCCGGTAGCGGAGTTCCTCGGCGACCCGGCAATGGGAACGCTGATTGGCGTGCTGa                                                         >  2:2679508/1‑93 (MQ=255)

                                 tGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGAt                                                        <  1:2683551/93‑1 (MQ=255)

                                  gCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGAtt                                                       >  2:3297164/1‑93 (MQ=255)

                                    cGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGc                                                     <  2:1841873/93‑1 (MQ=255)

                                     gAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCg                                                    >  2:1589695/1‑93 (MQ=255)

                                     gAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCg                                                    >  2:1047933/1‑93 (MQ=255)

                                      aaaTGATGCTGCCGAAAGGTCACGTTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  1:352390/93‑1 (MQ=255)

                                      aaaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  2:164851/93‑1 (MQ=255)

                                      aaaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  2:479964/93‑1 (MQ=255)

                                      aaaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  2:2435697/93‑1 (MQ=255)

                                      aaaTGACCCTGCCGAAAGGTAACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  1:352389/93‑1 (MQ=255)

                                       aaTGATCTTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGAt                                                  <  1:4454914/93‑1 (MQ=255)

                                        aTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATg                                                 >  1:2077753/1‑93 (MQ=255)

                                         tGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGt                                                >  1:2867312/1‑93 (MQ=255)

                                         tGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGAGATGt                                                >  2:181925/1‑93 (MQ=255)

                                          gacccTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGtt                                               <  1:3753968/90‑1 (MQ=255)

                                          gATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGtt                                               >  2:2937965/1‑93 (MQ=255)

                                          gATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGtt                                               >  1:4473231/1‑93 (MQ=255)

                                          gATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGtt                                               >  2:3603803/1‑93 (MQ=255)

                                           aTCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTc                                              <  1:4522945/93‑1 (MQ=255)

                                           aTCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTc                                              >  2:83001/1‑93 (MQ=255)

                                           aTCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTc                                              >  1:110262/1‑93 (MQ=255)

                                           aTCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTc                                              >  1:1280202/1‑93 (MQ=255)

                                              cTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCAcc                                           >  2:4510774/1‑93 (MQ=255)

                                                gCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTTACCtt                                         >  2:59182/1‑93 (MQ=255)

                                                 ccGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCttt                                        <  2:2199367/93‑1 (MQ=255)

                                                    aaaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGt                                     >  2:2621723/1‑93 (MQ=255)

                                                     aaGGTCGCTCATACCTGCCGGTAGCGGAGTAACTAGGCGACCCGGTAGTAGATACGCTGATTGCCGTGCTAATTGAGATATTCACCTTTGGTc                                    <  2:4202401/93‑1 (MQ=255)

                                                      aGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCt                                   >  1:4717223/1‑93 (MQ=255)

                                                       ggTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTa                                  <  2:2717788/93‑1 (MQ=255)

                                                         tCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTaaa                                >  2:537505/1‑93 (MQ=255)

                                                          cACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTTATTGCGATGTTCACCTTTGGTCTaagc                               <  2:4655515/93‑3 (MQ=255)

                                                          cACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAc                               >  2:3844506/1‑93 (MQ=255)

                                                          cACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAc                               <  1:4155843/93‑1 (MQ=255)

                                                           aCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAcc                              >  2:2167941/1‑93 (MQ=255)

                                                           aCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAcc                              >  2:3689630/1‑93 (MQ=255)

                                                           aCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAcc                              >  1:299506/1‑93 (MQ=255)

                                                            cGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCg                             <  1:2787188/93‑1 (MQ=255)

                                                             gCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGt                            >  2:2521402/1‑93 (MQ=255)

                                                             gCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGt                            <  2:3785824/93‑1 (MQ=255)

                                                             gCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGt                            <  1:687110/93‑1 (MQ=255)

                                                              cTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTg                           <  2:1073044/93‑1 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgg                          <  1:4022565/93‑1 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgg                          >  2:3079041/1‑93 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgg                          <  2:115868/93‑1 (MQ=255)

                                                                 tCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTc                        <  2:810492/93‑1 (MQ=255)

                                                                 tCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGAAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTc                        >  2:1469069/1‑93 (MQ=255)

                                                                   cTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACcgtcgtcgt                      >  2:2665175/1‑93 (MQ=255)

                                                                    tGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGtt                     >  1:2289419/1‑93 (MQ=255)

                                                                       cGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCaa                  <  2:2532729/93‑1 (MQ=255)

                                                                        ggTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGGTCAAt                 >  2:1727236/1‑93 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                >  1:3605793/1‑93 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                <  2:3558492/93‑1 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                >  2:3591027/1‑93 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                >  1:773841/1‑93 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                <  2:4521390/93‑1 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                <  1:2367737/93‑1 (MQ=255)

                                                                               gAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAg          <  1:1420751/93‑1 (MQ=255)

                                                                               gAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGGGCTGATTGCGATGGTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGGTCAg          >  1:1312421/1‑93 (MQ=255)

                                                                               gAGTGCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAg          <  2:257042/93‑1 (MQ=255)

                                                                                ggtgccTCGTCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGa         <  2:1553511/89‑1 (MQ=255)

                                                                                 gTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGAt        <  1:2179422/93‑1 (MQ=255)

                                                                                 gTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGAt        <  1:433384/93‑1 (MQ=255)

                                                                                 gTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGAt        >  2:1009649/1‑93 (MQ=255)

                                                                                   tCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTa      <  2:4701861/93‑1 (MQ=255)

                                                                                    tcTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTaa     <  1:4103331/92‑1 (MQ=255)

                                                                                    ccTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTaa     <  1:2044873/93‑1 (MQ=255)

                                                                                      tCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACg   >  2:2465444/1‑93 (MQ=255)

                                                                                       cGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGa  >  1:873556/1‑93 (MQ=255)

GGTGCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGA  >  NC_000913/3544881‑3545060

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