BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,110,991	A→G	D40G (GAT→GGT)	yiiT →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,110,991	0	A	G	84.4%	729.9 / 127.8	315	D40G (GAT→GGT)	yiiT	putative regulator
Reads supporting (aligned to +/- strand): ref base A (20/29); new base G (143/123); total (163/152)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 338 total aligned reads displayed.

AATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTG  >  NC_000913/4110899‑4111082

aaTTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATtga                                                                                             <  1:6266845/93‑1 (MQ=255)

  tttCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATtggtg                                                                                           >  1:4430217/1‑93 (MQ=255)

  tttCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATtggtg                                                                                           >  2:3739754/1‑93 (MQ=255)

   ttCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATtggtga                                                                                          <  1:655147/93‑1 (MQ=255)

   ttCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATtggtga                                                                                          <  2:4440081/93‑1 (MQ=255)

    tCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTggtgat                                                                                         >  2:6462969/1‑93 (MQ=255)

    tCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTggtgat                                                                                         <  1:2289259/93‑1 (MQ=255)

      cGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTCACGTTAATTCATATTGGTGATgg                                                                                       >  1:5034771/1‑93 (MQ=255)

       gggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGc                                                                                      <  1:1214774/93‑1 (MQ=255)

         gAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCtt                                                                                    >  1:6650360/1‑93 (MQ=255)

           aTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTaa                                                                                  <  2:2147182/93‑1 (MQ=255)

           aTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTaa                                                                                  <  2:436214/93‑1 (MQ=255)

            tGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGTCTTAAg                                                                                 <  2:6990279/93‑1 (MQ=255)

            tGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAg                                                                                 <  1:5579635/93‑1 (MQ=255)

            tGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAg                                                                                 >  1:930521/1‑93 (MQ=255)

             gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGc                                                                                >  1:1990046/1‑93 (MQ=255)

             gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGc                                                                                >  2:2540511/1‑93 (MQ=255)

             gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGc                                                                                >  2:405003/1‑93 (MQ=255)

               agaagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGa                                                                              <  2:4956450/93‑1 (MQ=255)

               agaagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGa                                                                              >  2:5580618/1‑93 (MQ=255)

                gaagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAg                                                                             >  1:4692836/1‑93 (MQ=255)

                 aagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGt                                                                            <  1:6932190/93‑1 (MQ=255)

                   gaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTg                                                                          >  2:1624555/1‑93 (MQ=255)

                     tGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTa                                                                        >  2:6200565/1‑93 (MQ=255)

                     tGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTa                                                                        <  2:7316482/93‑1 (MQ=255)

                      gCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTAt                                                                       <  1:4981565/93‑1 (MQ=255)

                        cTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATcc                                                                     <  1:4300459/93‑1 (MQ=255)

                         ttACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCg                                                                    >  2:5500970/1‑93 (MQ=255)

                         ttACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCg                                                                    >  2:2858099/1‑93 (MQ=255)

                         ttACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCg                                                                    >  2:4289842/1‑93 (MQ=255)

                          tACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCgg                                                                   <  2:1408013/93‑1 (MQ=255)

                            cTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGt                                                                 <  2:3658694/93‑1 (MQ=255)

                             tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTa                                                                >  2:2422015/1‑93 (MQ=255)

                             tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTa                                                                >  2:3891112/1‑93 (MQ=255)

                             tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCAAATTGGTGATGGCTTAAGCGAGTTGTATCCgtgca                                                                >  1:1607129/1‑89 (MQ=255)

                               gTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATc                                                              <  1:7464069/93‑1 (MQ=255)

                               gTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATc                                                              >  1:3605445/1‑93 (MQ=255)

                                 gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTa                                                            >  2:835919/1‑93 (MQ=255)

                                  aaTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTAc                                                           <  2:2248322/93‑1 (MQ=255)

                                   aTAAAGCCCTGGAGCTCACCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACt                                                          <  1:4149657/93‑1 (MQ=255)

                                    tAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACtt                                                         >  2:6273665/1‑93 (MQ=255)

                                    tAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACtt                                                         >  1:2454765/1‑93 (MQ=255)

                                    tAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACtt                                                         >  1:4885033/1‑93 (MQ=255)

                                    tAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACtt                                                         >  1:575638/1‑93 (MQ=255)

                                    tAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACtt                                                         >  1:455987/1‑93 (MQ=255)

                                     aaaGCCCTGGAGCTCGCCAGACATAATCACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTc                                                        >  2:5140186/1‑93 (MQ=255)

                                      aaGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTcc                                                       >  2:2652123/1‑93 (MQ=255)

                                       aGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTccc                                                      <  1:2930775/93‑1 (MQ=255)

                                       aGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTccc                                                      <  1:6320696/93‑1 (MQ=255)

                                            tGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGTGTTGTATCCGGGTATCTACTTCCCTGCaa                                                 >  1:7083757/1‑93 (MQ=255)

                                              gAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACa                                               >  1:774164/1‑93 (MQ=255)

                                              gAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACa                                               <  2:174877/93‑1 (MQ=255)

                                                gCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACaga                                             >  1:6807996/1‑93 (MQ=255)

                                                 cTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACagaa                                            >  2:5086867/1‑93 (MQ=255)

                                                 cTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACagaa                                            <  1:3829510/93‑1 (MQ=255)

                                                    gCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGat                                         >  1:5633776/1‑93 (MQ=255)

                                                     ccAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGata                                        <  2:4101637/93‑1 (MQ=255)

                                                      cAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGGGTTGTATCCGGGAATCTACTTCCCTGCAACAGAAGatat                                       >  1:1838296/1‑93 (MQ=255)

                                                      cAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGatat                                       >  2:2005172/1‑93 (MQ=255)

                                                      cAGACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGatat                                       <  2:5434987/93‑1 (MQ=255)

                                                        gACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAttc                                     >  2:5891647/1‑93 (MQ=255)

                                                        gACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAttc                                     >  2:1011354/1‑93 (MQ=255)

                                                        gACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAttc                                     <  2:1789394/93‑1 (MQ=255)

                                                        gACATAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAtt                                      <  2:2976519/92‑1 (MQ=255)

                                                           aTAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAttcttc                                  <  1:1313478/93‑1 (MQ=255)

                                                           aTAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAttcttc                                  <  1:2371593/93‑1 (MQ=255)

                                                           aTAATGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATAttcttc                                  >  1:3943662/1‑93 (MQ=255)

                                                            tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCCACAGAAGATATTCTTCa                                 <  2:4904086/93‑1 (MQ=255)

                                                               tGACGCTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                              >  2:419348/1‑93 (MQ=255)

                                                                   gCTCACCTGACGTTTATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                          >  1:6081881/1‑93 (MQ=255)

                                                                    cTCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                         <  2:6699141/93‑1 (MQ=255)

                                                                    cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                         <  1:2045259/93‑1 (MQ=255)

                                                                     tCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                        >  2:7431585/1‑93 (MQ=255)

                                                                     tCACCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                        >  1:7267141/1‑93 (MQ=255)

                                                                     tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                        <  2:3334178/93‑1 (MQ=255)

                                                                       aCCTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaaga                      >  2:6845359/1‑93 (MQ=255)

                                                                       aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaaga                      <  2:7340502/93‑1 (MQ=255)

                                                                        ccTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaagaa                     <  2:6405133/93‑1 (MQ=255)

                                                                        ccTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaagaa                     >  1:3015208/1‑93 (MQ=255)

                                                                         cTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                    >  2:3491351/1‑93 (MQ=255)

                                                                         cTGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                    >  1:2606862/1‑93 (MQ=255)

                                                                         cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                    >  1:6453494/1‑93 (MQ=255)

                                                                          tGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                   <  1:1004859/93‑1 (MQ=255)

                                                                          tGACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                   <  2:5267434/93‑1 (MQ=255)

                                                                           gACGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATaa                  <  1:1696786/93‑1 (MQ=255)

                                                                            aCGTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                 <  2:6093488/93‑1 (MQ=255)

                                                                             cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                <  1:3483382/93‑1 (MQ=255)

                                                                             cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                <  1:523086/93‑1 (MQ=255)

                                                                              gTTAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc               <  2:7516979/93‑1 (MQ=255)

                                                                              gTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc               <  2:3157827/93‑1 (MQ=255)

                                                                               ttAATTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg              >  1:7382295/1‑93 (MQ=255)

                                                                                  aTTCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt           >  1:6193288/1‑93 (MQ=255)

                                                                                   ttCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa          >  1:4444265/1‑93 (MQ=255)

                                                                                    tCATATTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa         <  2:1961618/93‑1 (MQ=255)

                                                                                      atatTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa       <  2:502793/93‑1 (MQ=255)

                                                                                      atatTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa       >  2:1536174/1‑93 (MQ=255)

                                                                                       tatTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACaa      <  1:5251309/93‑1 (MQ=255)

                                                                                       tatTGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACaa      >  2:4741497/1‑93 (MQ=255)

                                                                                         ttGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc    <  2:14474/93‑1 (MQ=255)

                                                                                         ttGGTGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc    >  2:6001000/1‑93 (MQ=255)

AATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTG  >  NC_000913/4110899‑4111082

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