BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,116,479	A→G	A150A (GCA→GCG)	yceA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,116,479	0	A	G	98.7%	464.8 / ‑4.2	148	A150A (GCA→GCG)	yceA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (64/82); total (66/82)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 154 total aligned reads displayed.

GTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGAGCA  >  NC_000913/1116391‑1116568

gTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTAt                                                                                       >  2:2711508/1‑93 (MQ=255)

gTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCACGATTGACGATCCCGATGCGCTAt                                                                                       >  2:372310/1‑93 (MQ=255)

 tATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTAtt                                                                                      <  1:782506/93‑1 (MQ=255)

   ttGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTa                                                                                    <  1:3400327/93‑1 (MQ=255)

   ttGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTa                                                                                    <  2:3673828/93‑1 (MQ=255)

   ttGACGATCCTCACTTTGATGCCAGCAATGGTGGTGAGTATCTGCAAGCGGCGGAAGTGTACGCCACGCTTGATTATCCCGTTGCGCTATTTa                                                                                    >  1:3615682/1‑93 (MQ=255)

     gACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATc                                                                                  >  1:3110668/1‑93 (MQ=255)

       cGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGa                                                                                >  2:1623013/1‑93 (MQ=255)

        gATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGAc                                                                               >  1:2333802/1‑93 (MQ=255)

         aTCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACa                                                                              >  1:2086124/1‑93 (MQ=255)

         aTCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACa                                                                              >  1:2086123/1‑93 (MQ=255)

         aTCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGGGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACa                                                                              >  2:4523626/1‑93 (MQ=255)

            cTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGc                                                                           >  1:1412425/1‑93 (MQ=255)

            cTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGc                                                                           <  1:3080525/93‑1 (MQ=255)

              cACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGt                                                                         <  2:623457/93‑1 (MQ=255)

               aCTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTa                                                                        <  1:1818700/93‑1 (MQ=255)

                cTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTaa                                                                       >  1:3840826/1‑93 (MQ=255)

                 tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTagc                                                                      >  1:3685365/1‑91 (MQ=255)

                 tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAAc                                                                      >  1:1982709/1‑93 (MQ=255)

                 tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAAc                                                                      >  2:705875/1‑93 (MQ=255)

                 tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGACGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAAc                                                                      >  1:709555/1‑93 (MQ=255)

                   tGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCa                                                                    <  1:1187366/93‑1 (MQ=255)

                    gATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCAc                                                                   <  2:4122016/93‑1 (MQ=255)

                     aTGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACt                                                                  <  1:1673097/93‑1 (MQ=255)

                      tGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTa                                                                 <  1:585244/93‑1 (MQ=255)

                       gCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTAt                                                                <  1:1030399/93‑1 (MQ=255)

                         cAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTACCCACTATGa                                                              >  2:2398236/1‑93 (MQ=255)

                         cAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGa                                                              >  1:2978052/1‑93 (MQ=255)

                         cAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGa                                                              <  1:4313505/93‑1 (MQ=255)

                         cAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGa                                                              >  1:4420455/1‑93 (MQ=255)

                          aGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAg                                                             <  1:537796/93‑1 (MQ=255)

                          aGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAg                                                             >  2:3766604/1‑93 (MQ=255)

                          aGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAg                                                             <  1:1216901/93‑1 (MQ=255)

                           gCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGt                                                            >  2:1138455/1‑93 (MQ=255)

                           gCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGt                                                            >  1:888009/1‑93 (MQ=255)

                           gCAATGTTGGGGAGGATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGt                                                            <  1:3353646/93‑1 (MQ=255)

                            gaaTGTTGGTGAGTATCTGCAAGCGGAGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTa                                                           <  2:1859903/92‑1 (MQ=255)

                               tGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGa                                                        >  2:1190295/1‑93 (MQ=255)

                               tGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGa                                                        <  2:3809513/93‑1 (MQ=255)

                                gTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGCAACCACTATGAGTATGaa                                                       >  1:816693/1‑93 (MQ=255)

                                 ttGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAg                                                      <  2:1966128/93‑1 (MQ=255)

                                    gTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTgg                                                   <  1:3348933/93‑1 (MQ=255)

                                     tGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTggg                                                  >  2:1217005/1‑93 (MQ=255)

                                      gAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTgggg                                                 <  1:189033/93‑1 (MQ=255)

                                      gAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTgggg                                                 >  1:3491259/1‑93 (MQ=255)

                                       aGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGc                                                <  1:2299381/93‑1 (MQ=255)

                                       aGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGc                                                >  1:3190151/1‑93 (MQ=255)

                                           tCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACtt                                            <  2:3496499/93‑1 (MQ=255)

                                            cTGCAAGCGGCGGAAGTGAACGCTATGCTTGACGATCGCGATGCGCTATTTATCGACAGGCGTAACCACTATGAGTATGAAGTGGGGCACttt                                           <  1:1298920/93‑1 (MQ=255)

                                              gCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGa                                         >  2:4326236/1‑93 (MQ=255)

                                                aaGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaaa                                       <  1:3345486/93‑1 (MQ=255)

                                                 aGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaaaa                                      <  2:1783721/93‑1 (MQ=255)

                                                 aGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaaaa                                      <  1:1161409/93‑1 (MQ=255)

                                                   cggcggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcg                                    <  1:1041843/93‑1 (MQ=255)

                                                    ggcggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgc                                   >  1:1403703/1‑93 (MQ=255)

                                                     gcggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgcg                                  >  2:482125/1‑93 (MQ=255)

                                                     gcggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATCAAGTGGGGCACTTTGAAAAcgcg                                  >  2:2590182/1‑93 (MQ=255)

                                                      cggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACgcgc                                 <  1:791570/93‑1 (MQ=255)

                                                      cggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACgcgc                                 <  2:3202946/93‑1 (MQ=255)

                                                      cggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACgcgc                                 >  2:1848223/1‑93 (MQ=255)

                                                       ggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCt                                >  1:4118455/1‑93 (MQ=255)

                                                       ggAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCt                                >  2:2773948/1‑93 (MQ=255)

                                                        gAAGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTc                               >  2:2195647/1‑93 (MQ=255)

                                                         aaGTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCg                              >  1:4719928/1‑93 (MQ=255)

                                                           gTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGaa                            >  2:4373988/1‑93 (MQ=255)

                                                           gTGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGaa                            <  2:2134260/93‑1 (MQ=255)

                                                            tGAACGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGaaa                           <  1:28157/93‑1 (MQ=255)

                                                              aaCGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAAtt                         >  2:1225646/1‑93 (MQ=255)

                                                                cGCCATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTcc                       <  1:4402199/93‑1 (MQ=255)

                                                                 cgcATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCg                      <  1:1674601/91‑1 (MQ=255)

                                                                  ttaTGCTTGATGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACCTTGAAAACGCGCTCGAAATTCCgg                     >  1:4194930/3‑93 (MQ=255)

                                                                  ccATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCgg                     <  1:2545060/93‑1 (MQ=255)

                                                                  ccATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCgg                     <  1:2545184/93‑1 (MQ=255)

                                                                  ccATGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCgg                     >  2:2679991/1‑93 (MQ=255)

                                                                    aTGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa                   >  1:1710479/1‑93 (MQ=255)

                                                                    aTGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa                   >  2:414247/1‑93 (MQ=255)

                                                                    aTGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa                   <  2:1933920/93‑1 (MQ=255)

                                                                     tGCTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAg                  >  1:3595467/1‑93 (MQ=255)

                                                                       cTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGAt                <  2:4687110/93‑1 (MQ=255)

                                                                       cTTGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGAt                <  2:3049993/93‑1 (MQ=255)

                                                                         tGACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAc              <  1:1218496/93‑1 (MQ=255)

                                                                          gACGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAcc             <  2:3248408/93‑1 (MQ=255)

                                                                           aCGATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCt            <  2:3929599/93‑1 (MQ=255)

                                                                            cGATCCCGATGCGCTATTGATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCtt           <  1:2007657/93‑1 (MQ=255)

                                                                             gATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTc          <  2:2447470/93‑1 (MQ=255)

                                                                             gATCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTc          <  1:3814382/93‑1 (MQ=255)

                                                                              aTCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTcc         >  1:922970/1‑93 (MQ=255)

                                                                               tCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCg        >  1:407926/1‑93 (MQ=255)

                                                                               tCCCGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCg        >  1:4000605/1‑93 (MQ=255)

                                                                                 ccGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTg      <  1:2142966/93‑1 (MQ=255)

                                                                                  cGATGCGCTATTTATCGACATTCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATCCCTTCCGTGa     <  2:4586397/93‑1 (MQ=255)

                                                                                  cGATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGa     <  2:3001072/93‑1 (MQ=255)

                                                                                   gATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGag    >  2:3719872/1‑93 (MQ=255)

                                                                                   gATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGag    >  1:1282164/1‑93 (MQ=255)

                                                                                   gATGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGag    >  1:426623/1‑93 (MQ=255)

                                                                                    aTGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGagc   <  2:2100812/93‑1 (MQ=255)

                                                                                    aTGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGagc   >  2:1462625/1‑93 (MQ=255)

                                                                                     tGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGagca  >  2:4299215/1‑93 (MQ=255)

                                                                                     tGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGagca  <  1:2144183/93‑1 (MQ=255)

                                                                                     tGCGCTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGagca  <  1:3128352/93‑1 (MQ=255)

GTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATACCTTCCGTGAGCA  >  NC_000913/1116391‑1116568

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