BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,158,977	T→C	V111A (GTC→GCC)	yijC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,158,977	0	T	C	99.4%	979.9 / ‑4.1	319	V111A (GTC→GCC)	yijC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (157/160); total (157/162)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 331 total aligned reads displayed.

TGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTG  >  NC_000913/4158887‑4159066

   ttGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGAGGGATCCGCACCTCGGCCTCca                                                                                      >  2:6066762/1‑93 (MQ=255)

    tgatgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCcac                                                                                     <  1:6700669/93‑1 (MQ=255)

    tgatgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCcac                                                                                     <  2:1551364/93‑1 (MQ=255)

      atgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAt                                                                                   >  1:1155596/1‑93 (MQ=255)

      atgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAt                                                                                   >  1:3690203/1‑93 (MQ=255)

       tgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACAtt                                                                                  >  2:2515734/1‑93 (MQ=255)

         agagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTa                                                                                <  1:3736729/93‑1 (MQ=255)

           agCGTTTTAATGCTACGCCACCTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATg                                                                              <  1:5352400/93‑1 (MQ=255)

           agCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATg                                                                              <  1:264064/93‑1 (MQ=255)

            gCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATgg                                                                             <  2:2226275/93‑1 (MQ=255)

            gCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATgg                                                                             >  1:3316255/1‑93 (MQ=255)

             cGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGa                                                                            <  2:2044567/93‑1 (MQ=255)

               gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGt                                                                          >  1:5049875/1‑93 (MQ=255)

               gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGt                                                                          <  2:3935337/93‑1 (MQ=255)

               gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGt                                                                          <  1:3467086/93‑1 (MQ=255)

                  tAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCa                                                                       <  1:1254725/93‑1 (MQ=255)

                  tAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCa                                                                       <  2:4115995/93‑1 (MQ=255)

                  tAATGCTACGCCAACTCATGCCCCAGGCGCGTCAGCTTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGACCTCCACATTTATGGAGTTCa                                                                       <  1:314734/93‑1 (MQ=255)

                   aaTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCAt                                                                      >  1:6236393/1‑93 (MQ=255)

                    aTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATc                                                                     <  1:3413556/93‑1 (MQ=255)

                     tGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTGCACATTTATGGAGGTCATCg                                                                    >  1:250017/1‑93 (MQ=255)

                     agCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCg                                                                    <  1:6796824/92‑1 (MQ=255)

                          cGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaat                                                               <  2:7210576/93‑1 (MQ=255)

                            ccAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaataa                                                             <  1:7628017/93‑1 (MQ=255)

                            ccAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGtaataa                                                             >  2:7129113/1‑93 (MQ=255)

                             cAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTaataat                                                            <  2:5470399/93‑1 (MQ=255)

                             cAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTaataat                                                            <  2:5887386/93‑1 (MQ=255)

                              aaCTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATc                                                           >  1:4384227/1‑93 (MQ=255)

                                cTCATGCGCCAGGCGCGTCAGGGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCt                                                         >  2:5012950/1‑93 (MQ=255)

                                cTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCt                                                         <  1:5714438/93‑1 (MQ=255)

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                       <  2:7054002/93‑1 (MQ=255)

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                       <  1:752664/93‑1 (MQ=255)

                                  cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTaa                                                       <  1:7406394/93‑1 (MQ=255)

                                   aTGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAAc                                                      <  1:7362610/93‑1 (MQ=255)

                                       gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCt                                                  >  2:100990/1‑93 (MQ=255)

                                       gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCt                                                  >  1:1870959/1‑93 (MQ=255)

                                       gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCt                                                  <  2:3569053/93‑1 (MQ=255)

                                        ccAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCtt                                                 >  2:2422786/1‑93 (MQ=255)

                                        ccAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCtt                                                 >  2:5172873/1‑93 (MQ=255)

                                        ccAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCtt                                                 <  2:7497307/93‑1 (MQ=255)

                                         cAGGCGCGTCAGCGTATCGCCAAGGGCGGGAGTGTGATCCGCGCCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                                <  2:858212/93‑1 (MQ=255)

                                         cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                                <  2:2804486/93‑1 (MQ=255)

                                          aGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTcc                                               <  2:3069578/93‑1 (MQ=255)

                                             cgcgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGt                                            <  1:2882889/93‑1 (MQ=255)

                                               cgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGtta                                          <  1:6461433/93‑1 (MQ=255)

                                               cgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGtta                                          <  1:1488616/93‑1 (MQ=255)

                                                gTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttat                                         >  1:4084482/1‑93 (MQ=255)

                                                 tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                        >  1:2010429/1‑93 (MQ=255)

                                                 tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                        >  1:2076056/1‑93 (MQ=255)

                                                 tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttatt                                        >  2:446046/1‑93 (MQ=255)

                                                   aGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                      <  1:5149461/93‑1 (MQ=255)

                                                   aGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                      >  2:1455803/1‑93 (MQ=255)

                                                   aGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                      <  2:5271883/93‑1 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                     <  2:5530515/93‑1 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                     <  2:3244318/93‑1 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                     >  2:563062/1‑93 (MQ=255)

                                                    gCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                     <  2:2839930/93‑1 (MQ=255)

                                                     cGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTg                                    <  1:907698/93‑1 (MQ=255)

                                                     cGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTg                                    <  2:7803817/93‑1 (MQ=255)

                                                      gTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGc                                   <  1:5658558/93‑1 (MQ=255)

                                                      gTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGc                                   <  1:6031654/93‑1 (MQ=255)

                                                      gTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGc                                   <  1:6866458/93‑1 (MQ=255)

                                                       tATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCg                                  <  2:667014/93‑1 (MQ=255)

                                                       tATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCg                                  >  1:6565022/1‑93 (MQ=255)

                                                       tATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCg                                  <  1:195987/93‑1 (MQ=255)

                                                       tATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACACATTACGGTTATTATTGCg                                  >  2:3508709/1‑93 (MQ=255)

                                                        aTCGCCAAAGTCGGGAGTGTGATCCGCACCTCGGCCTCCCCATTTATGGAGTTCAGCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCgg                                 <  2:3814570/93‑1 (MQ=255)

                                                        aTCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCgg                                 <  1:1754766/93‑1 (MQ=255)

                                                         tCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCggg                                >  2:4685124/1‑93 (MQ=255)

                                                          cGCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGa                               <  1:2988404/93‑1 (MQ=255)

                                                           gCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGaa                              <  1:2415639/93‑1 (MQ=255)

                                                           gCCAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGaa                              >  1:5899599/1‑93 (MQ=255)

                                                            ccAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAAc                             <  1:4471030/93‑1 (MQ=255)

                                                            ccAAAGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAAc                             <  1:7008383/93‑1 (MQ=255)

                                                             cAAAGGCGGGAGGGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACg                            >  2:1150115/1‑93 (MQ=255)

                                                             cAAAGGCGCGAGTGTGATCCACACCTCGGCCTGTACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTGTTATTGCGGGAACg                            <  2:6927339/93‑1 (MQ=255)

                                                              aaaGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGc                           <  2:5697156/93‑1 (MQ=255)

                                                              aaaGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGc                           >  2:7276443/1‑93 (MQ=255)

                                                                aGGCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTc                         <  1:55883/93‑1 (MQ=255)

                                                                 ggCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTcc                        <  2:945296/93‑1 (MQ=255)

                                                                 ggCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTcc                        <  2:1768246/93‑1 (MQ=255)

                                                                  gCGGGAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCg                       <  1:471624/93‑1 (MQ=255)

                                                                    gggAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGc                     >  2:4332280/1‑93 (MQ=255)

                                                                      gAGTGTGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCAc                   <  2:4013224/93‑1 (MQ=255)

                                                                        gtgtGATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCt                 >  2:6137658/1‑93 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc                >  1:3589270/1‑93 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc                >  2:1217522/1‑93 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc                <  1:5993043/93‑1 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc                >  1:4808338/1‑93 (MQ=255)

                                                                         tgtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGAGGGAACGCTCCGGCACCTc                >  2:6321888/1‑93 (MQ=255)

                                                                         ggtgATCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc                <  1:4186153/92‑1 (MQ=255)

                                                                              tCCGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTg           >  2:6737797/1‑93 (MQ=255)

                                                                                cGCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCg         <  1:829621/93‑1 (MQ=255)

                                                                                 gCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGt        <  1:1441207/93‑1 (MQ=255)

                                                                                 gCACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGt        >  1:4171603/1‑93 (MQ=255)

                                                                                  cACCTCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGtt       >  1:2328664/1‑93 (MQ=255)

                                                                                    ccTCGACCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTc     >  2:4867570/1‑93 (MQ=255)

                                                                                      tCGGCCTCCACATTTATGGAGTTCATCGGTAATAGTCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGt   >  2:4955029/1‑93 (MQ=255)

                                                                                      tCGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGt   >  2:3172560/1‑93 (MQ=255)

                                                                                       cGGCCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTg  >  1:2204619/1‑93 (MQ=255)

TGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTG  >  NC_000913/4158887‑4159066

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