BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,595,883	A→G	T76T (ACT→ACC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,595,883	0	A	G	90.0%	654.2 / 54.9	221	T76T (ACT→ACC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (11/11); new base G (93/106); total (104/117)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 223 total aligned reads displayed.

AATAGTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCT  >  NC_000913/1595793‑1595975

 aTAGTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTa                                                                                           >  1:3519796/1‑93 (MQ=255)

 aTAGTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTa                                                                                           <  2:5424467/93‑1 (MQ=255)

 aTAGTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTa                                                                                           <  2:1653455/93‑1 (MQ=255)

   aGTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATc                                                                                         <  1:5627240/93‑1 (MQ=255)

    gTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCa                                                                                        >  2:3946786/1‑93 (MQ=255)

     tGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCaa                                                                                       <  2:3129911/93‑1 (MQ=255)

     tGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCaa                                                                                       <  2:2850038/93‑1 (MQ=255)

     tGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCaa                                                                                       <  1:3137192/93‑1 (MQ=255)

      ggTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAat                                                                                      >  1:4083766/1‑93 (MQ=255)

       gTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAata                                                                                     >  1:6553017/1‑93 (MQ=255)

       gTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAata                                                                                     <  1:745971/93‑1 (MQ=255)

       gTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAata                                                                                     >  1:1289162/1‑93 (MQ=255)

        ttAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAatat                                                                                    >  2:1283706/1‑93 (MQ=255)

         tAAGGAGGCATTACGGCCAGCCAAAATATTAAGCACTCCAGTGCCATACTAGCCAACCAGGTAGGCATCATAATCAACATCGGTATCAATAtt                                                                                   <  2:3778879/93‑1 (MQ=255)

         tAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATAtt                                                                                   <  2:4170834/93‑1 (MQ=255)

          aaGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTg                                                                                  <  1:5769606/93‑1 (MQ=255)

           aggaggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTgg                                                                                 >  2:4852365/1‑93 (MQ=255)

           aggaggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTgg                                                                                 >  1:214852/1‑93 (MQ=255)

           aggaggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTgg                                                                                 <  2:6497329/93‑1 (MQ=255)

             gaggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTa                                                                               >  2:1209197/1‑93 (MQ=255)

             gaggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTa                                                                               <  2:5615014/93‑1 (MQ=255)

              aggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTaa                                                                              <  1:163386/93‑1 (MQ=255)

               ggCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAAt                                                                             >  1:1570583/1‑93 (MQ=255)

                 cATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCt                                                                           >  2:2796857/1‑93 (MQ=255)

                 cATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCt                                                                           <  2:3348243/93‑1 (MQ=255)

                  aTTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTg                                                                          >  1:6297349/1‑93 (MQ=255)

                    tACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGAt                                                                        <  1:1784467/93‑1 (MQ=255)

                    tACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGAt                                                                        >  2:1429741/1‑93 (MQ=255)

                     aCCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATc                                                                       <  1:1032186/93‑1 (MQ=255)

                       cGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCAt                                                                     >  1:3154486/1‑93 (MQ=255)

                         ccagccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTa                                                                   <  2:4019837/93‑1 (MQ=255)

                           agccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATc                                                                 >  1:2039061/1‑93 (MQ=255)

                            gccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCa                                                                >  2:376112/1‑93 (MQ=255)

                            gccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCa                                                                <  1:1148747/93‑1 (MQ=255)

                            gccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCa                                                                >  2:2947677/1‑93 (MQ=255)

                              caAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAAc                                                              <  1:5572498/93‑1 (MQ=255)

                              caAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAAc                                                              >  2:4103717/1‑93 (MQ=255)

                               aaaaTATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACt                                                             <  1:5084169/93‑1 (MQ=255)

                                aaaTATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACtt                                                            <  2:5128939/93‑1 (MQ=255)

                                  atatTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTcc                                                          >  1:91388/1‑93 (MQ=255)

                                  atatTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTcc                                                          <  2:467446/93‑1 (MQ=255)

                                   tatTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCa                                                         >  1:2468422/1‑93 (MQ=255)

                                   tatTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCa                                                         >  2:2270779/1‑93 (MQ=255)

                                     ttAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGa                                                       <  2:2547021/93‑1 (MQ=255)

                                     ttAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGa                                                       <  2:1106020/93‑1 (MQ=255)

                                     ttAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGa                                                       >  1:1888653/1‑93 (MQ=255)

                                     ttAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGa                                                       >  1:3724314/1‑93 (MQ=255)

                                      tAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGAc                                                      <  1:4657233/93‑1 (MQ=255)

                                      tAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGAc                                                      >  2:6198163/1‑93 (MQ=255)

                                        aGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACtt                                                    <  2:855167/93‑1 (MQ=255)

                                          cACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGc                                                  <  2:2999828/93‑1 (MQ=255)

                                           ccTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCa                                                 <  1:3372639/92‑1 (MQ=255)

                                           aCTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCa                                                 >  2:3720209/1‑93 (MQ=255)

                                            cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCAc                                                <  1:1937776/93‑1 (MQ=255)

                                             tCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCAcc                                               <  1:778711/93‑1 (MQ=255)

                                                aGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGa                                            <  2:747656/93‑1 (MQ=255)

                                                   gCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGc                                         <  2:4231264/93‑1 (MQ=255)

                                                    ccATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCt                                        >  2:5510934/1‑93 (MQ=255)

                                                    ccATACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCt                                        <  2:3496949/93‑1 (MQ=255)

                                                      aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGa                                      >  1:6148909/1‑93 (MQ=255)

                                                      aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGa                                      >  1:4354411/1‑93 (MQ=255)

                                                       tACCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGaa                                     >  1:4404719/1‑93 (MQ=255)

                                                        aCCAGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAAc                                    <  1:5647499/93‑1 (MQ=255)

                                                         ccAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAAcc                                   <  2:5295499/93‑1 (MQ=255)

                                                           aGCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATAATGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGc                                 <  1:3526601/93‑1 (MQ=255)

                                                            gCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCt                                >  2:4497279/1‑93 (MQ=255)

                                                            gCCAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCt                                >  2:2181785/1‑93 (MQ=255)

                                                              cAACCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTTAACCGCTTa                              <  2:6645900/93‑1 (MQ=255)

                                                               aaCCAGGTAAGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAg                             >  2:4142944/1‑93 (MQ=255)

                                                                aCCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGt                            <  2:4939436/93‑1 (MQ=255)

                                                                aCCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGt                            >  1:2650023/1‑93 (MQ=255)

                                                                aCCAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGCAATCTGATCATTATCAACTTCCAGACTTGCACCAGACTCTGAACCGCTTAGt                            >  2:1129912/1‑93 (MQ=255)

                                                                  cAGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGc                          <  2:6538738/93‑1 (MQ=255)

                                                                   aGGTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGcc                         <  2:1479284/93‑1 (MQ=255)

                                                                    ggTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCCTAGTGCCa                        <  1:6485663/93‑1 (MQ=255)

                                                                     gTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAg                       <  1:3426164/93‑1 (MQ=255)

                                                                     gTAGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAg                       >  1:5398235/1‑93 (MQ=255)

                                                                       aGGCATCATAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCa                     >  1:5335180/1‑93 (MQ=255)

                                                                          catcatAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCaaaa                  >  1:1518440/1‑93 (MQ=255)

                                                                           atcatAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAac                 >  2:2017042/1‑93 (MQ=255)

                                                                           atcatAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAac                 >  1:3879547/1‑93 (MQ=255)

                                                                            tcatAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAaca                <  2:5064971/93‑1 (MQ=255)

                                                                            tcatAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAaca                <  1:5295645/93‑1 (MQ=255)

                                                                            tcatAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAaca                <  1:3954042/93‑1 (MQ=255)

                                                                             catAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAacac               >  2:2759778/1‑93 (MQ=255)

                                                                              atAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACAcc              <  1:1220718/93‑1 (MQ=255)

                                                                               tAAGCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACAccc             <  1:1482718/93‑1 (MQ=255)

                                                                                  gCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc          >  1:4918288/1‑93 (MQ=255)

                                                                                  gCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc          <  2:2394602/93‑1 (MQ=255)

                                                                                  gCAACATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc          >  2:6027577/1‑93 (MQ=255)

                                                                                    aaCATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGt        <  1:411441/93‑1 (MQ=255)

                                                                                    aaCATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGt        <  2:6314218/93‑1 (MQ=255)

                                                                                    aaCATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGt        <  2:5213731/93‑1 (MQ=255)

                                                                                     aCATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTc       >  2:5846272/1‑93 (MQ=255)

                                                                                     aCATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTc       <  1:3715751/93‑1 (MQ=255)

                                                                                      cATCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCa      <  1:2796996/93‑1 (MQ=255)

                                                                                       aTCGGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCaa     >  2:5659951/1‑93 (MQ=255)

                                                                                       aTCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCaa     <  1:3582824/93‑1 (MQ=255)

                                                                                          aGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCt  <  1:1234148/93‑1 (MQ=255)

                                                                                          aGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCt  <  1:2025097/93‑1 (MQ=255)

AATAGTGGTTAAGGAGGCATTACCGCCAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCT  >  NC_000913/1595793‑1595975

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