BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,403,088	T→C	S108S (AGT→AGC)	ydaK →	putative transcriptional regulator LYSR‑type

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,403,088	0	T	C	100.0%	624.2 / NA	198	S108S (AGT→AGC)	ydaK	putative transcriptional regulator LYSR‑type
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (87/111); total (87/111)

Only 100 of 204 total aligned reads displayed.

GCTGCGCGCAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGTCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCGATGATTAA  >  NC_000913/1402998‑1403177

gCTGCGCGCAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCt                                                                                         >  1:6216172/1‑93 (MQ=255)

gCTGCGCGCAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACGGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCt                                                                                         <  1:1096121/93‑1 (MQ=255)

  tGCGCGCAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGa                                                                                       >  2:1334476/1‑93 (MQ=255)

     gcgcAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGc                                                                                    >  1:7603431/1‑93 (MQ=255)

      cgcAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGcc                                                                                   <  2:4147514/93‑1 (MQ=255)

       gcAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCa                                                                                  <  2:3686516/93‑1 (MQ=255)

        cAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAg                                                                                 <  1:7786774/93‑1 (MQ=255)

        cAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAg                                                                                 <  1:4601754/93‑1 (MQ=255)

          gCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCt                                                                               >  2:3584615/1‑93 (MQ=255)

           cccAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTg                                                                              <  2:23837/93‑1 (MQ=255)

            ccAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGt                                                                             >  2:1053606/1‑93 (MQ=255)

            ccAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGt                                                                             >  1:3771530/1‑93 (MQ=255)

             cAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTc                                                                            >  2:3350993/1‑93 (MQ=255)

              aaaagGATATTCCCCAACGACATGAGCAACTGGCCGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCCGCCTGATGCCAGCTGTct                                                                           >  2:512564/4‑92 (MQ=255)

               agagGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCat                                                                          >  2:527818/1‑93 (MQ=255)

                gagGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCata                                                                         >  1:6593349/1‑93 (MQ=255)

                gagGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCata                                                                         <  2:5789227/93‑1 (MQ=255)

                gagGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCata                                                                         >  2:4252591/1‑93 (MQ=255)

                  ggATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATAtc                                                                       <  1:2727441/93‑1 (MQ=255)

                  ggATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATAtc                                                                       >  1:6996956/1‑93 (MQ=255)

                   gATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATAtct                                                                      <  2:7754734/93‑1 (MQ=255)

                    agatTCGCCACCGACAAGGGCACCGGGCAGGGCAGATTAATATCGACATGGGGGCCAGTATTCCCCGCAGGCTGATGCCAGCCGTCATAtctc                                                                     <  2:1536327/91‑1 (MQ=255)

                     tatTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCg                                                                    >  2:5429088/1‑93 (MQ=255)

                      atTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGt                                                                   >  2:511032/1‑93 (MQ=255)

                        tCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGttt                                                                 <  1:340609/93‑1 (MQ=255)

                          gCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTtc                                                               >  1:5015129/1‑93 (MQ=255)

                           ccAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTtca                                                              <  2:4642024/93‑1 (MQ=255)

                           ccAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTtca                                                              <  1:612908/93‑1 (MQ=255)

                            cAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTtcat                                                             >  1:5210258/1‑93 (MQ=255)

                            cAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTtcat                                                             <  2:2221635/93‑1 (MQ=255)

                             aaCGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTtcata                                                            >  1:1735991/1‑92 (MQ=255)

                                gACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATcagc                                                         >  1:4118215/1‑93 (MQ=255)

                                gACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATcagc                                                         <  1:4979450/93‑1 (MQ=255)

                                 aCAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATcagca                                                        >  2:3433195/1‑93 (MQ=255)

                                 aCAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATcagca                                                        <  2:4547375/93‑1 (MQ=255)

                                 aCAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATcagca                                                        <  1:2180425/93‑1 (MQ=255)

                                    aGGGCAACTGGCAGGGCAGGTTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTATCATATCTCGTTTTCATCAgcagca                                                     >  1:5401177/1‑93 (MQ=255)

                                    aGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAgcagca                                                     >  2:1939627/1‑93 (MQ=255)

                                    aGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAgcagca                                                     <  2:6699701/93‑1 (MQ=255)

                                    aGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAgcagca                                                     >  1:7657504/1‑93 (MQ=255)

                                     gggCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCAt                                                    <  1:246333/93‑1 (MQ=255)

                                     gggCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCAt                                                    <  2:7666194/93‑1 (MQ=255)

                                      ggCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATc                                                   <  2:7785447/93‑1 (MQ=255)

                                      ggCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATc                                                   >  2:1728566/1‑93 (MQ=255)

                                      ggCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATc                                                   <  2:7291572/93‑1 (MQ=255)

                                       gCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATcc                                                  <  1:714920/93‑1 (MQ=255)

                                         aaCTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGc                                                >  2:3267568/1‑93 (MQ=255)

                                         aaCTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGc                                                >  1:163779/1‑93 (MQ=255)

                                          aCTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCa                                               <  1:660921/93‑1 (MQ=255)

                                          aCTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCa                                               <  2:1391719/93‑1 (MQ=255)

                                           cTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAg                                              >  2:6840983/1‑93 (MQ=255)

                                            tGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAgg                                             <  2:639740/93‑1 (MQ=255)

                                            tGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAgg                                             <  2:4386799/93‑1 (MQ=255)

                                             ggcagggcagATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGt                                            <  2:6127341/93‑1 (MQ=255)

                                                agggcagATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTaaa                                         >  2:7367542/1‑93 (MQ=255)

                                                 gggcagATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTaaaa                                        >  1:2918869/1‑93 (MQ=255)

                                                 gggcagATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTaaaa                                        <  1:7363856/93‑1 (MQ=255)

                                                  ggcagATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAg                                       <  2:7598875/93‑1 (MQ=255)

                                                   gcagATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGt                                      <  2:6309242/93‑1 (MQ=255)

                                                     agATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGATTTCATCAGCACCATCCGCAGGTAAAAGTAc                                    >  1:1612490/1‑93 (MQ=255)

                                                      gATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACg                                   >  1:3559977/1‑93 (MQ=255)

                                                        ttAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCa                                 >  2:5620266/1‑93 (MQ=255)

                                                        ttAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCa                                 <  1:4291942/93‑1 (MQ=255)

                                                        ttAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCa                                 <  1:6850046/93‑1 (MQ=255)

                                                         tAATATCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCAt                                <  2:7248651/93‑1 (MQ=255)

                                                           atatCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTa                              <  1:5595686/93‑1 (MQ=255)

                                                           atatCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTa                              >  1:2002866/1‑93 (MQ=255)

                                                            tatCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTAt                             <  2:635907/93‑1 (MQ=255)

                                                             atCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATg                            <  2:6934622/93‑1 (MQ=255)

                                                              tCGGCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATgg                           <  1:4855649/93‑1 (MQ=255)

                                                                ggCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGaa                         <  2:2991762/93‑1 (MQ=255)

                                                                ggCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGaa                         <  2:4780520/93‑1 (MQ=255)

                                                                ggCATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGaa                         <  1:1906219/93‑1 (MQ=255)

                                                                  cATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAgg                       >  1:4396301/1‑93 (MQ=255)

                                                                  cATGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAgg                       >  1:921602/1‑93 (MQ=255)

                                                                   aTGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAggg                      >  1:2783996/1‑93 (MQ=255)

                                                                    tGGGGGCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAggga                     >  1:6374756/1‑92 (MQ=255)

                                                                     gggggCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCa                    >  2:3892014/1‑93 (MQ=255)

                                                                      ggggCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCaa                   >  2:7922742/1‑93 (MQ=255)

                                                                      ggggCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCaa                   <  2:5371303/93‑1 (MQ=255)

                                                                      ggggCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCaa                   <  1:1627055/93‑1 (MQ=255)

                                                                       gggCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAAc                  <  2:6280051/93‑1 (MQ=255)

                                                                        ggCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACt                 <  1:4154677/93‑1 (MQ=255)

                                                                         gCCAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTg                <  2:2656459/93‑1 (MQ=255)

                                                                          ccAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTgg               <  2:732379/93‑1 (MQ=255)

                                                                           cAGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGgt              <  1:7286358/93‑1 (MQ=255)

                                                                            aGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGgtg             >  2:722039/1‑93 (MQ=255)

                                                                            aGTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGgtg             <  1:7874396/93‑1 (MQ=255)

                                                                             gTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGgtgt            >  1:1240118/1‑93 (MQ=255)

                                                                             gTATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGgtgt            >  2:4908670/1‑93 (MQ=255)

                                                                              tATTTCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTc           >  1:794375/1‑93 (MQ=255)

                                                                                tttCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTgg               >  2:3127235/1‑87 (MQ=255)

                                                                                tttCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTgg               <  1:3127236/87‑1 (MQ=255)

                                                                                tttCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCga         <  1:755828/93‑1 (MQ=255)

                                                                                tttCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCga         <  2:652803/93‑1 (MQ=255)

                                                                                tttCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCga         >  1:4070823/1‑93 (MQ=255)

                                                                                  tCCCGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCgatg       <  1:2512361/93‑1 (MQ=255)

                                                                                    ccGCAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCCACTGGTGTCgatgat     >  2:6059967/1‑93 (MQ=255)

                                                                                       cAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCGATGATTaa  >  1:6066182/1‑93 (MQ=255)

                                                                                       cAGCCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCGATGATTaa  >  1:5343031/1‑93 (MQ=255)

GCTGCGCGCAGCCCAAGAGGATATTCGCCAACGACAAGGGCAACTGGCAGGGCAGATTAATATCGGCATGGGGGCCAGTATTTCCCGCAGTCTGATGCCAGCTGTCATATCTCGTTTTCATCAGCAGCATCCGCAGGTAAAAGTACGCATTATGGAAGGGCAACTGGTGTCGATGATTAA  >  NC_000913/1402998‑1403177

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