BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,262,748	A→G	T286A (ACC→GCC)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,262,748	0	A	G	87.1%	838.6 / 104.7	342	T286A (ACC→GCC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand): ref base A (27/17); new base G (134/164); total (161/181)
Fisher's exact test for biased strand distribution p-value = 5.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 343 total aligned reads displayed.

GGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGATGACTCCTCC  >  NC_000913/4262657‑4262840

 gTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGcc                                                                                            >  2:4609669/1‑93 (MQ=255)

  tACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCg                                                                                           >  2:250323/1‑93 (MQ=255)

   aCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCgg                                                                                          <  1:4702208/93‑1 (MQ=255)

   aCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCgg                                                                                          >  1:1989946/1‑93 (MQ=255)

   aCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCgg                                                                                          >  2:1745375/1‑93 (MQ=255)

    cTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCggg                                                                                         >  1:6106001/1‑93 (MQ=255)

     ttATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGc                                                                                        <  2:4967185/93‑1 (MQ=255)

     ttATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGc                                                                                        >  1:5875884/1‑93 (MQ=255)

     ttATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGc                                                                                        >  2:4137376/1‑93 (MQ=255)

     ttATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGc                                                                                        >  1:4563679/1‑93 (MQ=255)

        tcttctCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGc                                                                                     >  2:2282024/1‑93 (MQ=255)

        tcttctCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGc                                                                                     <  2:7526106/93‑1 (MQ=255)

        tcttctCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGc                                                                                     <  1:5498933/93‑1 (MQ=255)

        tcttctCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGc                                                                                     >  1:1412076/1‑93 (MQ=255)

          ttctCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTc                                                                                   >  1:2148263/1‑93 (MQ=255)

           tctcGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCg                                                                                  >  2:3759531/1‑93 (MQ=255)

             tcGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGAt                                                                                >  1:7230238/1‑93 (MQ=255)

             tcGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGAt                                                                                <  2:4502897/93‑1 (MQ=255)

               gCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGa                                                                              <  2:2710963/93‑1 (MQ=255)

               gCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGa                                                                              >  2:4289543/1‑93 (MQ=255)

                cTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGAc                                                                             <  1:2480841/93‑1 (MQ=255)

                 tGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACg                                                                            >  1:2550267/1‑93 (MQ=255)

                   gagaTGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                                                          >  2:950705/1‑93 (MQ=255)

                     gaTGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGa                                                                        <  1:1192748/93‑1 (MQ=255)

                     gaTGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGa                                                                        >  1:3405757/1‑93 (MQ=255)

                      aTGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGAc                                                                       <  2:6833551/93‑1 (MQ=255)

                       tGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACt                                                                      >  2:5171281/1‑93 (MQ=255)

                         ccATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTgg                                                                    <  1:5752162/93‑1 (MQ=255)

                          cATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTggg                                                                   <  2:3668385/93‑1 (MQ=255)

                           aTCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgc                                                                  >  1:711132/1‑93 (MQ=255)

                           aTCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATCACGAAGACTGGgc                                                                  >  1:3202275/1‑93 (MQ=255)

                            tCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcg                                                                 >  1:238455/1‑93 (MQ=255)

                            tCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcg                                                                 >  2:3304533/1‑93 (MQ=255)

                              aGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcg                                                               >  1:3705993/1‑93 (MQ=255)

                              aGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcg                                                               >  1:2981684/1‑93 (MQ=255)

                              aGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcg                                                               >  2:912175/1‑93 (MQ=255)

                              aGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcg                                                               >  2:7283763/1‑93 (MQ=255)

                              aGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcg                                                               >  1:262434/1‑93 (MQ=255)

                              aGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcg                                                               >  1:1564082/1‑93 (MQ=255)

                               gAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcgc                                                              >  1:6478404/1‑93 (MQ=255)

                               gAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcgc                                                              >  1:4461803/1‑93 (MQ=255)

                               gAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGgcgcgc                                                              >  2:1628962/1‑93 (MQ=255)

                               gAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcgc                                                              >  2:6575631/1‑93 (MQ=255)

                                  cAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCAtt                                                           >  2:239569/1‑93 (MQ=255)

                                   aGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCAttt                                                          >  1:6681825/1‑93 (MQ=255)

                                    gATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                                                         >  1:693178/1‑93 (MQ=255)

                                    gATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                                                         <  2:7001451/93‑1 (MQ=255)

                                    gATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                                                         <  2:6150396/93‑1 (MQ=255)

                                    gATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                                                         >  2:577569/1‑93 (MQ=255)

                                     aTCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTcc                                                        >  1:6400583/1‑93 (MQ=255)

                                      tCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                                                       <  2:2901864/93‑1 (MQ=255)

                                      tCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                                                       <  1:54136/93‑1 (MQ=255)

                                         tgatgCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCa                                                    >  1:3506361/1‑93 (MQ=255)

                                          gatgCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCAc                                                   >  2:1313652/1‑93 (MQ=255)

                                             gCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCAt                                                <  1:7373157/93‑1 (MQ=255)

                                              cGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATg                                               <  1:187031/93‑1 (MQ=255)

                                               gTTCTCTGGAGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATgg                                              <  2:1006312/93‑1 (MQ=255)

                                                ttCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATggg                                             >  1:7007762/1‑93 (MQ=255)

                                                ttCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATggg                                             >  1:5215079/1‑93 (MQ=255)

                                                  ctctGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa                                           >  1:5351852/1‑93 (MQ=255)

                                                   tctGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt                                          >  2:457325/1‑93 (MQ=255)

                                                   tctGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt                                          <  1:4456439/93‑1 (MQ=255)

                                                   tctGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt                                          >  1:6624275/1‑93 (MQ=255)

                                                   tctGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt                                          >  2:3507969/1‑93 (MQ=255)

                                                     tGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt                                        >  2:4866046/1‑93 (MQ=255)

                                                     tGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt                                        >  1:2879181/1‑93 (MQ=255)

                                                      ggCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtttt                                       <  2:6384609/93‑1 (MQ=255)

                                                        cGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGc                                     <  1:4797760/93‑1 (MQ=255)

                                                        cGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGc                                     <  2:357931/93‑1 (MQ=255)

                                                        cGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGc                                     <  2:4398720/93‑1 (MQ=255)

                                                         gTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCt                                    <  2:1859113/93‑1 (MQ=255)

                                                          tCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTc                                   <  2:1463074/93‑1 (MQ=255)

                                                          tCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTc                                   >  2:3340155/1‑93 (MQ=255)

                                                          tCGCTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTc                                   <  1:5600737/93‑1 (MQ=255)

                                                             cTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGaa                                >  1:5000070/1‑93 (MQ=255)

                                                             cTGTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGaa                                >  1:7516153/1‑93 (MQ=255)

                                                               gTCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGaaaa                              >  2:3862831/1‑93 (MQ=255)

                                                                tCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGaaaaa                             <  1:3296047/93‑1 (MQ=255)

                                                                tCGCGCGTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGaaaaa                             <  2:5385415/93‑1 (MQ=255)

                                                                 cgcgcgTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAAc                            <  1:846014/93‑1 (MQ=255)

                                                                    gcgTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCa                         >  2:2204945/1‑93 (MQ=255)

                                                                     cgTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCaa                        <  2:3001992/93‑1 (MQ=255)

                                                                      gTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAat                       <  2:6187614/93‑1 (MQ=255)

                                                                      gTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAat                       >  1:1105329/1‑93 (MQ=255)

                                                                      gTTGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAat                       >  1:1012420/1‑93 (MQ=255)

                                                                       ttGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAata                      >  1:406474/1‑93 (MQ=255)

                                                                        tGACCAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAatat                     <  2:3993303/93‑1 (MQ=255)

                                                                           ccAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCta                  <  1:5815850/93‑1 (MQ=255)

                                                                           ccAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCta                  <  1:4021709/93‑1 (MQ=255)

                                                                           ccAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCta                  <  1:388377/93‑1 (MQ=255)

                                                                            cAGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCtat                 <  2:1565975/93‑1 (MQ=255)

                                                                             aGACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCtata                <  1:5960919/93‑1 (MQ=255)

                                                                              gACTAAAATCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTatat               <  1:6586147/93‑1 (MQ=255)

                                                                                  aaaaTCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGAt           <  2:5268054/93‑1 (MQ=255)

                                                                                  aaaaTCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGAt           <  1:739918/93‑1 (MQ=255)

                                                                                   aaaTCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGATg          <  2:812540/93‑1 (MQ=255)

                                                                                    aaTCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGATGa         <  2:5612937/93‑1 (MQ=255)

                                                                                    aaTCCGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGATGa         <  1:3459365/93‑1 (MQ=255)

                                                                                        cGTGCCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGATGActcc     <  1:2842238/93‑1 (MQ=255)

                                                                                          tACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATAGGTATTTTGCTCGAAAAACGCAATATCTATATCGATGActcctc   <  2:1349494/93‑1 (MQ=255)

GGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTCGAAAAACGCAATATCTATATCGATGACTCCTCC  >  NC_000913/4262657‑4262840

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