BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,673,594	T→C	Y74C (TAC→TGC)	yphC ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,673,594	0	T	C	98.5%	422.4 / ‑4.1	132	Y74C (TAC→TGC)	yphC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (76/54); total (76/56)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 143 total aligned reads displayed.

GCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATT  >  NC_000913/2673505‑2673683

gCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGcata                                                                                        >  1:4358738/1‑93 (MQ=255)

gCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGcata                                                                                        <  1:2903342/93‑1 (MQ=255)

gCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGcata                                                                                        >  2:1193771/1‑93 (MQ=255)

  aCGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAAc                                                                                      >  2:1314560/1‑93 (MQ=255)

  aCGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAAc                                                                                      <  1:730868/93‑1 (MQ=255)

   cgcgGTCGCCCTCTTTAAAATGGAGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCACTGGCATAACg                                                                                     >  2:4596960/1‑93 (MQ=255)

    gcgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACgg                                                                                    >  2:4479956/1‑93 (MQ=255)

      ggTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGtt                                                                                  >  2:1106591/1‑93 (MQ=255)

       gTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGttt                                                                                 >  1:2679468/1‑93 (MQ=255)

       gTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGttt                                                                                 >  2:3433204/1‑93 (MQ=255)

        tCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTa                                                                                >  1:1148436/1‑93 (MQ=255)

        tCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTa                                                                                >  1:2438190/1‑93 (MQ=255)

        tCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTa                                                                                <  1:2388162/93‑1 (MQ=255)

        tCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTa                                                                                <  1:2333100/93‑1 (MQ=255)

          gCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATc                                                                              <  2:1412767/93‑1 (MQ=255)

          gCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATc                                                                              >  1:1871163/1‑93 (MQ=255)

          gCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATc                                                                              <  2:4341383/93‑1 (MQ=255)

           cccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATACCGGTTTATCg                                                                             >  2:1084781/1‑93 (MQ=255)

             ctctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCggg                                                                           >  2:1242956/1‑93 (MQ=255)

              tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          >  1:1037924/1‑93 (MQ=255)

              tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          <  1:3489635/93‑1 (MQ=255)

              tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          >  2:2844792/1‑93 (MQ=255)

               ctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTg                                                                         <  1:1430618/93‑1 (MQ=255)

                tttAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgc                                                                        <  1:3495579/93‑1 (MQ=255)

                    aaaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgccgcc                                                                    <  1:2288889/93‑1 (MQ=255)

                      aTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCcgccgc                                                                  <  2:1934523/93‑1 (MQ=255)

                      aTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCcgccgc                                                                  >  2:908870/1‑93 (MQ=255)

                        ggcggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCtg                                                                <  2:2970037/93‑1 (MQ=255)

                         gcggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACTGTTTATCGGGTGCCGCCGCtgt                                                               <  2:2989733/93‑1 (MQ=255)

                         gcggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCtgt                                                               <  1:597402/93‑1 (MQ=255)

                         gcggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCtgt                                                               <  2:4453991/93‑1 (MQ=255)

                           ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTgggg                                                             >  1:4255987/1‑93 (MQ=255)

                            gcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGc                                                            >  2:4116613/1‑93 (MQ=255)

                              atccTTACCCCATCGCCACAATCTGCCTGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTCCCGCCGCTGTGGctc                                                          <  2:922002/91‑3 (MQ=255)

                              aGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCAc                                                          >  2:3043243/1‑93 (MQ=255)

                              aGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCAc                                                          >  1:580712/1‑93 (MQ=255)

                               gCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACg                                                         >  1:91986/1‑93 (MQ=255)

                               gCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACg                                                         >  2:407261/1‑93 (MQ=255)

                               gCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACg                                                         >  2:1206098/1‑93 (MQ=255)

                                ccTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACgg                                                        <  1:426985/93‑1 (MQ=255)

                                 cTTGCCCCATCGCCACAATCTGCACGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgt                                                       <  2:2598352/93‑1 (MQ=255)

                                  ttGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtg                                                      >  2:1560842/1‑93 (MQ=255)

                                  ttGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtg                                                      <  1:4707189/93‑1 (MQ=255)

                                   tGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtgt                                                     <  2:1877692/93‑1 (MQ=255)

                                   tGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtgt                                                     >  1:2206720/1‑93 (MQ=255)

                                    tccccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGtt                                                    <  1:4023633/92‑1 (MQ=255)

                                    gCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGtt                                                    >  1:2734085/1‑93 (MQ=255)

                                      cccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtga                                                  >  2:3456683/1‑93 (MQ=255)

                                       ccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCGGCCGCTGTGGCACGGTGTtgat                                                 <  1:2237204/93‑1 (MQ=255)

                                        cATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatg                                                >  1:2100025/1‑93 (MQ=255)

                                         aTCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatga                                               <  1:3958427/93‑1 (MQ=255)

                                          tCGCCACAATCTGCCCGCACGTTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTgatgat                                              >  1:846044/1‑93 (MQ=255)

                                          tCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTgatgat                                              >  1:846045/1‑93 (MQ=255)

                                          tCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTgatgat                                              >  1:1251679/1‑93 (MQ=255)

                                           cGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgata                                             <  1:3780785/93‑1 (MQ=255)

                                            gCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAATGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatag                                            >  1:3940502/1‑93 (MQ=255)

                                             ccACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgataga                                           <  2:662891/93‑1 (MQ=255)

                                              cacaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagat                                          >  1:3031847/1‑93 (MQ=255)

                                              cacaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagat                                          >  2:2146402/1‑93 (MQ=255)

                                               acaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagata                                         >  2:4695354/1‑93 (MQ=255)

                                               acaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGCGCCGCCGCTGTGGCACGGTGTTGATgatagata                                         >  1:1775232/1‑93 (MQ=255)

                                                 aaTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGAtata                                       <  1:1623375/93‑1 (MQ=255)

                                                 aaTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGAtata                                       <  1:3470762/93‑1 (MQ=255)

                                                  aTCTGCCCGCACTGTTCATGACCGTTGATAAAGCCCGGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAg                                      <  2:2763762/93‑1 (MQ=255)

                                                     tGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTgg                                   <  1:4664798/93‑1 (MQ=255)

                                                      tcccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGa                                  <  2:4285282/92‑1 (MQ=255)

                                                       cccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGAc                                 <  1:761249/93‑1 (MQ=255)

                                                        ccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACa                                <  2:4441758/93‑1 (MQ=255)

                                                        ccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACa                                >  1:1231924/1‑93 (MQ=255)

                                                          gcccGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTCCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATc                              <  1:3212614/90‑1 (MQ=255)

                                                          gCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATc                              <  1:3656823/93‑1 (MQ=255)

                                                          gCACGGTTCATGACCGTTGATAAAGCCCTGGCAGAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATc                              <  1:4569522/93‑1 (MQ=255)

                                                            aCGTTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGc                            <  2:1279361/93‑1 (MQ=255)

                                                              ggTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCtt                          >  2:2611980/1‑93 (MQ=255)

                                                               gTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCAAGGTGTTGATGATAGATATAGTGGAAATCGCTTc                         >  1:4039754/1‑93 (MQ=255)

                                                                ttCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTcc                        >  1:1353059/1‑93 (MQ=255)

                                                                   aTGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGGGTTGATGATAGATATAGTGGACATCGCTTCCGCa                     >  1:4041237/1‑93 (MQ=255)

                                                                    tGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCaa                    >  1:379218/1‑93 (MQ=255)

                                                                    tGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCaa                    >  2:1704155/1‑93 (MQ=255)

                                                                     gACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCaaa                   <  2:3367810/93‑1 (MQ=255)

                                                                       ccGTTGATAAAGCCCTGGCATAACGGTTTATCGGGGGCCGCCGCTGTGGCACGGGGTGGAGGATGGATATAGTGGACATGGCTTCGGCAAATc                 >  1:2425230/1‑93 (MQ=255)

                                                                        cGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATcc                >  2:2199226/1‑93 (MQ=255)

                                                                         gtggATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATccc               <  2:2871011/90‑1 (MQ=255)

                                                                         gTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATccc               <  2:1795981/93‑1 (MQ=255)

                                                                          ttGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCg              <  1:1134270/93‑1 (MQ=255)

                                                                           tGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgg             >  2:3279969/1‑93 (MQ=255)

                                                                            gATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgga            <  1:975735/93‑1 (MQ=255)

                                                                              tAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagg          <  1:3518724/93‑1 (MQ=255)

                                                                              tAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagg          >  1:2394280/1‑93 (MQ=255)

                                                                              tAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagg          >  1:906998/1‑93 (MQ=255)

                                                                               aaaGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagga         <  1:1588873/93‑1 (MQ=255)

                                                                               aaaGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagga         >  2:3824963/1‑93 (MQ=255)

                                                                                 aGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAtt       <  2:2646690/93‑1 (MQ=255)

                                                                                 aGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAtt       >  2:546051/1‑93 (MQ=255)

                                                                                 aGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAtt       <  1:4558493/93‑1 (MQ=255)

                                                                                   cccTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTc     <  2:4009845/93‑1 (MQ=255)

                                                                                   cccTGGCATAACGCTTTATGGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAACCCCGGAGGATTTc     <  2:2555608/93‑1 (MQ=255)

                                                                                    ccTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCa    <  1:1583752/93‑1 (MQ=255)

                                                                                     cTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCAt   >  2:2684603/1‑93 (MQ=255)

                                                                                      tGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCAtt  >  2:2323545/1‑93 (MQ=255)

GCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATT  >  NC_000913/2673505‑2673683

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