BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,531,089	T→C	I88M (ATA→ATG)	yjhQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,531,089	0	T	C	89.3%	791.0 / 68.5	272	I88M (ATA→ATG)	yjhQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (11/18); new base C (109/134); total (120/152)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 281 total aligned reads displayed.

AGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGC  >  NC_000913/4530997‑4531178

agaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGt                                                                                           <  1:1286375/93‑1 (MQ=255)

  aaCAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCAt                                                                                         >  2:6415956/1‑93 (MQ=255)

  aaCAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCAt                                                                                         <  2:2806465/93‑1 (MQ=255)

    cAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGa                                                                                       >  2:2004260/1‑93 (MQ=255)

     aaaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGAc                                                                                      <  2:5238437/93‑1 (MQ=255)

      aaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGAcc                                                                                     >  1:6327137/1‑93 (MQ=255)

        aaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGc                                                                                   >  2:5799758/1‑93 (MQ=255)

        aaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGc                                                                                   >  1:4519758/1‑93 (MQ=255)

         aCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGCGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCTGCATGACCGcc                                                                                  >  1:2498414/1‑93 (MQ=255)

          ccGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCa                                                                                 >  1:3486961/1‑93 (MQ=255)

            gTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAgg                                                                               >  1:3364587/1‑93 (MQ=255)

             tCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAggg                                                                              <  2:903710/93‑1 (MQ=255)

             tCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAggg                                                                              <  1:2231308/93‑1 (MQ=255)

              cTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAgggg                                                                             >  2:2475934/1‑93 (MQ=255)

              cTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAgggg                                                                             <  2:1100784/93‑1 (MQ=255)

              cTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAgggg                                                                             <  1:1139220/93‑1 (MQ=255)

               tGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGa                                                                            <  2:385299/93‑1 (MQ=255)

                ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAg                                                                           <  2:1811624/93‑1 (MQ=255)

                ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAg                                                                           >  1:3149187/1‑93 (MQ=255)

                 gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGc                                                                          <  1:1038914/93‑1 (MQ=255)

                 gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGc                                                                          >  1:3687475/1‑93 (MQ=255)

                 gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGc                                                                          <  1:1031430/93‑1 (MQ=255)

                    gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAg                                                                       >  1:1907528/1‑93 (MQ=255)

                    gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAg                                                                       <  2:3562531/93‑1 (MQ=255)

                    gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAg                                                                       >  1:711247/1‑93 (MQ=255)

                     cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAgg                                                                      <  2:5330164/93‑1 (MQ=255)

                      ccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCCTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGa                                                                     >  1:211513/1‑93 (MQ=255)

                       cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGAt                                                                    >  2:3024241/1‑93 (MQ=255)

                       cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGAt                                                                    >  1:963937/1‑93 (MQ=255)

                        aTCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGAtg                                                                   >  1:1715292/1‑93 (MQ=255)

                         tCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCTGCATGACCGCCAGGGGCGCGAGGAtgt                                                                  >  1:4620079/1‑93 (MQ=255)

                         tCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGAtgt                                                                  >  2:2300908/1‑93 (MQ=255)

                           aaCCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGc                                                                >  1:2844133/1‑93 (MQ=255)

                              cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATc                                                             <  2:6494892/93‑1 (MQ=255)

                               tCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCa                                                            >  2:4315597/1‑93 (MQ=255)

                                 aGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAAc                                                          >  1:5893454/1‑93 (MQ=255)

                                  ggTTCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACg                                                         <  1:4113100/93‑1 (MQ=255)

                                  ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACg                                                         >  1:5504399/1‑93 (MQ=255)

                                   gTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACgg                                                        <  2:1441373/93‑1 (MQ=255)

                                    tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGt                                                       <  2:101590/93‑1 (MQ=255)

                                     gCTCGATGCCAGTGCGTATCTACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTg                                                      <  2:4502325/93‑1 (MQ=255)

                                       tCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGaa                                                    <  1:158695/93‑1 (MQ=255)

                                        cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAAt                                                   <  2:5628302/93‑1 (MQ=255)

                                        cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAAt                                                   >  1:5268371/1‑93 (MQ=255)

                                        cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAAt                                                   <  1:3558742/93‑1 (MQ=255)

                                          aTGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATcc                                                 >  2:2987599/1‑93 (MQ=255)

                                          aTGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTCAATcc                                                 >  1:854929/1‑93 (MQ=255)

                                               aGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAttta                                            >  2:6382386/1‑92 (MQ=255)

                                               aGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtc                                            >  2:4669219/1‑93 (MQ=255)

                                               aGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtc                                            >  1:2701426/1‑93 (MQ=255)

                                                 tGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtctc                                          >  2:2632202/1‑93 (MQ=255)

                                                   cGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCCAGGATGTGCATCAACGGTGAATCCATTTCTCCt                                        <  1:1110497/93‑1 (MQ=255)

                                                    gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCtt                                       <  2:494796/93‑1 (MQ=255)

                                                     tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                                      <  1:4761936/93‑1 (MQ=255)

                                                     tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                                      >  1:4206277/1‑93 (MQ=255)

                                                     tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                                      <  2:1897640/93‑1 (MQ=255)

                                                       tCACCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                                    <  2:2504189/93‑1 (MQ=255)

                                                       tCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGAGAATCCATTTCTCCTTTGa                                    >  2:1539333/1‑93 (MQ=255)

                                                        cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGaa                                   <  1:5691190/93‑1 (MQ=255)

                                                         aaCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAg                                  <  2:1723969/93‑1 (MQ=255)

                                                           ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                <  2:1950873/93‑1 (MQ=255)

                                                           ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                <  2:4359989/93‑1 (MQ=255)

                                                           ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                <  2:2957172/93‑1 (MQ=255)

                                                             gTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTg                              >  2:6631833/1‑93 (MQ=255)

                                                              tcctccCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGc                             >  2:3232321/1‑93 (MQ=255)

                                                               cctccCACGCCCATCCCGTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                            <  2:2170748/93‑1 (MQ=255)

                                                                ctccCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                           >  1:4022965/1‑93 (MQ=255)

                                                                ctccCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                           >  2:4165443/1‑93 (MQ=255)

                                                                  cccACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCgg                         >  2:6123752/1‑93 (MQ=255)

                                                                   ccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCtgg                        >  2:598032/1‑93 (MQ=255)

                                                                   ccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCggg                        >  2:1342760/1‑93 (MQ=255)

                                                                   ccACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCggg                        >  2:3726147/1‑93 (MQ=255)

                                                                    cACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGt                       <  2:1190855/93‑1 (MQ=255)

                                                                     aCGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                      >  2:4624035/1‑93 (MQ=255)

                                                                      cGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaa                     <  2:5093960/93‑1 (MQ=255)

                                                                       gCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa                    <  1:2579912/93‑1 (MQ=255)

                                                                        cccATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat                   >  1:5192543/1‑93 (MQ=255)

                                                                         ccATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaata                  <  1:6620175/93‑1 (MQ=255)

                                                                          cATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataa                 >  2:606260/1‑93 (MQ=255)

                                                                          cATCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataa                 >  2:3347118/1‑93 (MQ=255)

                                                                            tCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaa               >  1:6467316/1‑93 (MQ=255)

                                                                            tCCCTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaa               <  2:5725218/93‑1 (MQ=255)

                                                                             cccTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaat              >  1:3077397/1‑93 (MQ=255)

                                                                              ccTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaata             >  2:271245/1‑93 (MQ=255)

                                                                              ccTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaata             >  1:4416381/1‑93 (MQ=255)

                                                                              ccTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaata             <  2:6049859/93‑1 (MQ=255)

                                                                               cTTGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAatat            >  1:3970264/1‑93 (MQ=255)

                                                                                ttGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATg           <  1:5755192/93‑1 (MQ=255)

                                                                                ttGGTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATg           >  2:5482306/1‑93 (MQ=255)

                                                                                  ggTATTCCGGTATGACCGCCAAGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGcc         <  2:739638/93‑1 (MQ=255)

                                                                                  ggTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGcc         >  1:5088841/1‑93 (MQ=255)

                                                                                   gTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGccc        >  1:6152240/1‑93 (MQ=255)

                                                                                   gTATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGccc        <  1:5569774/93‑1 (MQ=255)

                                                                                    tATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg       <  1:6221754/93‑1 (MQ=255)

                                                                                    tATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg       <  2:61978/93‑1 (MQ=255)

                                                                                    tATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg       <  1:2799120/93‑1 (MQ=255)

                                                                                    tATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg       <  1:2614190/93‑1 (MQ=255)

                                                                                    tATTCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg       <  2:2105922/93‑1 (MQ=255)

                                                                                       tCCGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACg    >  2:807211/1‑93 (MQ=255)

                                                                                        ccGGCATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACgg   >  2:3211170/1‑93 (MQ=255)

AGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGC  >  NC_000913/4530997‑4531178

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