BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,215,147	A→C	100%	N290H (AAC→CAC)	metA →	homoserine O‑transsuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,147	0	A	C	100.0%	379.7 / ‑6.7	110	N290H (AAC→CAC)	metA	homoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand): ref base A (1/0); major base C (55/53); minor base T (0/1); total (56/54)

Only 100 of 110 total aligned reads displayed.

GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGC  >  NC_000913/4215050‑4215243

ggATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCAc                                                                                                <  1:1184314/100‑1 (MQ=255)

ggATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCAc                                                                                                <  2:1560882/100‑1 (MQ=255)

 gATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACt                                                                                               <  1:905594/100‑1 (MQ=255)

  aTGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTa                                                                                              <  1:787316/100‑1 (MQ=255)

  aTGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTCGCGTAGTCACCGTAATTTACTGTTTACCTACTGGCTCCACTa                                                                                              <  1:1611530/100‑1 (MQ=255)

   tGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTAt                                                                                             <  1:1697580/100‑1 (MQ=255)

    gTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTAtt                                                                                            <  2:1011841/100‑1 (MQ=255)

     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                           <  2:89349/100‑1 (MQ=255)

     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                           >  2:1474219/1‑100 (MQ=255)

      aCCGTATAACTCTTCCCCGTACACTGTTCCGCAACATACCCCGCTTCCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                           <  2:368519/99‑1 (MQ=255)

      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                          >  2:254453/1‑100 (MQ=255)

      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                          >  1:1529769/1‑100 (MQ=255)

       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                         <  1:1812506/100‑1 (MQ=255)

       ccGTATAACTATTTCCCGCACAAGGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                         <  1:373377/100‑1 (MQ=255)

            taACTATTTCCCGCACAATGATCCGCAAAATACACAGCGAGCGAGCTGGCGTAGTCACGGTCATTTACTGTTTCCCAACTGGCTCCACTACTACATCAAc                                                                                    >  2:2029453/1‑100 (MQ=255)

             aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAcc                                                                                   <  1:1556105/100‑1 (MQ=255)

             aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAcc                                                                                   <  1:1225756/100‑1 (MQ=255)

                 aTTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGAt                                                                               >  1:428673/1‑100 (MQ=255)

                 aTTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGAt                                                                               >  1:1206142/1‑100 (MQ=255)

                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGCTCACg                                                                           >  1:1176573/1‑100 (MQ=255)

                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  2:2325868/1‑100 (MQ=255)

                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  1:654436/1‑100 (MQ=255)

                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          <  1:2165528/100‑1 (MQ=255)

                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          >  1:1278956/1‑100 (MQ=255)

                       cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACgac                                                                         <  1:1622924/100‑3 (MQ=255)

                       cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGcc                                                                         >  1:919785/1‑100 (MQ=255)

                       cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGcc                                                                         <  2:1697020/100‑1 (MQ=255)

                        gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          >  2:1443330/1‑98 (MQ=255)

                        gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCCCCACTATTACGTCTACCAGATCAc                                                                            >  1:603356/1‑96 (MQ=255)

                         cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCAt                                                                       <  1:447867/100‑1 (MQ=255)

                         cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCAt                                                                       <  2:2368747/100‑1 (MQ=255)

                          acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCa                                                                        <  2:2082/98‑1 (MQ=255)

                          acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATa                                                                      <  1:718065/100‑1 (MQ=255)

                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATAc                                                                     >  2:2211173/1‑100 (MQ=255)

                            aaTGATCCGCAAAATACTCCGCGAGCGAGCTGACGTAGTGACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACg                                                                    <  2:1676668/100‑1 (MQ=255)

                             aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCTATACGa                                                                   >  2:481697/1‑100 (MQ=255)

                              tGATCCGCAAAATACACCGCGATCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGAt                                                                  <  2:1639081/100‑1 (MQ=255)

                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 <  2:1801238/100‑1 (MQ=255)

                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 <  2:2104236/100‑1 (MQ=255)

                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCAAGCCATAcaa                                                                   >  2:283533/1‑96 (MQ=255)

                                aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCt                                                                >  1:386476/1‑100 (MQ=255)

                                 tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                              >  1:2168725/1‑100 (MQ=255)

                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                             >  2:721508/1‑100 (MQ=255)

                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                             <  2:16433/100‑1 (MQ=255)

                                    gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACgg                                                            <  1:659777/100‑1 (MQ=255)

                                      aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                                                          <  1:2017808/100‑1 (MQ=255)

                                      aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                                                          >  1:1965221/1‑100 (MQ=255)

                                       aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGCTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGcac                                                         <  2:1529759/100‑1 (MQ=255)

                                         aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGcacan                                                       <  1:596063/100‑2 (MQ=255)

                                            cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCCCATCaa                                                    >  2:869481/1‑100 (MQ=255)

                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                   >  2:892545/1‑99 (MQ=255)

                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  >  1:986933/1‑100 (MQ=255)

                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  >  1:205954/1‑100 (MQ=255)

                                                 cgagcgagCTGCCGTAATCACCTTAATTTACTGCTTCTCAAATGGCTTCACTTTTACTTCTCCCAGAGGACGCCATACGATCTACGGCACATGAATcc                                                 <  2:2367871/98‑1 (MQ=255)

                                                     cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCt                                           <  1:216680/100‑1 (MQ=255)

                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                       <  2:1682459/100‑1 (MQ=255)

                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                       <  1:1345027/100‑1 (MQ=255)

                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                       <  1:1640487/100‑1 (MQ=255)

                                                          tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                      >  1:2164356/1‑100 (MQ=255)

                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                      >  1:2141602/1‑99 (MQ=255)

                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                     >  2:2281564/1‑100 (MQ=255)

                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                     >  1:1394971/1‑100 (MQ=255)

                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                     <  1:185648/100‑1 (MQ=255)

                                                             cGTAGACACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCAGGATTAAt                                   >  1:1755058/1‑100 (MQ=255)

                                                              gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc                                  <  2:1112576/100‑1 (MQ=255)

                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                               <  1:1456770/100‑1 (MQ=255)

                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                               >  1:1250616/1‑100 (MQ=255)

                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGAACTACGGCACATGAATACAACGCTGGATTAAtcttcag                             >  2:2322945/1‑98 (MQ=255)

                                                                    aCGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                            >  1:505410/1‑100 (MQ=255)

                                                                     cGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                           >  1:1331066/1‑100 (MQ=255)

                                                                      ggTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                          >  2:1465838/1‑100 (MQ=255)

                                                                        tAATTTACTGTTTATCAACTGGTTCCATTATTACGTCTACCAGATCATGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                         <  2:1181211/99‑1 (MQ=255)

                                                                        tAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                         >  1:15881/1‑99 (MQ=255)

                                                                        tAATTTACTGTTTAACAACTGGCTCTACCATCACGTGGACAAGAGCACGCTATACG‑TCTACGGCACATGAAACCAACGCTGGATTAACCGTCTGTGAt                         <  2:1833969/98‑1 (MQ=255)

                                                                        tAATTTACTGGTTACCAACTCGCTCCACTATTACATCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                            <  2:963695/96‑1 (MQ=255)

                                                                         aaTTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                       <  2:2346488/100‑1 (MQ=255)

                                                                         aaTTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                       <  1:1543827/100‑1 (MQ=255)

                                                                           tttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                     <  1:1226553/100‑1 (MQ=255)

                                                                            ttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                    <  1:1727569/100‑1 (MQ=255)

                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                 >  1:333777/1‑100 (MQ=255)

                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTATGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                <  2:1325372/100‑1 (MQ=255)

                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                >  2:174750/1‑100 (MQ=255)

                                                                                 gTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATctt               >  1:1651383/1‑98 (MQ=255)

                                                                                 gTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt               >  2:1058721/1‑100 (MQ=255)

                                                                                 gTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt               <  1:1942218/100‑1 (MQ=255)

                                                                                 gTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt               >  2:537625/1‑100 (MQ=255)

                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa             <  1:977577/100‑1 (MQ=255)

                                                                                    tACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt              >  1:1867904/1‑98 (MQ=255)

                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTATGc          >  1:450848/1‑100 (MQ=255)

                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          >  2:826122/1‑100 (MQ=255)

                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          <  1:266960/100‑1 (MQ=255)

                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          >  2:508188/1‑100 (MQ=255)

                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACACTGGATTAATCTTCTGTGATAGTCGTTCGTTa             >  1:1163461/1‑97 (MQ=255)

                                                                                       cAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          >  2:1587277/1‑99 (MQ=255)

                                                                                         acgggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt       <  2:1435220/97‑1 (MQ=255)

                                                                                         aCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt       >  2:336858/1‑100 (MQ=255)

                                                                                          cTGGCTCCACTATTACGTCTACCAGATCACGCCATACGCTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt       <  2:2301972/99‑1 (MQ=255)

                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa    <  1:1861091/100‑1 (MQ=255)

                                                                                             gCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg   >  1:853679/1‑100 (MQ=255)

                                                                                              cTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCatc  >  1:823993/1‑98 (MQ=255)

GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGC  >  NC_000913/4215050‑4215243

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