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.1%	348.0 / ‑4.7	110	D346G (GAC→GGC)	yehP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (53/56); total (54/56)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 114 total aligned reads displayed.

TACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATC  >  NC_000913/2206933‑2207112

tACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGcc                                                                                         >  1:3825052/1‑93 (MQ=255)

tACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGcc                                                                                         >  2:3582349/1‑93 (MQ=255)

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

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

   tGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCcgcg                                                                                      >  1:3981863/1‑93 (MQ=255)

   tGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCcgcg                                                                                      >  1:4408699/1‑93 (MQ=255)

   tGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGAGCTCGATAGCACCGCAACACCTTGCTATGGCcgcg                                                                                      >  1:3981864/1‑93 (MQ=255)

    gACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGa                                                                                     <  1:2127245/93‑1 (MQ=255)

    gACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGa                                                                                     >  2:1731748/1‑93 (MQ=255)

      cGCATCAGGTGAAAAAGTGTGTCCAGAGCGGGATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATa                                                                                   <  1:3537732/93‑1 (MQ=255)

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

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACTGCAACACCTTGCTATGGCCGCGATAc                                                                                  >  1:4323607/1‑93 (MQ=255)

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATAc                                                                                  >  2:3206037/1‑93 (MQ=255)

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

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

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGAGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATAc                                                                                  <  1:1334907/93‑1 (MQ=255)

        cATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACg                                                                                 >  1:190002/1‑93 (MQ=255)

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

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

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

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

             ggTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCa                                                                            <  1:1261688/93‑1 (MQ=255)

              gTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAg                                                                           <  1:2216670/93‑1 (MQ=255)

              gTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAg                                                                           <  2:391730/93‑1 (MQ=255)

                 aaaaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgcg                                                                        >  1:387992/1‑93 (MQ=255)

                    aaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTg                                                                     <  2:4423565/93‑1 (MQ=255)

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

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

                     aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTgg                                                                    >  2:2230427/1‑93 (MQ=255)

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

                     aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTgg                                                                    >  1:3164742/1‑93 (MQ=255)

                      gtgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGAGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGt                                                                   >  1:2221577/1‑93 (MQ=255)

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

                           tCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATg                                                              >  1:1324717/1‑93 (MQ=255)

                             cAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                            >  1:12942/1‑93 (MQ=255)

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

                              agagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCg                                                           >  2:19747/1‑93 (MQ=255)

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

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

                               gagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCgg                                                          >  1:4018259/1‑93 (MQ=255)

                                agCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCATCCCCTTGCTATGGCCGCGATACGGCCCAGGCTCTGGTTAATGTCGgc                                                         <  2:3324430/93‑1 (MQ=255)

                                 gCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGgcg                                                        >  2:4368553/1‑93 (MQ=255)

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

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

                                     cATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcaca                                                    <  2:2340478/93‑1 (MQ=255)

                                      aTCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTGCGCGCACaa                                                   >  2:3462301/1‑93 (MQ=255)

                                        cAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAAt                                                 <  1:3453486/93‑1 (MQ=255)

                                         aaaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAAtc                                                <  2:2881275/93‑2 (MQ=255)

                                         aaaGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATa                                                >  1:3220227/1‑93 (MQ=255)

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

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

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

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

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

                                                tGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCa                                         <  1:4434631/93‑1 (MQ=255)

                                                 gggACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCAt                                        >  2:2479947/1‑93 (MQ=255)

                                                  ggACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATg                                       <  2:4542368/93‑1 (MQ=255)

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

                                                    aCTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGAc                                     >  2:2574551/1‑93 (MQ=255)

                                                     cTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGAcc                                    >  2:1070674/1‑92 (MQ=255)

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

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

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

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

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

                                                       ggCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGcc                                  >  2:3336900/1‑93 (MQ=255)

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

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

                                                           gcgcTCGATAGCACCGCAACCCCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGc                              <  2:2288430/93‑1 (MQ=255)

                                                             gcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGa                            <  2:2214137/93‑1 (MQ=255)

                                                             gcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGa                            <  1:2084754/93‑1 (MQ=255)

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

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

                                                              cTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAg                           <  2:3267233/93‑1 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          >  1:364159/1‑93 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          >  1:3318604/1‑93 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          <  1:3275059/93‑1 (MQ=255)

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

                                                                 gATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTg                        <  2:4535901/93‑1 (MQ=255)

                                                                      cACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatc                   >  2:3098723/1‑93 (MQ=255)

                                                                      cACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatc                   >  2:4187180/1‑93 (MQ=255)

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

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

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

                                                                        ccGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatcat                 <  2:2747779/93‑1 (MQ=255)

                                                                        ccGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatcat                 <  1:4375504/93‑1 (MQ=255)

                                                                        ccGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatcat                 <  1:4503394/93‑1 (MQ=255)

                                                                        ccGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatcat                 <  1:2954237/93‑1 (MQ=255)

                                                                          gCAACACCTTCCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATgg               <  1:2821987/93‑1 (MQ=255)

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

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

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

                                                                             acacCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtt            >  1:1501696/1‑93 (MQ=255)

                                                                             acacCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtt            >  1:473887/1‑93 (MQ=255)

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

                                                                                ccTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTCg         >  1:906885/1‑93 (MQ=255)

                                                                                 cTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCgg        <  1:2965312/93‑1 (MQ=255)

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

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

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

TACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATC  >  NC_000913/2206933‑2207112

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