BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,494,883:1	(A)_7→8	100%	coding (261/564 nt)	idnK →	D‑gluconate kinase, thermosensitive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,494,876	1	.	A	92.0%	405.2 / 26.3	125	coding (254/564 nt)	idnK	D‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand): ref base . (3/7); new base A (56/59); total (59/66)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  NC_000913/4494778‑4494973

aTCTGATGAAGATCAACTTCCACGCCTGGAACGCTTAAATGATGCTTCATCCAGTCTTTATAACAAGCCTGAAACAGGATTTACTGTCTGTTCATCAT‑T‑‑a                                                                                                    >  1:2235695/1‑100 (MQ=255)

 tCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑a                                                                                                  <  1:1810587/100‑1 (MQ=255)

  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aa                                                                                                 <  2:3519382/100‑1 (MQ=255)

    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑a                                                                                              <  1:3777422/100‑1 (MQ=255)

    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑a                                                                                              >  1:1626812/1‑100 (MQ=255)

     atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aa                                                                                             >  2:3429657/1‑100 (MQ=255)

      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaa                                                                                            >  2:3256771/1‑100 (MQ=255)

        aaGATCGACTTCCCTGGCTGGAACGGTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaaac                                                                                          <  1:420700/100‑3 (MQ=255)

        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaaac                                                                                          <  2:218691/100‑3 (MQ=255)

        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaaac                                                                                          <  1:1365654/100‑3 (MQ=255)

        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaaac                                                                                          <  2:1143627/100‑3 (MQ=255)

         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaaaca                                                                                         >  1:3228676/1‑97 (MQ=255)

          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑a‑aaa‑aaaacag                                                                                        <  1:1348387/100‑5 (MQ=255)

             cGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTAt                                                                                      >  2:3146458/1‑100 (MQ=255)

              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATc                                                                                     >  2:2604980/1‑100 (MQ=255)

                     ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGatat                                                                              >  2:3698773/1‑100 (MQ=255)

                        gttggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATAtttt                                                                           <  1:1961133/98‑1 (MQ=255)

                         ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTa                                                                          <  1:868614/100‑1 (MQ=255)

                          tggAACGGTTAAATGATGCTTCATACAGTTTTTAGAAAAAGAAAGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑A‑AAAGAAACAGTATCGTGATATTTTAc                                                                         <  1:2207545/100‑1 (MQ=255)

                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTAc                                                                         <  2:1344727/100‑1 (MQ=255)

                           ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACg                                                                        <  2:885375/100‑1 (MQ=255)

                              aatcTTAAATGATTCTTTATACATTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTCTTCATCAT‑TAAA‑AAA‑AAACAGTATCGTGATATTTTa                                                                          <  1:3003788/93‑1 (MQ=255)

                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAg                                                                    <  2:46856/100‑1 (MQ=255)

                                gattAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAgg                                                                   <  1:2309188/98‑1 (MQ=255)

                                gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAgg                                                                   >  2:1634706/1‑100 (MQ=255)

                                gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAgg                                                                   >  1:1737074/1‑100 (MQ=255)

                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAggg                                                                  <  2:2047833/100‑1 (MQ=255)

                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAggg                                                                  >  2:802037/1‑100 (MQ=255)

                                 cTTAAATGATGATTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAg                                                                    <  1:1541937/98‑1 (MQ=255)

                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGt                                                                 <  2:3647946/100‑1 (MQ=255)

                                   tAAATGAAGCTTCATACAGTTTTTATAAAAAGAATGAAACAGGTTTTATTGTCTGTTAATCAT‑TA‑A‑AAA‑AAACAGTATCATGATGTTTTTCGTAAtgg                                                                  >  1:714016/1‑98 (MQ=255)

                                    aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAg                                                               >  1:1081190/1‑100 (MQ=255)

                                     aaggatgCTTCATACAGTCTTTATAAAAAGAATGAAAAAATATTTATTGTCTGTTTATCAA‑TAAA‑AAA‑AAACAGTATCGTGATATTTTACGCAAggg                                                                  <  1:3373086/94‑1 (MQ=255)

                                       tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            >  2:1368842/1‑100 (MQ=255)

                                       tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                             >  1:103359/1‑99 (MQ=255)

                                         attcattaTACAGTCTTTATAAAAAGAATGAAAAAGGATTTATTGTTTGTTCATAAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGcn                                                             <  1:3463805/90‑2 (MQ=255)

                                         atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                          >  1:3498450/1‑100 (MQ=255)

                                          tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                         <  2:3081290/100‑1 (MQ=255)

                                          tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                         <  2:2588745/100‑1 (MQ=255)

                                           gCTTCATACAGTTTTTATAAAAAGAATGAAACAGGATTTATTGTATGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        <  1:907702/100‑1 (MQ=255)

                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        <  2:3077545/100‑1 (MQ=255)

                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  1:1426428/1‑100 (MQ=255)

                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  1:2641517/1‑100 (MQ=255)

                                            cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAAGAGTATCGTGATATTTTACGCAAGGGTAGCNCCCa                                                         >  2:2870723/1‑98 (MQ=255)

                                            cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                       >  2:3231045/1‑100 (MQ=255)

                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGctccc                                                          >  2:2766463/1‑96 (MQ=255)

                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                       >  2:2091972/1‑99 (MQ=255)

                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                    <  2:2969428/100‑1 (MQ=255)

                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑ACAAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                    >  1:1288348/1‑100 (MQ=255)

                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                    >  2:2288051/1‑99 (MQ=255)

                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   >  2:3255680/1‑100 (MQ=255)

                                                aTACAGAATATATAAAAAAAATGAAACAGGATTTATTGTATGAAAAAAAAAA‑‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     <  1:760197/98‑1 (MQ=255)

                                                  aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                 <  1:2434101/100‑1 (MQ=255)

                                                   cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTAAGCAAGGGTAGCCCCCATGTTCAt                                                  >  1:2959470/1‑98 (MQ=255)

                                                    aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                                <  1:588221/99‑1 (MQ=255)

                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                              >  1:3214267/1‑100 (MQ=255)

                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAAT‑TA‑A‑AAA‑AAACAGTATAGTGATATTTTACGCAAGGGTAACACCCATGTTCAttt                                                >  1:844811/1‑98 (MQ=255)

                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                            <  2:1583983/100‑1 (MQ=255)

                                                         ttATAAAAAGAATGATACAGGATTTATAATCTGTGCAACAT‑T‑‑A‑AAACAAACAGTATCGTGATATTTTACGTAAGGGTAGCCCCCATGTTCATTTCCTCTg                                          <  1:2176148/100‑1 (MQ=255)

                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                          <  1:2334774/100‑1 (MQ=255)

                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                         <  1:1359825/100‑1 (MQ=255)

                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                         <  2:1126633/100‑1 (MQ=255)

                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                         <  1:1181369/100‑1 (MQ=255)

                                                          tataAAAAGAAAGAAACAGGATTTATTGTCTGTACATCATAA‑‑A‑AAA‑AAACAGTAGCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctcn                                           <  1:3028128/98‑2 (MQ=255)

                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        >  1:335734/1‑100 (MQ=255)

                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        <  2:3172634/100‑1 (MQ=255)

                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        >  2:490651/1‑100 (MQ=255)

                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        <  2:1180778/100‑1 (MQ=255)

                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       >  1:1274802/1‑100 (MQ=255)

                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCAATGTTCATTTCCTCTGGtt                                       >  2:2261859/1‑100 (MQ=255)

                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      >  2:2869164/1‑100 (MQ=255)

                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCATCTGGTTa                                      >  1:809437/1‑100 (MQ=255)

                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCACGGGTAGCCCCAAAGTTCATTTCCTCTGGTTAg                                     >  1:2013632/1‑100 (MQ=255)

                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      >  1:2614730/1‑99 (MQ=255)

                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCATGATATTTTACTCAACGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      >  1:2126331/1‑99 (MQ=255)

                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑A‑AAACAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATGTCCTCTGTTTa                                      >  1:1660575/1‑99 (MQ=255)

                                                               aaaGAATGAAACAGGATTAATTGTCTGTTGATCAT‑TA‑A‑AAA‑AAACATTATCATTATGTTTTACGTTACTGTACACACCATTTTCATTTCATCTGGTTAGa                                    >  1:856108/1‑100 (MQ=255)

                                                                 aGAATGAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                                  <  1:3741226/100‑1 (MQ=255)

                                                                      gAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                             >  1:3235203/1‑100 (MQ=255)

                                                                      gAAACAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                             <  2:3213565/100‑1 (MQ=255)

                                                                       aaaCAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                            >  2:777835/1‑100 (MQ=255)

                                                                        aaCAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                           >  2:1118805/1‑100 (MQ=255)

                                                                          cAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         <  2:884441/100‑1 (MQ=255)

                                                                          cAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         <  1:2870321/100‑1 (MQ=255)

                                                                          cAGGATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         >  2:2551276/1‑100 (MQ=255)

                                                                            ggATTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                       >  1:1228755/1‑100 (MQ=255)

                                                                              aTTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  2:103857/100‑1 (MQ=255)

                                                                              aTTTATTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  2:3396820/100‑1 (MQ=255)

                                                                                  aTTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt                 <  2:2004672/100‑1 (MQ=255)

                                                                                  aTTGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt                 >  1:2922435/1‑100 (MQ=255)

                                                                                    tGTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct               <  2:2074098/100‑1 (MQ=255)

                                                                                     gTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGACACTATtctc              >  1:2400009/1‑100 (MQ=255)

                                                                                     gTCTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc              >  2:1286121/1‑100 (MQ=255)

                                                                                       cTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc            <  1:3423203/100‑1 (MQ=255)

                                                                                       cTGTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc            >  1:3187178/1‑100 (MQ=255)

                                                                                         gTTCATCAT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc          <  2:1626359/100‑1 (MQ=255)

                                                                                           tcatcaT‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa        >  1:87921/1‑100 (MQ=255)

                                                                                           taatcat‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg         <  1:1113021/97‑1 (MQ=255)

                                                                                              tcat‑TA‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg     <  2:1690807/100‑1 (MQ=255)

                                                                                                 aat‑‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAg  <  1:2158995/98‑1 (MQ=255)

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAT‑T‑‑A‑AAA‑AAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  NC_000913/4494778‑4494973

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