BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	646,808	(C)_6→7	coding (702/879 nt)	citG ←	triphosphoribosyl‑dephospho‑CoA synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,802	1	.	C	97.6%	1008.5 / 14.1	250	coding (708/879 nt)	citG	triphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand): ref base . (1/5); new base C (99/145); total (100/150)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 274 total aligned reads displayed.

TACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646702‑646897

  cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑cc                                                                                                <  2:1294251/101‑1 (MQ=255)

    cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑cccc                                                                                              <  2:1293933/101‑1 (MQ=255)

       ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cc                                                                                          <  2:1294194/101‑2 (MQ=255)

       ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cc                                                                                          <  2:1294339/101‑2 (MQ=255)

        cTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cct                                                                                         <  1:1294622/101‑3 (MQ=255)

        cTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cct                                                                                         <  1:1294530/101‑3 (MQ=255)

         tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑T‑‑ccccc‑cctc                                                                                        <  1:1294324/101‑4 (MQ=255)

         tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cctc                                                                                        <  1:1294176/101‑4 (MQ=255)

         tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cctc                                                                                        <  1:1294369/101‑4 (MQ=255)

            cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGcc                                                                                      <  1:1294254/101‑1 (MQ=255)

              aaCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAc                                                                                    <  1:1294400/101‑1 (MQ=255)

                      gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGAt                                                                            <  1:1294286/101‑1 (MQ=255)

                         aTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCa                                                                         <  1:1294249/101‑1 (MQ=255)

                          tAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCaa                                                                        <  2:1294304/101‑1 (MQ=255)

                           aaTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAAc                                                                       <  2:1294423/101‑1 (MQ=255)

                            aTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACg                                                                      <  1:1294144/101‑1 (MQ=255)

                               gagaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTg                                                                   <  2:1293886/101‑1 (MQ=255)

                               gagaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTg                                                                   <  1:1294499/101‑1 (MQ=255)

                               gagaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTg                                                                   <  1:1294375/101‑1 (MQ=255)

                                    cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATc                                                              <  1:1294488/101‑1 (MQ=255)

                                    cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATc                                                              <  2:1294199/101‑1 (MQ=255)

                                     ggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCg                                                             >  1:1294935/1‑101 (MQ=255)

                                      gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGc                                                            <  1:1294444/101‑1 (MQ=255)

                                      gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTAGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGc                                                            >  2:1294937/1‑101 (MQ=255)

                                       cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGcc                                                           >  2:1294939/1‑101 (MQ=255)

                                        gggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                          <  2:1294681/101‑1 (MQ=255)

                                        gggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                          <  2:1294046/101‑1 (MQ=255)

                                         ggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                         <  1:1294873/101‑1 (MQ=255)

                                          gggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                        <  2:1294126/101‑1 (MQ=255)

                                          gggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                        >  1:1294942/1‑101 (MQ=255)

                                          gggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                        <  1:1294208/101‑1 (MQ=255)

                                           ggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑CCCCCGCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                       <  1:1294347/101‑1 (MQ=255)

                                            gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAcc                                                                                   >  1:1294947/1‑72 (MQ=255)

                                            gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGAt                                                                            >  2:1294945/1‑79 (MQ=255)

                                            gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                      <  2:1294279/101‑1 (MQ=255)

                                            gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                      <  2:1294399/101‑1 (MQ=255)

                                             ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAcc                                                                                   >  2:1294949/1‑71 (MQ=255)

                                                gAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAcc                                                                                   >  1:1294954/1‑68 (MQ=255)

                                                gAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGc                                                  <  1:1293875/101‑1 (MQ=255)

                                                 aaTGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATc                                                              >  2:1294958/1‑88 (MQ=255)

                                                 aaTGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGcc                                                 >  1:1294957/1‑101 (MQ=255)

                                                  aTGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                                <  1:1294380/101‑1 (MQ=255)

                                                  aTGCCCCCGTTTTGCAATAAGGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAATGTTGGTATCGCCGTTGATCGCCa                                                >  1:1294962/1‑101 (MQ=255)

                                                    gCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATc                                              <  1:1294762/101‑1 (MQ=255)

                                                    gCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑ccccc‑cc                                                                                          >  2:1294966/1‑55 (MQ=255)

                                                     cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                             >  1:1294967/1‑101 (MQ=255)

                                                     cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                             <  2:1294393/101‑1 (MQ=255)

                                                       cccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACcgc                                                                                 >  1:1294972/1‑63 (MQ=255)

                                                          tttttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                        >  2:1294974/1‑101 (MQ=255)

                                                           ttttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGAt                                                                            >  1:1294976/1‑64 (MQ=255)

                                                            tttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCCTCAGTAGGAGc                                     >  1:1294978/1‑101 (MQ=255)

                                                             ttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGc                                     <  1:1294307/101‑1 (MQ=255)

                                                              tGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAAc                                                                       >  1:1294982/1‑66 (MQ=255)

                                                               gCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                   >  2:1294985/1‑101 (MQ=255)

                                                                cAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAg                                  >  1:1294988/1‑101 (MQ=255)

                                                                 aataatGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACcgc                                                                                 >  2:1294994/1‑53 (MQ=255)

                                                                 aataatGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                                 >  2:1294992/1‑101 (MQ=255)

                                                                    aatGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                              <  2:1294412/101‑1 (MQ=255)

                                                                    aatGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                              >  1:1294999/1‑101 (MQ=255)

                                                                     atGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                             <  2:1294627/101‑1 (MQ=255)

                                                                     atGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                             >  2:1295001/1‑101 (MQ=255)

                                                                      tGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAcc                                                                                   >  2:1295004/1‑46 (MQ=255)

                                                                       gTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                           <  1:1294318/101‑1 (MQ=255)

                                                                        tttGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAg                          >  1:1295007/1‑101 (MQ=255)

                                                                         ttGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                         <  2:1294705/101‑1 (MQ=255)

                                                                         ttGCGCCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                         <  2:1294456/101‑1 (MQ=255)

                                                                           gcgcCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCaa                       <  1:1294465/101‑1 (MQ=255)

                                                                            cgcCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAcc                                                                                   >  2:1295010/1‑40 (MQ=255)

                                                                            cgcCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                          <  2:1294217/65‑1 (MQ=255)

                                                                             gcCTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATg                     <  1:1294664/101‑1 (MQ=255)

                                                                              ccTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCAcc                                                                                   >  2:1295012/1‑38 (MQ=255)

                                                                               cTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc                   <  1:1294569/101‑1 (MQ=255)

                                                                               cTCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc                   >  1:1295015/1‑101 (MQ=255)

                                                                                tCGCGCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa                  <  2:1294749/101‑1 (MQ=255)

                                                                                 cgcgCTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAg                 >  2:1295017/1‑101 (MQ=255)

                                                                                   cgcTGTAGCCAGCGCAGTC‑‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt               <  1:1294408/101‑1 (MQ=255)

                                                                                   cgcTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt               >  2:1295021/1‑101 (MQ=255)

                                                                                    gcTGTAGCCAGCGCAG‑GCCCCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt               <  2:1294482/101‑1 (MQ=255)

                                                                                    gcTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACcgcgcg                                                                              >  1:1295026/1‑37 (MQ=255)

                                                                                    gcTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGt                                                                 >  2:1295025/1‑50 (MQ=255)

                                                                                    gcTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc              <  1:1294001/101‑1 (MQ=255)

                                                                                     cTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa             <  1:1294361/101‑1 (MQ=255)

                                                                                     cTGTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa             >  2:1295028/1‑101 (MQ=255)

                                                                                       gTAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg           <  2:1294564/101‑1 (MQ=255)

                                                                                        tAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCATCGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa          >  1:1295031/1‑101 (MQ=255)

                                                                                        tAGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa          <  2:1294533/101‑1 (MQ=255)

                                                                                         aGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         <  1:1294718/101‑1 (MQ=255)

                                                                                         aGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         <  1:1294421/101‑1 (MQ=255)

                                                                                         aGCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         >  2:1295034/1‑101 (MQ=255)

                                                                                          gCCAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc        >  1:1295037/1‑101 (MQ=255)

                                                                                           ccAGCGCAGTC‑‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       <  1:1294548/101‑1 (MQ=255)

                                                                                           ccAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       <  2:1294439/101‑1 (MQ=255)

                                                                                           ccAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       <  1:1294349/101‑1 (MQ=255)

                                                                                           ccAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt       >  2:1295040/1‑101 (MQ=255)

                                                                                            cAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa      <  2:1294520/101‑1 (MQ=255)

                                                                                            cAGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa      <  2:1294460/101‑1 (MQ=255)

                                                                                             aGCGCAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTaa     <  1:1294706/101‑1 (MQ=255)

                                                                                              gcgcAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAc    <  2:1294827/101‑1 (MQ=255)

                                                                                               cgcAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc   <  1:1294391/101‑1 (MQ=255)

                                                                                                gcAG‑GC‑CCCCC‑CTCGCCACCGCGCGATGCAACGGTGGTATCGCCGtt                                                        >  1:1295044/1‑47 (MQ=38)

TACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAG‑G‑‑CCCCC‑CTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646702‑646897

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