BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,283,731	A→G	N43S (AAC→AGC)	narH →	nitrate reductase A subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,283,731	0	A	G	100.0%	1276.9 / NA	306	N43S (AAC→AGC)	narH	nitrate reductase A subunit beta
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (154/152); total (154/152)

Only 100 of 317 total aligned reads displayed.

GCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCAAACTGCAG  >  NC_000913/1283633‑1283828

    aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGa                                                                                             <  2:2358097/101‑1 (MQ=255)

    aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGa                                                                                             >  2:2358802/1‑101 (MQ=255)

    aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGa                                                                                             >  1:2358801/1‑101 (MQ=255)

          gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCaa                                                                                       >  2:2358807/1‑101 (MQ=255)

          gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCaa                                                                                       <  1:2358521/101‑1 (MQ=255)

          gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCaa                                                                                       <  1:2358395/101‑1 (MQ=255)

           aTAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAg                                                                                      >  1:2358809/1‑101 (MQ=255)

            tAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGc                                                                                     <  2:2358241/101‑1 (MQ=255)

            tAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGc                                                                                     >  2:2358812/1‑101 (MQ=255)

            tAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGc                                                                                     <  1:2358329/101‑1 (MQ=255)

              aGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCg                                                                                   >  2:2358816/1‑101 (MQ=255)

                 gCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGc                                                                                >  2:2358820/1‑101 (MQ=255)

                  cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATACGTGGTTCAACAGCGTGGAAACCAAGCCGGGcc                                                                               <  1:2358551/101‑1 (MQ=255)

                   aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCa                                                                              >  1:2358824/1‑101 (MQ=255)

                    tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAg                                                                             >  1:2358828/1‑101 (MQ=255)

                    tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAg                                                                             >  2:2358830/1‑101 (MQ=255)

                     cGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAgg                                                                            >  1:2358832/1‑101 (MQ=255)

                      ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCggg                                                                                 <  1:2358351/95‑1 (MQ=255)

                      ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAggg                                                                           <  1:2358365/101‑1 (MQ=255)

                        cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCt                                                                         >  2:2358836/1‑101 (MQ=255)

                          gCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCg                                                                                                 <  2:2358363/75‑1 (MQ=255)

                          gCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAggg                                                                           >  1:2358840/1‑97 (MQ=255)

                          gCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTc                                                                       <  1:2358470/101‑1 (MQ=255)

                           ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTcc                                                                      <  1:2358150/101‑1 (MQ=255)

                            cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTccc                                                                     >  2:2358842/1‑101 (MQ=255)

                            cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTccc                                                                     >  2:2358845/1‑101 (MQ=255)

                             acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAg                                                                                      <  1:2358316/83‑1 (MQ=255)

                             acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGcc                                                                                    <  2:2358473/85‑1 (MQ=255)

                             acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCg                                                                    >  2:2358848/1‑101 (MQ=255)

                               acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCgac                                                                  >  1:2358852/1‑101 (MQ=255)

                                ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCgac                                                                  <  1:2358491/100‑1 (MQ=255)

                                ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGAtt                                                                        <  1:2358186/94‑1 (MQ=255)

                                  tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactga                                                               >  2:2358854/1‑101 (MQ=255)

                                   gTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                                                              <  1:2358207/101‑1 (MQ=255)

                                    ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact                                                             <  1:2358747/101‑1 (MQ=255)

                                     tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGactgactg                                                            >  1:2358858/1‑101 (MQ=255)

                                      cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTgg                                                           <  2:2358349/101‑1 (MQ=255)

                                       aGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTggg                                                          >  2:2358862/1‑101 (MQ=255)

                                         ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaa                                                        <  1:2358421/101‑1 (MQ=255)

                                         ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaa                                                        <  1:2358450/101‑1 (MQ=255)

                                          tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaa                                                       <  1:2358387/101‑1 (MQ=255)

                                          tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaa                                                       >  1:2358866/1‑101 (MQ=255)

                                           aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaaa                                                      >  2:2358869/1‑101 (MQ=255)

                                            ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAAc                                                     >  1:2358871/1‑101 (MQ=255)

                                             cTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAg                                                                             <  1:2358374/76‑1 (MQ=255)

                                              tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCa                                                   <  1:2358110/101‑1 (MQ=255)

                                               gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAg                                                  >  1:2358874/1‑101 (MQ=255)

                                                 aaaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGa                                                >  1:2358878/1‑101 (MQ=255)

                                                   aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaa                                              >  1:2358882/1‑101 (MQ=255)

                                                    aaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaa                                             >  2:2358884/1‑101 (MQ=255)

                                                    aaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaa                                             >  1:2358887/1‑101 (MQ=255)

                                                     aCGTCTGGCCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGcc                                                                               <  2:2358234/66‑1 (MQ=255)

                                                     aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaaa                                            >  2:2358890/1‑101 (MQ=255)

                                                      cGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactg                                                                <  2:2358386/80‑1 (MQ=255)

                                                       gTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATa                                          >  1:2358894/1‑101 (MQ=255)

                                                       gTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATa                                          >  2:2358898/1‑101 (MQ=255)

                                                        tCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATAc                                         <  1:2358710/101‑1 (MQ=255)

                                                         cTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACa                                        >  1:2358900/1‑101 (MQ=255)

                                                          tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACaa                                       >  2:2358904/1‑101 (MQ=255)

                                                           ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaa                                                       <  2:2358317/84‑1 (MQ=255)

                                                           ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACaaa                                      >  1:2358905/1‑101 (MQ=255)

                                                           ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACaaa                                      <  2:2358431/101‑1 (MQ=255)

                                                             aCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAgg                                    >  1:2358908/1‑101 (MQ=255)

                                                              ccAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACaa                                       <  2:2358455/97‑1 (MQ=255)

                                                              ccAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAggc                                   <  1:2358261/101‑1 (MQ=255)

                                                                aGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAggcgg                                 >  2:2358913/1‑101 (MQ=255)

                                                                aGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAg                                     >  1:2358914/1‑97 (MQ=255)

                                                                 gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAggcggc                                <  2:2358649/101‑1 (MQ=255)

                                                                   cGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAAt                                           <  2:2358292/88‑1 (MQ=255)

                                                                   cGTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTg                              <  1:2358310/101‑1 (MQ=255)

                                                                    gTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGa                                                         <  1:2358411/73‑1 (MQ=255)

                                                                    gTGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTgg                             >  2:2358919/1‑101 (MQ=255)

                                                                     tGAAGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGa                            >  1:2358920/1‑101 (MQ=255)

                                                                       aaGGCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATc                          <  1:2358239/101‑1 (MQ=255)

                                                                         ggCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCg                        <  1:2358565/101‑1 (MQ=255)

                                                                         ggCGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCg                        >  1:2358926/1‑101 (MQ=255)

                                                                           cGTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAg                                     >  2:2358929/1‑86 (MQ=255)

                                                                            gTGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTaa                     <  2:2358481/101‑1 (MQ=255)

                                                                             tGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTaaa                    <  2:2358791/101‑1 (MQ=255)

                                                                             tGGAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTaaa                    >  2:2358933/1‑101 (MQ=255)

                                                                              ggAATATGCGTAGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACaa                                       <  2:2358379/81‑1 (MQ=255)

                                                                               gAATATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTaa                     <  2:2358405/98‑1 (MQ=255)

                                                                                aaTATGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATc                 >  1:2358936/1‑101 (MQ=255)

                                                                                 atatGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCa                <  2:2358550/101‑1 (MQ=255)

                                                                                 atatGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCa                <  1:2358604/101‑1 (MQ=255)

                                                                                 atatGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCa                <  2:2358524/101‑1 (MQ=255)

                                                                                 atatGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCa                >  2:2358939/1‑101 (MQ=255)

                                                                                 atatGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCa                <  1:2358289/101‑1 (MQ=255)

                                                                                  tatGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCaa               <  2:2358391/101‑1 (MQ=255)

                                                                                   atGCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaa                                                       >  1:2358941/1‑60 (MQ=255)

                                                                                     gCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACgg            <  1:2358631/101‑1 (MQ=255)

                                                                                     gCGTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACgg            <  1:2358501/101‑1 (MQ=255)

                                                                                       gTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCa          <  2:2358254/101‑1 (MQ=255)

                                                                                       gTGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCa          >  2:2358945/1‑101 (MQ=255)

                                                                                        tGGTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCaa         >  1:2358946/1‑101 (MQ=255)

                                                                                         ggTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCaaa        <  2:2358534/101‑1 (MQ=255)

                                                                                          gTTCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCAAAc       >  2:2358949/1‑101 (MQ=255)

                                                                                            tCAACAGCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCAAACTg     <  2:2358438/101‑1 (MQ=255)

                                                                                             caacagCGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCAAACTGc    >  1:2358955/1‑101 (MQ=255)

                                                                                               acagcGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCg                        <  2:2358265/79‑1 (MQ=255)

GCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAACGGCAAACTGCAG  >  NC_000913/1283633‑1283828

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