BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,207,023	A→G	D346G (GAC→GGC)	yehP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,207,023	0	A	G	99.5%	656.1 / ‑4.7	203	D346G (GAC→GGC)	yehP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (114/87); total (114/89)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 206 total aligned reads displayed.

TACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATC  >  NC_000913/2206933‑2207112

 aCTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCcg                                                                                        >  1:565020/1‑93 (MQ=255)

  cTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCcgc                                                                                       >  2:1200805/1‑93 (MQ=255)

      cGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATa                                                                                   <  1:3955043/93‑1 (MQ=255)

      cGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATa                                                                                   >  2:4882009/1‑93 (MQ=255)

      cGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATa                                                                                   <  2:145131/93‑1 (MQ=255)

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATAc                                                                                  <  2:2053360/93‑1 (MQ=255)

        cATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACg                                                                                 <  1:5911376/93‑1 (MQ=255)

         aTCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACgg                                                                                >  2:5164388/1‑93 (MQ=255)

          tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGc                                                                               >  1:5168857/1‑93 (MQ=255)

          tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGc                                                                               >  1:242198/1‑93 (MQ=255)

            agatGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGccc                                                                             <  1:2120637/90‑1 (MQ=255)

            aGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGccc                                                                             >  1:1770816/1‑93 (MQ=255)

            aGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGccc                                                                             <  2:4960489/93‑1 (MQ=255)

            aGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGccc                                                                             <  2:2297811/93‑1 (MQ=255)

             ggTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCa                                                                            >  2:1857019/1‑93 (MQ=255)

             ggTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCa                                                                            >  1:877116/1‑93 (MQ=255)

               tGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAgg                                                                          >  2:5085883/1‑93 (MQ=255)

                gAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgc                                                                         >  2:2187672/1‑93 (MQ=255)

                 aaaaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgcg                                                                        <  2:905139/93‑1 (MQ=255)

                  aaaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgcgc                                                                       >  2:2973195/1‑93 (MQ=255)

                   aaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCt                                                                      <  1:1085837/93‑1 (MQ=255)

                   aaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCt                                                                      <  1:561063/93‑1 (MQ=255)

                   aaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCt                                                                      >  2:4428432/1‑93 (MQ=255)

                    aaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTg                                                                     >  2:546860/1‑93 (MQ=255)

                     aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTgg                                                                    <  2:5316197/93‑1 (MQ=255)

                      gtgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGt                                                                   <  2:3290405/93‑1 (MQ=255)

                      gtgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGt                                                                   <  1:3019427/93‑1 (MQ=255)

                      gtgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGt                                                                   <  1:5395713/93‑1 (MQ=255)

                      atgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGt                                                                   <  1:2824651/92‑1 (MQ=255)

                         tgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTaa                                                                >  2:642333/1‑93 (MQ=255)

                         tgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTaa                                                                <  2:2021774/93‑1 (MQ=255)

                          gtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAAt                                                               >  1:3163452/1‑93 (MQ=255)

                             cAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                            <  2:5535725/93‑1 (MQ=255)

                             cAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                            <  1:4909468/93‑1 (MQ=255)

                              agagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCg                                                           <  1:2583520/93‑1 (MQ=255)

                               gagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCgg                                                          <  2:4807989/93‑1 (MQ=255)

                               gagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCgg                                                          >  2:3614060/1‑93 (MQ=255)

                                agCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGgc                                                         >  1:3326404/1‑93 (MQ=255)

                                   ggCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGca                                                      >  2:5873205/1‑93 (MQ=255)

                                    gCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcac                                                     >  2:1346860/1‑93 (MQ=255)

                                     cATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcaca                                                    >  1:6691818/1‑93 (MQ=255)

                                      aacAAAGTGCTGGGACGGTCAGCTCTCGATAGCACCGCAACACCTTGCTATGTCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaa                                                   <  1:4732913/91‑1 (MQ=255)

                                      aTCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaa                                                   >  2:6047601/1‑93 (MQ=255)

                                      aTCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaa                                                   >  1:2723094/1‑93 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  >  1:4021049/1‑93 (MQ=255)

                                         aaaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACCCCTTGCTATGGCCGCGATACGGCCCCGGCGCTGGTTAATGGCGGCGCACAAATa                                                >  1:1594694/1‑93 (MQ=255)

                                         aaaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATa                                                <  1:3670024/93‑1 (MQ=255)

                                          aaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAg                                               >  1:2297898/1‑93 (MQ=255)

                                          aaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAg                                               >  2:4079244/1‑93 (MQ=255)

                                          aaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAg                                               <  2:2542110/93‑1 (MQ=255)

                                           aGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAgc                                              <  1:314106/93‑1 (MQ=255)

                                              gCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAgccgc                                           <  1:1613866/93‑1 (MQ=255)

                                                tGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCa                                         >  2:2363808/1‑93 (MQ=255)

                                                tGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCa                                         >  2:2934516/1‑93 (MQ=255)

                                                tGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCa                                         >  1:2865871/1‑93 (MQ=255)

                                                tGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCa                                         >  2:6225621/1‑93 (MQ=255)

                                                 gggACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCAt                                        <  2:4621180/93‑1 (MQ=255)

                                                   gACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                                      <  2:3192905/93‑1 (MQ=255)

                                                   gACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                                      >  2:4484015/1‑93 (MQ=255)

                                                   gACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                                      <  2:4525216/93‑1 (MQ=255)

                                                   gACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                                      <  2:67469/93‑1 (MQ=255)

                                                    gcTGGCACGGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGAc                                     <  2:2684525/92‑1 (MQ=255)

                                                      tGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGc                                   >  2:150652/1‑93 (MQ=255)

                                                      tGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGc                                   <  1:6115165/93‑1 (MQ=255)

                                                       ggCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGcc                                  <  1:4521598/93‑1 (MQ=255)

                                                         cAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCgg                                <  1:6498368/93‑1 (MQ=255)

                                                          aGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCggg                               <  1:1891695/93‑1 (MQ=255)

                                                           gcgcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGc                              >  2:70189/1‑93 (MQ=255)

                                                           gcgcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGc                              >  2:1592880/1‑93 (MQ=255)

                                                           gcgcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGc                              >  2:56544/1‑93 (MQ=255)

                                                            cgcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCg                             >  1:222533/1‑93 (MQ=255)

                                                              cTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAg                           >  2:432906/1‑93 (MQ=255)

                                                                cGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCt                         <  1:2452668/93‑1 (MQ=255)

                                                                 gATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTg                        >  1:4092081/1‑93 (MQ=255)

                                                                   tAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGc                      >  1:2173855/1‑93 (MQ=255)

                                                                   tAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGc                      >  1:2023884/1‑93 (MQ=255)

                                                                    aGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGca                     <  2:6638912/93‑1 (MQ=255)

                                                                     gCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcat                    >  1:1130884/1‑93 (MQ=255)

                                                                      cACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatc                   <  1:3282876/93‑1 (MQ=255)

                                                                      cACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatc                   <  1:6430042/93‑1 (MQ=255)

                                                                       cccGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatca                  <  1:5064296/92‑1 (MQ=255)

                                                                       cccGCAACACCTTGCGATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatca                  <  1:4355051/92‑1 (MQ=255)

                                                                       aCCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatca                  <  2:1294936/93‑1 (MQ=255)

                                                                       aCCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatca                  >  1:6235607/1‑93 (MQ=255)

                                                                         cGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATg                >  1:2200259/1‑93 (MQ=255)

                                                                         cGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATg                <  2:5795208/93‑1 (MQ=255)

                                                                          gCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCCTCATgg               >  1:1494161/1‑93 (MQ=255)

                                                                           cAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGc              <  2:2780474/93‑1 (MQ=255)

                                                                           cAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGc              <  1:4581647/93‑1 (MQ=255)

                                                                            aaCACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCt             <  2:6510197/93‑1 (MQ=255)

                                                                            aaCACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCt             >  2:6179326/1‑93 (MQ=255)

                                                                               acCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtttt          >  2:3437067/1‑91 (MQ=255)

                                                                               acCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtttt          >  1:6088064/1‑91 (MQ=255)

                                                                               acCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGc          <  2:1596110/93‑1 (MQ=255)

                                                                                ccTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCttttg         >  1:2523698/1‑90 (MQ=255)

                                                                                 cTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCttttgg        >  1:5439506/1‑89 (MQ=255)

                                                                                  ttGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGga       <  2:6464913/93‑1 (MQ=255)

                                                                                     cTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGaa    <  2:3637154/93‑1 (MQ=255)

                                                                                       aTGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATc  <  2:1093048/93‑1 (MQ=255)

                                                                                       aTGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATc  <  2:2431592/93‑1 (MQ=255)

TACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATC  >  NC_000913/2206933‑2207112

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