BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,302,652	A→G	Y513H (TAT→CAT)	yjcS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,302,652	0	A	G	99.5%	609.9 / NA	196	Y513H (TAT→CAT)	yjcS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (78/117); minor base T (1/0); total (79/117)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 207 total aligned reads displayed.

TGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCG  >  NC_000913/4302562‑4302741

tggtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGc                                                                                         <  2:3267405/93‑1 (MQ=255)

 ggtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGcc                                                                                        <  1:949326/93‑1 (MQ=255)

     ccGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGc                                                                                    >  2:3231548/1‑93 (MQ=255)

     ccGGGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGc                                                                                    <  2:2890425/93‑1 (MQ=255)

      cGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCt                                                                                   <  1:3815311/93‑1 (MQ=255)

       gTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTg                                                                                  >  1:4677097/1‑93 (MQ=255)

        tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGt                                                                                 <  2:649400/93‑1 (MQ=255)

        tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGt                                                                                 <  2:1418488/93‑1 (MQ=255)

        tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGCAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGt                                                                                 >  2:1118068/1‑93 (MQ=255)

         ggCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGtt                                                                                >  1:562154/1‑93 (MQ=255)

         ggCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGtt                                                                                <  1:2222257/93‑1 (MQ=255)

          gCTGAACTTATGCACCCCTTCGCGCTGCTCTTTTGCGCCGTTCTGGTAGAAATCGCGCCATGTGTCGGACTCGGCCTGATGGCCCAGCTGTTc                                                                               >  2:4686887/1‑93 (MQ=255)

          gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTc                                                                               >  1:3655770/1‑93 (MQ=255)

          gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTc                                                                               <  1:1564646/93‑1 (MQ=255)

          gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTc                                                                               >  1:579635/1‑93 (MQ=255)

           cTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCa                                                                              <  2:2135456/93‑1 (MQ=255)

            tGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaa                                                                             >  2:4594862/1‑93 (MQ=255)

            tGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaa                                                                             >  1:1629911/1‑93 (MQ=255)

             gAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaaa                                                                            >  2:4531901/1‑93 (MQ=255)

             gAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaaa                                                                            >  2:267271/1‑93 (MQ=255)

             gAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaaa                                                                            <  2:3905789/93‑1 (MQ=255)

             gAACTTATGCACCCCTTCGCGCAGCTCTTTCACGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaaa                                                                            <  1:368523/93‑1 (MQ=255)

              aaCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAg                                                                           <  2:2612819/93‑1 (MQ=255)

               aCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGt                                                                          >  2:4157718/1‑93 (MQ=255)

                cTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGtt                                                                         <  1:1200203/93‑1 (MQ=255)

                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGCCTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  2:1818091/93‑1 (MQ=255)

                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  1:4025906/93‑1 (MQ=255)

                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  1:4180266/93‑1 (MQ=255)

                  tATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTAt                                                                       <  1:3205273/93‑1 (MQ=255)

                    tGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTc                                                                     <  2:2346132/93‑1 (MQ=255)

                       aCCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCt                                                                  <  2:1395902/93‑1 (MQ=255)

                        ccccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCtt                                                                 >  2:531830/1‑93 (MQ=255)

                        ccccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCtt                                                                 <  2:1130051/93‑1 (MQ=255)

                          ccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGc                                                               <  1:4514423/93‑1 (MQ=255)

                          ccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGc                                                               <  2:3718979/93‑1 (MQ=255)

                          ccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTAGCTTGc                                                               >  1:1744984/1‑93 (MQ=255)

                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              <  2:757860/93‑1 (MQ=255)

                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              >  2:49630/1‑93 (MQ=255)

                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              <  2:2818526/93‑1 (MQ=255)

                            ttCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAg                                                             >  2:3550715/1‑93 (MQ=255)

                              cgcgCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAt                                                           >  2:3640049/1‑93 (MQ=255)

                               gcgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAtt                                                          >  2:1195833/1‑93 (MQ=255)

                               gcgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAtt                                                          <  1:1873517/93‑1 (MQ=255)

                                cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttc                                                         <  1:3039297/93‑1 (MQ=255)

                                cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttc                                                         <  2:202167/93‑1 (MQ=255)

                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGTCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        <  2:4234486/93‑1 (MQ=255)

                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        <  1:2848780/93‑1 (MQ=255)

                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        <  2:3066643/93‑1 (MQ=255)

                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        >  1:1947870/1‑93 (MQ=255)

                                   aactctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      >  2:2959950/3‑93 (MQ=255)

                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  2:339542/93‑1 (MQ=255)

                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  1:680017/93‑1 (MQ=255)

                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  1:3083498/93‑1 (MQ=255)

                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      >  2:4097755/1‑93 (MQ=255)

                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  2:3991362/93‑1 (MQ=255)

                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     <  1:3466967/93‑1 (MQ=255)

                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGGGCAAAGTGATTCGCTTGTAGAttcttgg                                                     >  2:2497337/1‑91 (MQ=255)

                                     ctctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcgc                                                    >  1:4378078/1‑93 (MQ=255)

                                       ctTTCGCGCCGGTCATGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGTCCAGCTGTTCAAAGTTATTCGCTTTGAGATTCTTCACGa                                                  >  2:4670808/1‑93 (MQ=255)

                                       ctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcgcgc                                                  >  2:163980/1‑92 (MQ=255)

                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  2:2213584/93‑1 (MQ=255)

                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  1:2750568/93‑1 (MQ=255)

                                            gcgcCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTgg                                             <  2:1703724/93‑1 (MQ=255)

                                                cGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCAc                                         >  1:2587623/1‑93 (MQ=255)

                                                    cAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCgg                                     >  1:763322/1‑93 (MQ=255)

                                                      ggTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGt                                   <  2:1535311/93‑1 (MQ=255)

                                                      ggTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGt                                   <  1:4336043/93‑1 (MQ=255)

                                                      ggTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGt                                   >  2:3330133/1‑93 (MQ=255)

                                                       gTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGtt                                  >  1:2480518/1‑93 (MQ=255)

                                                            aaCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcgg                             <  1:2013053/93‑1 (MQ=255)

                                                            aaCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcgg                             >  1:2122821/1‑93 (MQ=255)

                                                            aaCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcgg                             >  2:3435290/1‑93 (MQ=255)

                                                            aaCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcgg                             >  1:2805232/1‑93 (MQ=255)

                                                             aCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcggc                            >  2:1908076/1‑93 (MQ=255)

                                                             aCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcggc                            <  1:2811735/93‑1 (MQ=255)

                                                               cgcgCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCaa                          <  1:1810362/93‑1 (MQ=255)

                                                                gcgcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAAt                         <  1:2171361/93‑1 (MQ=255)

                                                                 cgcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATc                        <  1:698345/93‑1 (MQ=255)

                                                                 cgcCATGTGCCGGACTCGGCCTGATGCCCAAGCTGTCAAAACTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCGGGGTTGGCGGCAATc                        <  2:2981216/93‑1 (MQ=255)

                                                                  gcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCa                       <  2:2539017/93‑1 (MQ=255)

                                                                   ccATGTGGCGGACTCGTCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAc                      <  2:596436/93‑1 (MQ=255)

                                                                     aTGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCt                    >  2:2451464/1‑93 (MQ=255)

                                                                          gCGGACTCGGCCTGATGGCCCAGCTGTTCAAATTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTc               <  2:159954/93‑1 (MQ=255)

                                                                            ggACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTACAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcag             <  2:641657/93‑1 (MQ=255)

                                                                               cTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagcag          >  2:3531648/1‑93 (MQ=255)

                                                                               cTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagcag          >  2:2385112/1‑93 (MQ=255)

                                                                               cTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagcag          <  1:4403294/93‑1 (MQ=255)

                                                                                tCGTCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGt         <  2:4422717/93‑1 (MQ=255)

                                                                                tCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGGGACCTCATCACACGGGCTAGCGGCCATCACCTGTTTGAGCAGt         >  1:793618/1‑93 (MQ=255)

                                                                                tCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGt         <  1:3320204/93‑1 (MQ=255)

                                                                                 cGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGtt        >  2:3839801/1‑93 (MQ=255)

                                                                                   gCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCt      <  2:856175/93‑1 (MQ=255)

                                                                                   gCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCt      >  2:3993957/1‑93 (MQ=255)

                                                                                   gCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCt      <  1:3403247/93‑1 (MQ=255)

                                                                                     cTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGc    >  1:3613041/1‑93 (MQ=255)

                                                                                     cTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGc    >  2:4215150/1‑93 (MQ=255)

                                                                                      tGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGcc   <  1:1022406/93‑1 (MQ=255)

                                                                                      tGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGcc   >  1:2357863/1‑93 (MQ=255)

                                                                                      tGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGcc   >  2:379172/1‑93 (MQ=255)

                                                                                       gATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCg  <  2:1585034/93‑1 (MQ=255)

TGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCG  >  NC_000913/4302562‑4302741

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