BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	978,817	A→G	D1090G (GAC→GGC)	mukB →	MukB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	978,817	0	A	G	88.3%	584.9 / 58.3	230	D1090G (GAC→GGC)	mukB	MukB
Reads supporting (aligned to +/- strand): ref base A (20/7); new base G (118/85); total (138/92)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 232 total aligned reads displayed.

TCGCCGTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGGTGTGCGGT  >  NC_000913/978725‑978907

tCGCCGTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGa                                                                                            <  2:2597740/93‑1 (MQ=255)

  gCCGTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCa                                                                                          >  2:1676531/1‑93 (MQ=255)

  gCCGTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACa                                                                                          >  1:6722168/1‑93 (MQ=255)

     gTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAAcc                                                                                       >  1:6234324/1‑93 (MQ=255)

     gTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAAcc                                                                                       <  1:4362827/93‑1 (MQ=255)

     gTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAAcc                                                                                       >  2:4400556/1‑93 (MQ=255)

      tGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCt                                                                                      <  2:1389038/93‑1 (MQ=255)

          gAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGAcc                                                                                  >  2:3889879/1‑93 (MQ=255)

           aGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGAccc                                                                                 >  1:1484279/1‑93 (MQ=255)

              tGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCa                                                                              >  1:1381623/1‑93 (MQ=255)

              tGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCa                                                                              >  2:5352413/1‑93 (MQ=255)

                cATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCaaa                                                                            >  2:1227037/1‑93 (MQ=255)

                cATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCaaa                                                                            >  1:6811846/1‑93 (MQ=255)

                cATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCaaa                                                                            >  1:3616762/1‑93 (MQ=255)

                 aTGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAAc                                                                           <  1:2486052/93‑1 (MQ=255)

                 aTGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAAc                                                                           <  1:6955491/93‑1 (MQ=255)

                 aTGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAAc                                                                           <  1:4140646/93‑1 (MQ=255)

                  tGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACt                                                                          >  1:2581610/1‑93 (MQ=255)

                    cgcAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTgc                                                                        >  2:4446862/1‑93 (MQ=255)

                     gcAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTgcg                                                                       >  1:4942443/1‑93 (MQ=255)

                        aCTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCaa                                                                    <  1:7163903/93‑1 (MQ=255)

                        aCTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCaa                                                                    <  2:4339891/93‑1 (MQ=255)

                        aCTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCaa                                                                    >  1:4612818/1‑93 (MQ=255)

                         cTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAg                                                                   >  1:733475/1‑93 (MQ=255)

                         cTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAg                                                                   >  2:2991186/1‑93 (MQ=255)

                          tGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGc                                                                  <  1:692634/93‑1 (MQ=255)

                           gAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCt                                                                 <  1:29703/93‑1 (MQ=255)

                            aGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTg                                                                >  1:5073761/1‑93 (MQ=255)

                            aGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTg                                                                >  1:829853/1‑93 (MQ=255)

                             gCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTgg                                                               <  1:1841517/93‑1 (MQ=255)

                              cAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGa                                                              <  1:3192755/93‑1 (MQ=255)

                              cAATAACCGTTCACGCCGCAATCAACCGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGa                                                              >  2:2771571/1‑93 (MQ=255)

                                aTAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGc                                                            >  2:2411887/1‑93 (MQ=255)

                                 tAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCg                                                           >  2:3173616/1‑93 (MQ=255)

                                 tAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCg                                                           >  2:6190199/1‑93 (MQ=255)

                                 tAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCg                                                           >  2:3464246/1‑93 (MQ=255)

                                 tAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCg                                                           <  2:5761637/93‑1 (MQ=255)

                                 tAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCg                                                           <  1:2120641/93‑1 (MQ=255)

                                  aaCCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCgg                                                          <  1:1788833/93‑1 (MQ=255)

                                   aCCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCggg                                                         >  2:6243401/1‑93 (MQ=255)

                                      gTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGTTGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGAtt                                                      >  1:5302375/1‑93 (MQ=255)

                                      gTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGAtt                                                      >  1:6349839/1‑93 (MQ=255)

                                       ttCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTa                                                     >  2:5888623/1‑93 (MQ=255)

                                        tCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTAc                                                    >  2:1937279/1‑93 (MQ=255)

                                         cACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACt                                                   <  1:5800298/93‑1 (MQ=255)

                                          aCGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACtt                                                  >  1:2987690/1‑93 (MQ=255)

                                          aCGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACtt                                                  >  2:6937521/1‑93 (MQ=255)

                                            gccgcAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTg                                                >  2:2124964/1‑93 (MQ=255)

                                                cAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGAt                                            <  2:2816972/93‑1 (MQ=255)

                                                cAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGAt                                            >  1:945907/1‑93 (MQ=255)

                                                 aaTCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATg                                           <  2:4282264/93‑1 (MQ=255)

                                                   tCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATgcg                                         >  1:5954380/1‑93 (MQ=255)

                                                   tCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATgcg                                         >  1:5385200/1‑93 (MQ=255)

                                                   tCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATgcg                                         >  1:6503494/1‑93 (MQ=255)

                                                    cAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATgcgc                                        <  2:6585599/93‑1 (MQ=255)

                                                    cAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATgcgc                                        <  2:2887047/93‑1 (MQ=255)

                                                     aaCTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGcgcg                                       >  1:1183101/1‑93 (MQ=255)

                                                     aaCTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGcgcg                                       >  2:742369/1‑93 (MQ=255)

                                                     aaCTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGcgcg                                       >  1:7408481/1‑93 (MQ=255)

                                                        tGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGc                                    >  1:5045300/1‑93 (MQ=255)

                                                         ggAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCa                                   <  2:3347060/93‑1 (MQ=255)

                                                         ggAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCa                                   >  1:1292838/1‑93 (MQ=255)

                                                         ggAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCa                                   <  2:6880847/93‑1 (MQ=255)

                                                           aaaaaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAgg                                 >  2:888918/1‑93 (MQ=255)

                                                            aaaaGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGt                                >  1:1038288/1‑93 (MQ=255)

                                                             aaaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTgaga                                             >  1:5186957/1‑79 (MQ=255)

                                                             aaaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTa                               >  1:3917922/1‑93 (MQ=255)

                                                             aaaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTa                               <  1:2900445/93‑1 (MQ=255)

                                                             aaaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCCGGATTACTTTgaga                                             <  2:5186957/79‑1 (MQ=255)

                                                              aaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAg                              <  2:4955461/93‑1 (MQ=255)

                                                              aaGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAg                              >  2:2391709/1‑93 (MQ=255)

                                                               aGCGCTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGt                             <  1:4406067/93‑1 (MQ=255)

                                                                gcgcTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTg                            >  2:7272183/1‑93 (MQ=255)

                                                                 cgcTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGa                           <  2:6833335/93‑1 (MQ=255)

                                                                  gcTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGAc                          >  2:2446007/1‑93 (MQ=255)

                                                                   cTTACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGAcc                         >  2:1054895/1‑93 (MQ=255)

                                                                    ttACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCg                        >  1:2727228/1‑93 (MQ=255)

                                                                    ttACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCg                        >  1:3129527/1‑93 (MQ=255)

                                                                    ttACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCg                        >  2:1370107/1‑93 (MQ=255)

                                                                    ttACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCg                        >  2:367704/1‑93 (MQ=255)

                                                                     tACCTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGc                       >  1:263145/1‑93 (MQ=255)

                                                                        cTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCaa                    <  2:1566824/93‑1 (MQ=255)

                                                                        cTTCTGCGAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCaa                    >  1:3366856/1‑93 (MQ=255)

                                                                        cTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCaa                    >  2:5966787/1‑93 (MQ=255)

                                                                         ttCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCaaa                   >  2:5626814/1‑93 (MQ=255)

                                                                          tCTGCGAAGCGGAGATGGGCAACCGGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAg                  <  1:5614927/93‑1 (MQ=255)

                                                                           cTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGc                 <  1:6594463/93‑1 (MQ=255)

                                                                               gAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGc             >  2:3266634/1‑93 (MQ=255)

                                                                               gAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGc             <  1:2695646/93‑1 (MQ=255)

                                                                               gAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGc             >  2:6478915/1‑93 (MQ=255)

                                                                               gAAGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGc             <  2:4761234/93‑1 (MQ=255)

                                                                                aaGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCt            >  1:5555134/1‑93 (MQ=255)

                                                                                aaGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCt            >  1:6121688/1‑93 (MQ=255)

                                                                                 aGCGGAGATGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTg           >  2:4611419/1‑93 (MQ=255)

                                                                                    ggAGATGGGCACCCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGgtg        <  2:2519990/93‑1 (MQ=255)

                                                                                     gagaTGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGgtgt       <  1:3742586/93‑1 (MQ=255)

                                                                                     gagaTGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGgtgt       >  2:7152628/1‑93 (MQ=255)

                                                                                       gaTGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGGTGTGc     <  2:5822451/93‑1 (MQ=255)

                                                                                         tGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGGTGTGCgg   <  1:6102170/93‑1 (MQ=255)

                                                                                         tGGGCAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGGTGTGCgg   <  1:998412/93‑1 (MQ=255)

TCGCCGTGACGAGCTGCATGCGCAACTGAGCAATAACCGTTCACGCCGCAATCAACTGGAAAAAGCGCTTACCTTCTGCGAAGCGGAGATGGACAACCTGACCCGCAAACTGCGCAAGCTGGAGCGGGATTACTTTGAGATGCGCGAGCAGGTAGTGACCGCCAAAGCGGGCTGGTGTGCGGT  >  NC_000913/978725‑978907

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