BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,647,298	T→C	T1004A (ACG→GCG)	b2520 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,647,298	0	T	C	99.6%	713.2 / ‑5.4	222	T1004A (ACG→GCG)	b2520	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (133/88); total (133/89)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 236 total aligned reads displayed.

CACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCG  >  NC_000913/2647208‑2647387

cACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGccc                                                                                         >  1:2361175/1‑93 (MQ=255)

 ccTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCa                                                                                        <  1:294706/92‑1 (MQ=255)

 aCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCa                                                                                        >  2:5704669/1‑93 (MQ=255)

 aCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCa                                                                                        >  2:3201164/1‑93 (MQ=255)

  cTCGCAATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAg                                                                                       >  1:3408961/1‑93 (MQ=255)

    cGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGc                                                                                     <  2:1396818/93‑1 (MQ=255)

    cGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGc                                                                                     <  2:1813385/93‑1 (MQ=255)

     gCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCt                                                                                    <  2:6670659/93‑1 (MQ=255)

      ccATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCtt                                                                                   >  2:5756449/1‑93 (MQ=255)

      ccATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCtt                                                                                   >  2:5739160/1‑93 (MQ=255)

        aTAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTgc                                                                                 >  1:5569027/1‑93 (MQ=255)

        aTAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTgc                                                                                 >  1:302685/1‑93 (MQ=255)

         tAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTgcg                                                                                >  2:113978/1‑93 (MQ=255)

            aGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCa                                                                             >  2:4778011/1‑93 (MQ=255)

            aGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCa                                                                             <  2:3897262/93‑1 (MQ=255)

             gcgcGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCAt                                                                            <  1:2572815/93‑1 (MQ=255)

             gcgcGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCAt                                                                            >  2:6076237/1‑93 (MQ=255)

               gcgGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGa                                                                          <  1:5900769/93‑1 (MQ=255)

                 ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGAcg                                                                        <  1:350421/93‑1 (MQ=255)

                 ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGAcg                                                                        >  2:1042581/1‑93 (MQ=255)

                    aTGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACgcgc                                                                     >  1:6085756/1‑93 (MQ=255)

                      gttcagttcaGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCaa                                                                   >  2:2638778/1‑93 (MQ=255)

                      gttcagttcaGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCaa                                                                   >  1:648504/1‑93 (MQ=255)

                        tcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAtt                                                                 <  1:3151642/93‑1 (MQ=255)

                        tcagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAAtt                                                                 >  2:556071/1‑93 (MQ=255)

                         cagttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTc                                                                <  2:2324813/93‑1 (MQ=255)

                          agttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCa                                                               >  2:19711/1‑93 (MQ=255)

                           gttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCAc                                                              >  1:4911321/1‑93 (MQ=255)

                           gttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCAc                                                              <  2:6733212/93‑1 (MQ=255)

                            ttcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCAcc                                                             <  1:6652708/93‑1 (MQ=255)

                             tcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCg                                                            >  1:4149130/1‑93 (MQ=255)

                             tcagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCg                                                            <  1:2313530/93‑1 (MQ=255)

                              cagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGt                                                           <  2:4536129/93‑1 (MQ=255)

                              cagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGt                                                           <  1:1085249/93‑1 (MQ=255)

                              cagCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGt                                                           <  1:1515676/93‑1 (MQ=255)

                               agCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGtt                                                          <  2:3688795/93‑1 (MQ=255)

                                 cAATCACCGGTGCGGCAACTATCACTTTACTTTCATTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTaa                                                        <  2:1758123/93‑1 (MQ=255)

                                  aaTCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTaaa                                                       >  1:4572134/1‑93 (MQ=255)

                                   aTCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTaaaa                                                      <  1:4098169/93‑1 (MQ=255)

                                    tCACGGGTACGGCACCAATCACTTTGCTTTCGTTGCGACAGAAGTCATCTGCCGCCCAGGCTTGCGCCGTGACGCGCAATTCACCGTTAAAAt                                                     <  1:1954743/93‑1 (MQ=255)

                                      aCCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCg                                                   <  2:2933825/93‑1 (MQ=255)

                                          gTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCaa                                               <  2:4367592/93‑1 (MQ=255)

                                          gTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCaa                                               >  1:286023/1‑93 (MQ=255)

                                           tGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAt                                              >  2:5184258/1‑93 (MQ=255)

                                            gCGGTAACTATCACTTTACGTTCGTTGCTACCGAAGTCATCTTCCGCCCAGGCTTGCGCCATGACGCTCAATTCACCGTTAAAGGCGCCAATc                                             <  1:5516379/93‑1 (MQ=255)

                                            gCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATc                                             <  1:4877553/93‑1 (MQ=255)

                                            gCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATc                                             <  2:5202763/93‑1 (MQ=255)

                                            gCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATc                                             <  2:4562039/93‑1 (MQ=255)

                                             cGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCg                                            <  2:4866317/93‑1 (MQ=255)

                                              ggCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCgg                                           >  2:5813441/1‑93 (MQ=255)

                                               gCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGc                                          >  2:3604766/1‑93 (MQ=255)

                                               gCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGc                                          <  2:4629087/93‑1 (MQ=255)

                                                  aCTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAgt                                       <  2:3095301/93‑1 (MQ=255)

                                                    tATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCTATCGGCAgtgt                                     >  1:2476904/1‑93 (MQ=255)

                                                    tATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAgtgt                                     >  1:4736443/1‑93 (MQ=255)

                                                     aTCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTa                                    <  1:5725413/93‑1 (MQ=255)

                                                      tCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTaa                                   <  2:6322432/93‑1 (MQ=255)

                                                       cACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAc                                  <  1:2134573/93‑1 (MQ=255)

                                                       cACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAc                                  >  1:5002671/1‑93 (MQ=255)

                                                       cACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAc                                  >  1:677749/1‑93 (MQ=255)

                                                        actttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAcc                                 <  1:3786898/93‑1 (MQ=255)

                                                        actttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAcc                                 >  1:3590545/1‑93 (MQ=255)

                                                        actttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAAcc                                 >  2:54680/1‑93 (MQ=255)

                                                         ctttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCg                                >  1:4289666/1‑93 (MQ=255)

                                                         ctttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCg                                <  1:4986107/93‑1 (MQ=255)

                                                         ctttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCCGCAGTGTAAGCg                                >  1:209086/1‑93 (MQ=255)

                                                          tttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGa                               >  2:4921667/1‑93 (MQ=255)

                                                          tttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGa                               <  2:4308033/93‑1 (MQ=255)

                                                          tttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGa                               <  1:945343/93‑1 (MQ=255)

                                                          tttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGa                               >  2:6539816/1‑93 (MQ=255)

                                                          tttacgttCGTTGCTACCGAAGTCAGCTGCCGCCCGGGCTTGCGCCATGACGCGCAATTCACCGTTAAAAGCGCCAATCGGCAGTGTATCCGa                               <  1:1838911/93‑1 (MQ=255)

                                                           ttactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTTACCGAg                              >  1:5105909/1‑93 (MQ=255)

                                                            tactttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGc                             >  1:5729014/1‑93 (MQ=255)

                                                             tctttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGcc                            <  1:6222810/92‑1 (MQ=255)

                                                             actttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGcc                            >  2:1583564/1‑93 (MQ=255)

                                                             actttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGcc                            >  2:2780674/1‑93 (MQ=255)

                                                              ctttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCt                           >  1:4774453/1‑93 (MQ=255)

                                                              ctttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCt                           >  2:3873499/1‑93 (MQ=255)

                                                              ctttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCt                           >  2:806826/1‑93 (MQ=255)

                                                               tttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCtt                          >  1:2670364/1‑93 (MQ=255)

                                                                ttCGTTGCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTc                         >  2:3323692/1‑93 (MQ=255)

                                                                  cGTTGCTACCGAAGTCATCTGCCGCCCAGGATTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGACATCGc                       >  1:402942/1‑93 (MQ=255)

                                                                      gCTACCGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTg                   >  2:5117318/1‑93 (MQ=255)

                                                                        tACCGAAGTCATCTGCAGCCCAGGCTTGCTCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGtt                 <  1:1378634/93‑1 (MQ=255)

                                                                           cGAAGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGt              <  2:3002036/93‑1 (MQ=255)

                                                                             aaGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTg            >  1:4261165/1‑93 (MQ=255)

                                                                              aGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGa           <  1:529118/93‑1 (MQ=255)

                                                                              aGTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGa           <  1:4716696/93‑1 (MQ=255)

                                                                               gTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAg          <  1:3816473/93‑1 (MQ=255)

                                                                               gTCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAg          >  1:1624693/1‑93 (MQ=255)

                                                                                tCATCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGc         <  1:1912030/93‑1 (MQ=255)

                                                                                  aTCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGt       >  1:6184701/1‑93 (MQ=255)

                                                                                   tCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGtt      >  2:3129208/1‑93 (MQ=255)

                                                                                   tCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGtt      >  1:3926212/1‑93 (MQ=255)

                                                                                   tCTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGtt      <  1:5940460/93‑1 (MQ=255)

                                                                                    cTGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTa     >  2:3635076/1‑93 (MQ=255)

                                                                                     tGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTTAGCGTTAc    >  1:1474740/1‑93 (MQ=255)

                                                                                     tGCCGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTAc    <  2:4149253/93‑1 (MQ=255)

                                                                                      gccgccCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTAcc   <  1:105596/93‑1 (MQ=255)

                                                                                       ccGCCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCg  >  2:848666/1‑93 (MQ=255)

CACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCG  >  NC_000913/2647208‑2647387

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