BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,566	A→G	G515G (GGT→GGC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,566	0	A	G	99.5%	623.3 / NA	195	G515G (GGT→GGC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (96/98); total (96/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAG  >  NC_000913/1594476‑1594655

aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                                                                         <  1:3074850/93‑1 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                     >  2:2656409/1‑93 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                     >  1:1779331/1‑93 (MQ=255)

     agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTAc                                                                                    <  2:5913237/93‑1 (MQ=255)

     agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTAc                                                                                    <  1:2925459/93‑1 (MQ=255)

       agaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTg                                                                                  <  2:1832280/93‑1 (MQ=255)

        taCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 <  1:6047462/92‑1 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  2:1188771/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  2:4216800/1‑93 (MQ=255)

         aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                                                                >  1:1954808/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                               <  1:330020/93‑1 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                             <  2:4918717/93‑1 (MQ=255)

              ccGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                                                                           <  2:2104539/93‑1 (MQ=255)

              ccGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                                                                           >  1:5278817/1‑93 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                                                                          >  2:2820952/1‑93 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                         >  2:4262191/1‑93 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                         >  2:2237469/1‑93 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                         <  2:1502581/93‑1 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                        <  1:3710099/93‑1 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                        <  1:6207843/93‑1 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                        >  1:4680304/1‑93 (MQ=255)

                 ttGAAATATTAACTATCCCCCTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                        <  1:5218371/93‑1 (MQ=255)

                  tGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                                                                       >  1:4535978/1‑93 (MQ=255)

                   gAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                                                                      <  1:3550681/93‑1 (MQ=255)

                   gAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                                                                      >  2:3217814/1‑93 (MQ=255)

                    aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                                                                     >  1:5556301/1‑93 (MQ=255)

                     aaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAcc                                                                    <  1:992639/93‑1 (MQ=255)

                     aaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAcc                                                                    <  2:764703/93‑1 (MQ=255)

                             cTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaa                                                            >  2:634471/1‑93 (MQ=255)

                             cTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaa                                                            >  1:130030/1‑93 (MQ=255)

                              tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                           >  2:550141/1‑93 (MQ=255)

                              tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                           <  2:5628300/93‑1 (MQ=255)

                              tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                           >  1:5490209/1‑93 (MQ=255)

                                tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt                                                         <  2:1766946/93‑1 (MQ=255)

                                 ccccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt                                                        >  2:5720601/1‑93 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       >  1:6154521/1‑93 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       >  1:3651710/1‑93 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                      <  1:5050761/93‑1 (MQ=255)

                                    cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                     >  2:465010/1‑93 (MQ=255)

                                    cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                     <  1:2717459/93‑1 (MQ=255)

                                     ttACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcca                                                    >  2:1642458/1‑92 (MQ=255)

                                     ttACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg                                                    <  1:1533847/93‑1 (MQ=255)

                                      tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa                                                   <  2:3053548/93‑1 (MQ=255)

                                      tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa                                                   >  1:218187/1‑93 (MQ=255)

                                       aCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                  <  1:5431537/93‑1 (MQ=255)

                                            gggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                             >  2:4987883/2‑93 (MQ=255)

                                            gggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                             <  2:1509487/92‑1 (MQ=255)

                                              gCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATa                                           <  2:1711831/93‑1 (MQ=255)

                                               cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAtcc                                          <  2:4647229/93‑3 (MQ=255)

                                                aTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                         >  2:3442080/1‑93 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                        <  2:5200080/93‑1 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                        >  2:1472087/1‑93 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                        >  2:742024/1‑93 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                        <  2:4060111/93‑1 (MQ=255)

                                                  ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcccc                                       >  1:2309559/1‑92 (MQ=255)

                                                   tcctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                      <  2:5566369/93‑1 (MQ=255)

                                                   tcctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                      >  2:2996894/1‑93 (MQ=255)

                                                    cctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                     >  2:3718593/1‑93 (MQ=255)

                                                     ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                    >  1:5139007/1‑93 (MQ=255)

                                                      tccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                   >  2:135831/1‑93 (MQ=255)

                                                      tccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                   >  1:2409374/1‑93 (MQ=255)

                                                      tccGATAAATCCATAGCCCGTTGAATTAAGGGGGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                   >  2:389695/1‑93 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATCCCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 <  1:5776927/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 >  2:476936/1‑93 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 <  2:2409611/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 >  2:6276803/1‑93 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCAACAGTAt                                 <  1:935212/93‑1 (MQ=255)

                                                         gATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                >  1:3066191/1‑93 (MQ=255)

                                                         gATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                <  1:4976314/93‑1 (MQ=255)

                                                          aTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                               >  1:2372415/1‑93 (MQ=255)

                                                          aTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                               >  2:3160063/1‑93 (MQ=255)

                                                          aTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                               >  1:4063008/1‑93 (MQ=255)

                                                           tAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGa                              >  2:2042031/1‑93 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                           >  1:4417946/1‑93 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                           >  1:1683308/1‑93 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                           <  2:2172845/93‑1 (MQ=255)

                                                               tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                          <  1:1571806/93‑1 (MQ=255)

                                                                ccATAGCCCGTTTAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                         <  2:3273983/93‑1 (MQ=255)

                                                                ccATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                         <  1:1020642/93‑1 (MQ=255)

                                                                ccATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                         >  2:2448907/1‑93 (MQ=255)

                                                                ccATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                         <  2:262556/93‑1 (MQ=255)

                                                                 cATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                        <  2:356161/93‑1 (MQ=255)

                                                                 cATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                        <  2:3951235/93‑1 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                    <  2:908172/93‑1 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                    >  1:1401473/1‑93 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                    <  2:378788/93‑1 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGca                >  2:5337350/1‑93 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGca                >  2:1929163/1‑93 (MQ=255)

                                                                           tGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc              <  1:3760228/93‑1 (MQ=255)

                                                                           tGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc              <  1:4000119/93‑1 (MQ=255)

                                                                            gAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc             >  1:5932677/1‑93 (MQ=255)

                                                                            gAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc             >  2:5838549/1‑93 (MQ=255)

                                                                             aaTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca            >  1:2552880/1‑93 (MQ=255)

                                                                              aTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccat           >  2:392133/1‑93 (MQ=255)

                                                                               ttAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccatc          <  1:583980/93‑1 (MQ=255)

                                                                               ttAAGGGTGGCGCCATTACTGACATTGACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccatc          <  2:44096/93‑1 (MQ=255)

                                                                                tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa         <  1:337136/93‑1 (MQ=255)

                                                                                 aaGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg        <  1:6214674/93‑1 (MQ=255)

                                                                                  aGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg       <  1:4798255/93‑1 (MQ=255)

                                                                                       ggCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg  <  1:4269351/93‑1 (MQ=255)

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAG  >  NC_000913/1594476‑1594655

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