BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,949,359	C→T	I179I (ATC→ATT)	yecD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,359	0	C	T	76.5%	1926.0 / inf	903	I179I (ATC→ATT)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (106/106); new base T (391/299); total (498/405)
Fisher's exact test for biased strand distribution p-value = 9.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

Only 100 of 924 total aligned reads displayed.

ACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT  >  NC_000913/1949267‑1949451

   gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAc                                                                                           <  2:4615562/93‑1 (MQ=255)

   gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaat                                                                                           >  1:2466985/1‑92 (MQ=255)

      cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                                        >  2:5893710/1‑93 (MQ=255)

        cAATGCCTGGGAGCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                      <  1:2367178/93‑1 (MQ=255)

        cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                      >  2:1546695/1‑93 (MQ=255)

         aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                                     <  1:4094839/93‑1 (MQ=255)

         aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                                     <  2:1317214/93‑1 (MQ=255)

         aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                                     >  2:3636276/1‑93 (MQ=255)

           tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                   <  1:5282379/93‑1 (MQ=255)

              cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                >  1:4798124/1‑93 (MQ=255)

              cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCCCCAGGGCCGAGCAACACAATAACAGCATTATTCATTTCTAccc                                                                                >  1:4977802/1‑93 (MQ=255)

               tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                               >  1:2011535/1‑93 (MQ=255)

               tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                               >  1:2240427/1‑93 (MQ=255)

               tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                               >  2:3254120/1‑93 (MQ=255)

                gggAACTTGGCTTTAATCTGGTGATTCCCGAAGACGCCTGTAGTCCGGCCAGCGCCGAGCAGCACAATAAAAGTATCAACCATATCTACCcgc                                                                              >  2:2625724/1‑93 (MQ=255)

                  gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                                            >  2:2762547/1‑93 (MQ=255)

                  gAACTTGGATTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                                            >  1:262596/1‑93 (MQ=255)

                   aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCa                                                                           >  1:962932/1‑93 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCAt                                                                          >  1:5238051/1‑93 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                                          <  1:5793037/93‑1 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                                          <  2:919030/93‑1 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                                          <  1:3447450/93‑1 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATAAATCATATCTACCCGCGCAt                                                                          >  2:4246997/1‑93 (MQ=255)

                     cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                         <  1:1921657/93‑1 (MQ=255)

                     cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                         <  2:2914317/93‑1 (MQ=255)

                     cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATc                                                                         >  2:4827982/1‑93 (MQ=255)

                     cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATc                                                                         >  1:1812247/1‑93 (MQ=255)

                      ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                                        >  1:5163445/1‑93 (MQ=255)

                        ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                      <  1:1095881/93‑1 (MQ=255)

                         gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                                     >  1:4691159/1‑93 (MQ=255)

                         gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                                     <  1:5699021/93‑1 (MQ=255)

                          cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                                    >  2:4917829/1‑93 (MQ=255)

                          cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                                    >  2:349493/1‑93 (MQ=255)

                            ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                  <  2:5643397/93‑1 (MQ=255)

                            ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                  >  2:4029363/1‑93 (MQ=255)

                            ttAATCTGGTCATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                  <  2:60218/93‑1 (MQ=255)

                                tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCg                                                              >  2:199461/1‑93 (MQ=255)

                                 cTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGt                                                             >  1:768675/1‑93 (MQ=255)

                                   ggTGATGGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                           <  1:5427524/93‑1 (MQ=255)

                                    gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                          <  1:1378345/93‑1 (MQ=255)

                                     tGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCg                                                         >  1:3124618/1‑93 (MQ=255)

                                      gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                        >  2:1691020/1‑93 (MQ=255)

                                       aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                       >  2:2122377/1‑93 (MQ=255)

                                           ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                   <  2:3778476/93‑1 (MQ=255)

                                            cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                                  <  1:3652300/93‑1 (MQ=255)

                                             gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                                 <  2:2865646/93‑1 (MQ=255)

                                             gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                                 >  2:1826471/1‑93 (MQ=255)

                                             gAAGACGCCTGTAGGGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                                 <  1:2845356/93‑1 (MQ=255)

                                              aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                >  1:3990288/1‑93 (MQ=255)

                                               aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAt                                               <  1:3547650/93‑1 (MQ=255)

                                               aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAt                                               <  1:5880867/93‑1 (MQ=255)

                                                 aCGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                             <  1:655820/93‑1 (MQ=255)

                                                 aCGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                             >  2:2293281/1‑93 (MQ=255)

                                                     cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttca                                         <  2:4291887/93‑4 (MQ=255)

                                                     cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttca                                         >  2:5535945/1‑90 (MQ=255)

                                                     cTGTAGTGCCGCCAGCGCAGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                         >  1:1476325/1‑93 (MQ=255)

                                                      tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                        >  2:1394770/1‑93 (MQ=255)

                                                      tGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcac                                        >  2:4514268/1‑89 (MQ=255)

                                                          gTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                     <  2:5310994/76‑1 (MQ=255)

                                                           tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                   >  1:1210387/1‑93 (MQ=255)

                                                            gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTa                                  <  1:2546484/93‑1 (MQ=255)

                                                            gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                  >  1:4139230/1‑93 (MQ=255)

                                                              cgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATg                                <  2:1457204/93‑1 (MQ=255)

                                                              cgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATg                                <  2:175376/93‑1 (MQ=255)

                                                               gccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGa                               <  1:3731880/93‑1 (MQ=255)

                                                               gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                               >  1:1637459/1‑93 (MQ=255)

                                                                ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                              >  2:768405/1‑93 (MQ=255)

                                                                 cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt                             >  2:823686/1‑93 (MQ=255)

                                                                  aGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAttt                            <  1:5591134/93‑1 (MQ=255)

                                                                    cgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTAc                          >  1:2741122/1‑93 (MQ=255)

                                                                     gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                         >  2:4687889/1‑93 (MQ=255)

                                                                      ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                        <  1:5361774/93‑1 (MQ=255)

                                                                      ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                        <  2:1219882/93‑1 (MQ=255)

                                                                      ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt                        <  2:5856509/93‑1 (MQ=255)

                                                                      ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt                        >  1:2114462/1‑93 (MQ=255)

                                                                       cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATc                       >  2:3306977/1‑93 (MQ=255)

                                                                       cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATc                       >  1:508779/1‑93 (MQ=255)

                                                                        gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                      >  1:4523542/1‑93 (MQ=255)

                                                                           cagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTc                   <  2:5443688/93‑1 (MQ=255)

                                                                           caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                   <  1:1001413/90‑1 (MQ=255)

                                                                           caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                   <  1:4410911/90‑1 (MQ=255)

                                                                           caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                   <  1:4258947/90‑1 (MQ=255)

                                                                           caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                   <  2:51628/90‑1 (MQ=255)

                                                                            aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                  <  2:64237/91‑1 (MQ=255)

                                                                            aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                  >  1:6084918/3‑93 (MQ=255)

                                                                            aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                  >  2:507727/3‑93 (MQ=255)

                                                                              cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc                <  2:3902463/93‑3 (MQ=255)

                                                                              cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc                <  2:4198326/93‑3 (MQ=255)

                                                                              cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc                <  2:104976/93‑3 (MQ=255)

                                                                                caATAACAGTATCAACCATATCTACCCGCGCCTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg              >  1:2995530/1‑93 (MQ=255)

                                                                                caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCg              >  1:5979484/1‑93 (MQ=255)

                                                                                caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGCGTGCGGAGCGTGGAAGAGATTCTTCACGCGTTATGATATACATCGGGCTGCCg              >  1:5468167/1‑93 (MQ=255)

                                                                                  aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa            <  2:2529777/93‑1 (MQ=255)

                                                                                   tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa           <  2:5720040/93‑1 (MQ=255)

                                                                                     aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATg         <  2:4414209/93‑1 (MQ=255)

                                                                                     aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATg         <  2:1141047/93‑1 (MQ=255)

                                                                                      cAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg        <  1:3298978/93‑1 (MQ=255)

                                                                                      cAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg        <  2:313782/93‑1 (MQ=255)

                                                                                        gTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc      >  1:4916942/1‑93 (MQ=255)

                                                                                            cAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAt  >  1:3360244/1‑93 (MQ=255)

ACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT  >  NC_000913/1949267‑1949451

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