BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,584,684	G→A	Q37Q (CAG→CAA)	yhhA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,584,684	0	G	A	87.5%	853.1 / 71.8	289	Q37Q (CAG→CAA)	yhhA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (18/18); new base A (132/121); total (150/139)
Fisher's exact test for biased strand distribution p-value = 8.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 309 total aligned reads displayed.

ACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAG  >  NC_000913/3584595‑3584774

aCGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACgg                                                                                         >  1:4625406/1‑93 (MQ=255)

aCGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACgg                                                                                         <  2:6755574/93‑1 (MQ=255)

 cGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGa                                                                                        <  2:2205/93‑1 (MQ=255)

  ggCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGAt                                                                                       <  2:5321867/93‑1 (MQ=255)

     aCTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCa                                                                                    <  2:2554324/93‑1 (MQ=255)

      cTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCaa                                                                                   <  2:2432264/93‑1 (MQ=255)

      cTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCaa                                                                                   <  1:6890658/93‑1 (MQ=255)

       tCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCaaa                                                                                  >  1:1024576/1‑93 (MQ=255)

          tGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TGCCCAGCCAGCAACGGATGCAAAccc                                                                               <  2:4645705/93‑1 (MQ=255)

          tGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAAccc                                                                               >  1:3472963/1‑93 (MQ=255)

          tGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAAccc                                                                               <  2:1699801/93‑1 (MQ=255)

          tGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAAccc                                                                               >  2:3477926/1‑93 (MQ=255)

             cGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGa                                                                            <  2:5164662/93‑1 (MQ=255)

             cGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGa                                                                            >  2:5934253/1‑93 (MQ=255)

             cGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGa                                                                            >  2:1376843/1‑93 (MQ=255)

              gTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGAt                                                                           >  1:1992836/1‑93 (MQ=255)

              gTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGAt                                                                           <  2:712775/93‑1 (MQ=255)

               tttGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATg                                                                          >  2:2872793/1‑93 (MQ=255)

                 tGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCa                                                                        <  1:4882165/93‑1 (MQ=255)

                     ggCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACt                                                                    >  2:5871570/1‑93 (MQ=255)

                     ggCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACt                                                                    <  2:5472918/93‑1 (MQ=255)

                      gCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTc                                                                   >  1:137789/1‑93 (MQ=255)

                        tttGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAg                                                                 <  2:4613980/93‑1 (MQ=255)

                        tttGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAg                                                                 >  2:4477389/1‑93 (MQ=255)

                             acaGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAAt                                                            >  2:4806836/1‑93 (MQ=255)

                             acaGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAAt                                                            >  1:3223953/1‑93 (MQ=255)

                             acaGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAAt                                                            >  1:1771048/1‑93 (MQ=255)

                               aGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATcc                                                          >  2:4906102/1‑93 (MQ=255)

                                gCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCa                                                         <  1:4658452/93‑1 (MQ=255)

                                gCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCa                                                         >  2:4015868/1‑93 (MQ=255)

                                  ccATTAATGCTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGc                                                       <  1:5832321/93‑1 (MQ=255)

                                   cATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCa                                                      >  2:5250443/1‑93 (MQ=255)

                                     ttAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaa                                                    >  2:1590837/1‑93 (MQ=255)

                                       aaTACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaaa                                                  >  1:6678525/1‑93 (MQ=255)

                                       aaTACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaaa                                                  >  2:3934726/1‑93 (MQ=255)

                                       aaTACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaaa                                                  <  1:3739233/93‑1 (MQ=255)

                                       aaTACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaaa                                                  <  1:3119306/93‑1 (MQ=255)

                                       aaTACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaaa                                                  >  2:7365896/1‑93 (MQ=255)

                                        aTACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAg                                                 >  2:108480/1‑93 (MQ=255)

                                          aCTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCATACTCAGCAAATCCAGCAAAAAggg                                               <  1:5386601/93‑1 (MQ=255)

                                          aCTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAggg                                               >  2:2302316/1‑93 (MQ=255)

                                           ctctGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                              <  1:6780170/93‑1 (MQ=255)

                                           ctctGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                              <  1:3581710/93‑1 (MQ=255)

                                              tGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGc                                           >  1:5537063/1‑93 (MQ=255)

                                               gAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGGAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGACGCt                                          <  1:6970642/93‑1 (MQ=255)

                                               gAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCt                                          >  1:212876/1‑93 (MQ=255)

                                                   aacCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAAt                                      >  2:7285965/1‑93 (MQ=255)

                                                    acCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATc                                     >  2:6658938/1‑93 (MQ=255)

                                                     cccTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCa                                    <  1:6007368/93‑1 (MQ=255)

                                                      ccTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAg                                   >  2:5592107/1‑93 (MQ=255)

                                                       cTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGc                                  >  2:3869829/1‑93 (MQ=255)

                                                       cTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGc                                  >  2:268792/1‑93 (MQ=255)

                                                       cTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGc                                  <  1:1624902/93‑1 (MQ=255)

                                                       cTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGc                                  <  1:1649694/93‑1 (MQ=255)

                                                        tAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCa                                 <  1:7143549/93‑1 (MQ=255)

                                                          aCCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                               <  2:6379002/93‑1 (MQ=255)

                                                          aCCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                               >  1:377284/1‑93 (MQ=255)

                                                          aCCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                               <  2:6099195/93‑1 (MQ=255)

                                                          aCCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                               <  1:450518/93‑1 (MQ=255)

                                                          aCCAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                               <  1:7355327/93‑1 (MQ=255)

                                                            cAGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTg                             >  2:3742364/1‑93 (MQ=255)

                                                             aGCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGa                            <  2:2388242/93‑1 (MQ=255)

                                                              gCCGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaa                           >  1:5275102/1‑93 (MQ=255)

                                                               gcGGGGTATCCGATTCCCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaa                           <  1:1139681/92‑1 (MQ=255)

                                                                cGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaaaa                         <  2:3125023/93‑1 (MQ=255)

                                                                cGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaaaa                         >  2:7178702/1‑93 (MQ=255)

                                                                cGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaaaa                         <  1:2831335/93‑1 (MQ=255)

                                                                cGGGGTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaaaa                         <  2:7018447/93‑1 (MQ=255)

                                                                 ggggTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAAc                        >  1:5654671/1‑93 (MQ=255)

                                                                   ggTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGc                      <  1:2308456/93‑1 (MQ=255)

                                                                   ggTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGc                      <  1:5749354/93‑1 (MQ=255)

                                                                    gTATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCa                     <  1:6475964/93‑1 (MQ=255)

                                                                     tATCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCaa                    <  1:5122651/93‑1 (MQ=255)

                                                                      aTCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCaaa                   <  1:2721760/93‑1 (MQ=255)

                                                                      aTCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCaaa                   <  1:5013245/93‑1 (MQ=255)

                                                                       tCAGA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAAc                  >  1:3333837/1‑93 (MQ=255)

                                                                          gA‑TCCCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCa               <  1:3897295/93‑1 (MQ=255)

                                                                          gA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCa               <  1:4034602/93‑1 (MQ=255)

                                                                          gA‑TACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCa               >  2:4329127/1‑93 (MQ=255)

                                                                           a‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAg              >  2:1981208/1‑93 (MQ=255)

                                                                             tACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGt             >  2:1498377/1‑93 (MQ=255)

                                                                             tACCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGt             <  2:6458122/93‑1 (MQ=255)

                                                                              aCCCAGCCAGCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGtt            <  1:610217/93‑1 (MQ=255)

                                                                                ccagccagCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGc          <  1:4809440/93‑1 (MQ=255)

                                                                                ccagccagCAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGc          <  1:6560908/93‑1 (MQ=255)

                                                                                 cagccagcAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCa         >  2:797544/1‑93 (MQ=255)

                                                                                 cagccagcAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCa         >  2:5666447/1‑93 (MQ=255)

                                                                                   gccagcAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAc       <  2:2120016/93‑1 (MQ=255)

                                                                                   gccagcAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAc       <  1:3041788/93‑1 (MQ=255)

                                                                                   gccagcAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAc       <  1:2260059/93‑1 (MQ=255)

                                                                                    ccagcAACGGATGCAACCCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAca      <  1:6145926/93‑1 (MQ=255)

                                                                                    ccagcAACGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAca      >  1:7425537/1‑93 (MQ=255)

                                                                                     cagcaaCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcag     >  2:2037883/1‑93 (MQ=255)

                                                                                     cagcaaCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcag     <  1:845825/93‑1 (MQ=255)

                                                                                      agcagcGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcagc    >  2:6956829/1‑93 (MQ=255)

                                                                                      agcaacGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcagc    <  2:1311602/93‑1 (MQ=255)

                                                                                       gcagcGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcagca   <  1:7243177/93‑1 (MQ=255)

                                                                                       gcaacGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcagca   >  2:3542353/1‑93 (MQ=255)

                                                                                       gcaacGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcagca   <  1:6352943/93‑1 (MQ=255)

                                                                                        cagcGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAcagcag  <  1:2609048/93‑1 (MQ=255)

ACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGA‑TACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAG  >  NC_000913/3584595‑3584774

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