BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,178	Δ1 bp	coding (899/930 nt)	metA →	homoserine O‑transsuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,178	0	A	.	100.0%	554.3 / NA	120	coding (899/930 nt)	metA	homoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (56/64); total (56/64)

Only 100 of 120 total aligned reads displayed.

TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215083‑4215272

tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                                                           <  2:562305/100‑1 (MQ=255)

tCCGCAAAATACACCGCGAGCGAGCTAGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                                                           >  2:77816/1‑100 (MQ=255)

 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                          <  2:1556414/100‑1 (MQ=255)

  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACgg                                                                                         >  1:1516713/1‑100 (MQ=255)

   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACTAGATCACGCCATACG‑TCTACGGc                                                                                        >  2:253178/1‑100 (MQ=255)

   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                        <  1:2008450/100‑1 (MQ=255)

   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                        >  1:1825067/1‑100 (MQ=255)

    cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGca                                                                                       >  2:170773/1‑100 (MQ=255)

    cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGca                                                                                       <  2:383174/100‑1 (MQ=255)

     aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcac                                                                                      <  1:2587638/100‑1 (MQ=255)

      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcaca                                                                                     <  1:2635733/100‑1 (MQ=255)

      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                        >  2:832205/1‑97 (MQ=255)

       aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACAt                                                                                    <  1:16405/100‑1 (MQ=255)

       aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACAt                                                                                    <  1:1836870/100‑1 (MQ=255)

       aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACAt                                                                                    <  2:2477586/100‑1 (MQ=255)

        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                   <  2:2797721/100‑1 (MQ=255)

        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                   <  2:565054/100‑1 (MQ=255)

        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                   <  2:1757201/100‑1 (MQ=255)

        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                   <  2:6852/100‑1 (MQ=255)

         tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGa                                                                                  <  2:1860886/100‑1 (MQ=255)

           cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAAt                                                                                >  1:1016187/1‑100 (MQ=255)

           cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAAt                                                                                <  2:1986841/100‑1 (MQ=255)

             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              >  2:487762/1‑100 (MQ=255)

             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              >  1:710918/1‑100 (MQ=255)

             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              >  1:1515277/1‑100 (MQ=255)

               gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCaa                                                                            <  2:2012482/100‑1 (MQ=255)

                cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAc                                                                           >  1:1279786/1‑100 (MQ=255)

                cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAc                                                                           <  2:2788955/100‑1 (MQ=255)

                cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAc                                                                           >  2:123562/1‑100 (MQ=255)

                 gagcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACg                                                                          <  2:2410596/100‑1 (MQ=255)

                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                         <  1:2784172/100‑1 (MQ=255)

                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                         <  1:1325836/100‑1 (MQ=255)

                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                         >  1:106111/1‑100 (MQ=255)

                        cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGAtt                                                                   <  2:1173243/100‑1 (MQ=255)

                         tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTa                                                                  >  2:1541723/1‑100 (MQ=255)

                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               <  2:788103/100‑1 (MQ=255)

                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  1:1939996/1‑100 (MQ=255)

                             gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtct                                                              <  1:2229724/100‑1 (MQ=255)

                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                            <  2:794756/100‑1 (MQ=255)

                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                        <  1:762555/100‑1 (MQ=255)

                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                       <  2:1635546/100‑1 (MQ=255)

                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                       >  2:2551799/1‑100 (MQ=255)

                                     ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                      <  2:2558747/100‑1 (MQ=255)

                                       tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                    <  1:1614843/100‑1 (MQ=255)

                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                  >  2:149095/1‑100 (MQ=255)

                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                 >  2:2881997/1‑100 (MQ=255)

                                            tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                               <  2:120305/100‑1 (MQ=255)

                                             tcTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                                              <  1:2365298/99‑1 (MQ=255)

                                             aCTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                                              <  2:3038460/100‑1 (MQ=255)

                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                             <  1:1832232/100‑1 (MQ=255)

                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATcaga                                           >  2:2365275/1‑97 (MQ=255)

                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                           <  2:2863599/100‑1 (MQ=255)

                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                           <  2:876962/100‑1 (MQ=255)

                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTAGATTAATCTTCTGTGATAGTCGATCGtt                                           <  2:876879/100‑1 (MQ=255)

                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        >  2:736911/1‑100 (MQ=255)

                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        >  2:1587774/1‑100 (MQ=255)

                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        <  1:563794/100‑1 (MQ=255)

                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATAACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                         >  2:2033025/1‑99 (MQ=255)

                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                       >  1:776021/1‑100 (MQ=255)

                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                       >  2:1284665/1‑100 (MQ=255)

                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                      >  1:2740697/1‑100 (MQ=255)

                                                      cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                     <  2:1547488/100‑1 (MQ=255)

                                                        aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAtt                                   <  2:1850312/100‑1 (MQ=255)

                                                        aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑GCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                       <  2:2607447/96‑1 (MQ=255)

                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                  <  2:2207320/100‑1 (MQ=255)

                                                         cTGGCTCAACTATTACGTCTACAAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCGGTGATAGTCGATCGTTAACCGAtt                                   >  2:2176201/1‑99 (MQ=255)

                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                                >  1:510229/1‑100 (MQ=255)

                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                               <  1:2551822/100‑1 (MQ=255)

                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                               <  2:525974/100‑1 (MQ=255)

                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                              >  1:1319157/1‑100 (MQ=255)

                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                              >  2:2593309/1‑100 (MQ=255)

                                                              tCAACTATTATGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                             >  2:3013804/1‑100 (MQ=255)

                                                              tCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                             <  2:2289616/100‑1 (MQ=255)

                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAGGCGATTCAGCAcc                            <  1:2093745/100‑1 (MQ=255)

                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  1:337729/100‑1 (MQ=255)

                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  2:392879/100‑1 (MQ=255)

                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  2:3067325/100‑1 (MQ=255)

                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  2:7361/100‑1 (MQ=255)

                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  2:2717072/100‑1 (MQ=255)

                                                                    aTTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAcc                       >  2:597315/1‑100 (MQ=255)

                                                                      tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  2:2335237/1‑100 (MQ=255)

                                                                      tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  1:1339636/1‑100 (MQ=255)

                                                                           cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa                >  1:1996204/1‑100 (MQ=255)

                                                                           cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa                <  1:1853287/100‑1 (MQ=255)

                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc               >  1:1071857/1‑100 (MQ=255)

                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc               >  1:1142400/1‑100 (MQ=255)

                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc               >  1:2682750/1‑100 (MQ=255)

                                                                             aCCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  1:2183050/1‑100 (MQ=255)

                                                                             aCCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  2:1922287/1‑100 (MQ=255)

                                                                             aCCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  1:1564248/1‑100 (MQ=255)

                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc           >  1:2565699/1‑100 (MQ=255)

                                                                                 gATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg          >  1:65313/1‑100 (MQ=255)

                                                                                 gATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg          <  1:2756117/100‑1 (MQ=255)

                                                                                 gATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg          <  1:1117988/100‑1 (MQ=255)

                                                                                   tCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        <  1:326146/100‑1 (MQ=255)

                                                                                     aCGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg      >  1:1288538/1‑100 (MQ=255)

                                                                                      cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc     >  2:2747330/1‑100 (MQ=255)

                                                                                       gCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc    >  1:2398827/1‑100 (MQ=255)

                                                                                        ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt   <  1:736915/100‑1 (MQ=255)

                                                                                         cATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtt  >  1:738850/1‑100 (MQ=255)

TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215083‑4215272

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