BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,630,942	T→C	K383K (AAA→AAG)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,630,942	0	T	C	99.6%	783.5 / NA	239	K383K (AAA→AAG)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (131/107); total (132/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

CACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCC  >  NC_000913/2630852‑2631031

 aCCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTg                                                                                        >  1:2954415/1‑93 (MQ=255)

  cctgTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGt                                                                                       <  2:3036249/90‑1 (MQ=255)

  ccAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGt                                                                                       <  1:4484189/93‑1 (MQ=255)

  ccAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGt                                                                                       >  2:1337252/1‑93 (MQ=255)

  ccAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGt                                                                                       <  2:4721235/93‑1 (MQ=255)

    aGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTaa                                                                                     >  1:4942983/1‑93 (MQ=255)

    aGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTaa                                                                                     <  1:2517651/93‑1 (MQ=255)

    aGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTaa                                                                                     >  2:4857170/1‑93 (MQ=255)

     gTTTGTCGGCGGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCTTACCACGGTAAGACTTGTAAg                                                                                    <  2:4010815/93‑1 (MQ=255)

     gTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAg                                                                                    <  2:5484563/93‑1 (MQ=255)

      tttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGa                                                                                   >  2:5948199/1‑93 (MQ=255)

      tttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGa                                                                                   <  2:3642219/93‑1 (MQ=255)

      tttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACAGCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGa                                                                                   <  2:442820/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  1:2322092/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  2:2722646/93‑1 (MQ=255)

        tGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGCAc                                                                                 >  2:1889782/1‑93 (MQ=255)

        tGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAAc                                                                                 <  2:2925985/93‑1 (MQ=255)

         gTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACg                                                                                <  1:3094853/93‑1 (MQ=255)

          gcGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACgg                                                                               <  1:3808776/92‑1 (MQ=255)

             gCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTGTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGccc                                                                            <  1:4640259/93‑1 (MQ=255)

              cAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCt                                                                           >  2:3938070/1‑93 (MQ=255)

               aGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTg                                                                          <  1:3503123/93‑1 (MQ=255)

                gCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTgg                                                                         >  1:3939519/1‑93 (MQ=255)

                gCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTgg                                                                         >  2:3270219/1‑93 (MQ=255)

                   ttATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTag                                                                      >  1:3309011/1‑93 (MQ=255)

                     aTCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTagag                                                                    <  2:2587377/93‑1 (MQ=255)

                     aTCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTagag                                                                    >  1:5793040/1‑93 (MQ=255)

                       cGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGtt                                                                  <  1:4098328/93‑1 (MQ=255)

                        gCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTc                                                                 >  2:80954/1‑93 (MQ=255)

                        gCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTc                                                                 >  1:5180876/1‑93 (MQ=255)

                        gCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTc                                                                 >  2:3602517/1‑93 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               >  2:2680360/1‑93 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               >  1:1543195/1‑93 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               >  1:1782894/1‑93 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              >  1:4524292/1‑93 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              <  1:5843007/93‑1 (MQ=255)

                            tGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAtt                                                             >  1:5556911/1‑93 (MQ=255)

                            tGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAtt                                                             <  2:463509/93‑1 (MQ=255)

                             ggAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAttt                                                            >  2:6173832/1‑93 (MQ=255)

                             ggAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGACCCTGGGAGAGGTCGAttt                                                            >  2:1298572/1‑93 (MQ=255)

                              gAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTc                                                           <  1:4994343/93‑1 (MQ=255)

                               aaaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCa                                                          >  1:384161/1‑93 (MQ=255)

                                aaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAc                                                         <  1:5993389/93‑1 (MQ=255)

                                aaTAACGGTCAGAGCAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAc                                                         <  1:4264860/93‑1 (MQ=255)

                                 aTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAcc                                                        <  1:4059017/93‑1 (MQ=255)

                                 aTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAcc                                                        >  1:6123460/1‑93 (MQ=255)

                                  tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAccc                                                       <  2:4341072/93‑1 (MQ=255)

                                  aaaCGCTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGGAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCcccc                                                       <  2:2262400/92‑1 (MQ=255)

                                   aaCGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCg                                                      >  1:5655938/1‑93 (MQ=255)

                                    aCGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCgg                                                     >  1:1604436/1‑93 (MQ=255)

                                     cGGTCAGAGCAACCTTTGGCCATCGCGCTCAGGGAACCCATACCACGGTAAGACTTGTAAGACCGGCCCTGGTAGAGTTCGATTTCACCCCga                                                    <  2:3441245/93‑1 (MQ=255)

                                        tCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAt                                                 >  2:3819258/1‑93 (MQ=255)

                                        tCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAt                                                 >  1:1033549/1‑93 (MQ=255)

                                         cAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAtt                                                >  2:1758407/1‑93 (MQ=255)

                                          agagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttc                                               <  2:2985392/93‑1 (MQ=255)

                                          agagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttc                                               <  2:2438926/93‑1 (MQ=255)

                                           gagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                              >  2:391145/1‑93 (MQ=255)

                                               aaCCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAg                                          <  1:635614/93‑1 (MQ=255)

                                                aCCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGt                                         >  2:5287989/1‑93 (MQ=255)

                                                aCCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGt                                         <  2:5919139/93‑1 (MQ=255)

                                                 ccTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTa                                        >  1:2849493/1‑93 (MQ=255)

                                                  cTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAc                                       >  2:3183575/1‑93 (MQ=255)

                                                   tttGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAcc                                      >  1:2119964/1‑93 (MQ=255)

                                                    ttGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAccc                                     <  1:3279956/93‑1 (MQ=255)

                                                     tGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCg                                    >  2:4573106/1‑93 (MQ=255)

                                                      ggACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGc                                   <  1:2806550/93‑1 (MQ=255)

                                                        aCATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCa                                 <  2:1057140/93‑1 (MQ=255)

                                                        aCATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCa                                 <  2:5103363/93‑1 (MQ=255)

                                                           tCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCa                              <  2:1558134/93‑1 (MQ=255)

                                                           tCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCa                              <  1:4000571/93‑1 (MQ=255)

                                                           tCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCa                              <  1:6209024/93‑1 (MQ=255)

                                                            cgcgCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCAt                             >  2:2779889/1‑93 (MQ=255)

                                                             gcgcCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATg                            >  2:3105052/1‑93 (MQ=255)

                                                               gcCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGa                          <  1:2651253/93‑1 (MQ=255)

                                                                cccAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGaa                         >  2:2893210/1‑93 (MQ=255)

                                                                  cAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAcc                       >  1:1262023/1‑93 (MQ=255)

                                                                     ggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAc                    <  2:1127100/93‑1 (MQ=255)

                                                                     ggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCCAc                    >  2:3880350/1‑93 (MQ=255)

                                                                      gAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAcc                   >  2:584072/1‑93 (MQ=255)

                                                                       aaCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCa                  <  1:102050/93‑1 (MQ=255)

                                                                       aaCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCa                  <  1:1825631/93‑1 (MQ=255)

                                                                       aaCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCa                  <  1:5134419/93‑1 (MQ=255)

                                                                       aaCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCa                  <  1:5097348/93‑1 (MQ=255)

                                                                          ccATCCCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCa               <  2:4172184/93‑1 (MQ=255)

                                                                          ccATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCa               >  1:5443169/1‑93 (MQ=255)

                                                                          ccATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCa               >  1:4242392/1‑93 (MQ=255)

                                                                           cATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAc              >  1:2229298/1‑93 (MQ=255)

                                                                           cATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAc              <  2:4906639/93‑1 (MQ=255)

                                                                            aTACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAcc             >  2:1454742/1‑93 (MQ=255)

                                                                            aTACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAcc             <  1:2729270/93‑1 (MQ=255)

                                                                             tACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACcg            >  2:5349050/1‑93 (MQ=255)

                                                                             tACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCAGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACcg            <  1:532671/93‑1 (MQ=255)

                                                                              aCCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACcgc           <  1:3156294/93‑1 (MQ=255)

                                                                                 aCGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCt        <  1:6042862/93‑1 (MQ=255)

                                                                                    gTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgc     <  2:5204658/93‑1 (MQ=255)

                                                                                    gTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgc     <  1:5256720/93‑1 (MQ=255)

                                                                                    gTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgc     >  2:5556478/1‑93 (MQ=255)

                                                                                     tAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgcg    >  1:3694274/1‑93 (MQ=255)

                                                                                       aGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGcc  >  2:579491/1‑93 (MQ=255)

                                                                                       aGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGcc  <  2:2571451/93‑1 (MQ=255)

CACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCC  >  NC_000913/2630852‑2631031

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