BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,210,743	A→G	Q19R (CAA→CGA)	rpoD →	RNA polymerase sigma(70)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,210,743	0	A	G	100.0%	1191.9 / NA	372	Q19R (CAA→CGA)	rpoD	RNA polymerase sigma(70)
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (173/199); total (173/199)

Only 100 of 380 total aligned reads displayed.

CACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCA  >  NC_000913/3210653‑3210832

  ccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGc                                                                                       <  2:6496329/93‑1 (MQ=255)

  ccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGc                                                                                       >  1:2368537/1‑93 (MQ=255)

   cAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCt                                                                                      >  1:5500572/1‑93 (MQ=255)

    aaCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTAAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTa                                                                                     >  2:1170992/1‑93 (MQ=255)

     aCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTAt                                                                                    >  1:5717584/1‑93 (MQ=255)

     aCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTAt                                                                                    <  2:2861572/93‑1 (MQ=255)

      ccTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATc                                                                                   <  2:7397740/93‑1 (MQ=255)

       cTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCt                                                                                  >  2:1345500/1‑93 (MQ=255)

       cTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCt                                                                                  >  2:6073875/1‑93 (MQ=255)

        tCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTg                                                                                 >  2:7155367/1‑93 (MQ=255)

        tCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTg                                                                                 >  2:3300321/1‑93 (MQ=255)

            gAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCt                                                                             <  2:2991008/93‑1 (MQ=255)

             aaaTAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTa                                                                            >  2:5087821/1‑93 (MQ=255)

                tAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGc                                                                         <  1:1079584/93‑1 (MQ=255)

                tAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGc                                                                         <  1:6329179/93‑1 (MQ=255)

                tAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGc                                                                         >  2:4124001/1‑93 (MQ=255)

                 aaGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGcc                                                                        <  1:3364855/93‑1 (MQ=255)

                 aaGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGcc                                                                        >  2:4800561/1‑93 (MQ=255)

                 aaGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGcc                                                                        >  2:5789702/1‑93 (MQ=255)

                  aGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCg                                                                       >  2:1624476/1‑93 (MQ=255)

                   gtgtGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGa                                                                      <  1:204827/93‑1 (MQ=255)

                      tgGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGt                                                                   <  2:7670871/93‑1 (MQ=255)

                         aTACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCaa                                                                <  2:3709490/93‑1 (MQ=255)

                            ccGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGa                                                             <  1:3855568/93‑1 (MQ=255)

                             cGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGAc                                                            <  2:7805071/93‑1 (MQ=255)

                              gTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGAcc                                                           >  2:3070433/1‑93 (MQ=255)

                               tCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCa                                                          >  1:3400419/1‑93 (MQ=255)

                                cTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCAt                                                         <  2:2080637/93‑1 (MQ=255)

                                 ttATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATc                                                        <  1:1819039/93‑1 (MQ=255)

                                 ttATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATc                                                        <  2:2475420/93‑1 (MQ=255)

                                  tATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCt                                                       <  1:6646709/93‑1 (MQ=255)

                                  tATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCt                                                       >  1:3081977/1‑93 (MQ=255)

                                   aTGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTg                                                      <  1:7722780/93‑1 (MQ=255)

                                   aTGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTg                                                      <  1:649244/93‑1 (MQ=255)

                                    tGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGc                                                     <  2:2561356/93‑1 (MQ=255)

                                     ggAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGcc                                                    <  1:4789329/93‑1 (MQ=255)

                                     ggAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGcc                                                    >  2:3937906/1‑93 (MQ=255)

                                       aGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCgg                                                  <  1:97420/93‑1 (MQ=255)

                                        gCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGa                                                 <  1:2146977/93‑1 (MQ=255)

                                        gCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGGCAATGACCATCTGCCGGa                                                 >  1:5747332/1‑93 (MQ=255)

                                         cAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGaa                                                <  1:7137622/93‑1 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAg                                               <  1:4335584/93‑1 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAg                                               >  1:753607/1‑93 (MQ=255)

                                          aaaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATACCGAGGTCAATGACCATCTGCCGGaga                                               >  1:5385880/1‑91 (MQ=255)

                                          aaaaCCCGCAGGCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGGAg                                               >  2:1897993/1‑93 (MQ=255)

                                           aaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGa                                              <  2:6782282/93‑1 (MQ=255)

                                           aaaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGa                                              >  1:496313/1‑93 (MQ=255)

                                            aaCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGat                                             >  1:7581869/1‑93 (MQ=255)

                                             aCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGata                                            <  1:6105910/93‑1 (MQ=255)

                                             aCCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGata                                            <  1:100115/93‑1 (MQ=255)

                                              cccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGatat                                           <  1:1727764/93‑1 (MQ=255)

                                              cccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGatat                                           <  2:1965542/93‑1 (MQ=255)

                                               ccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtc                                          <  1:3690784/93‑1 (MQ=255)

                                               ccGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtc                                          >  1:1287789/1‑93 (MQ=255)

                                               ccGCAGTCACAGCTGAAACATCTTGTCCCCCGTGGTAAGGAGCGAGGCTATCTGCCCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtc                                          <  1:4640860/93‑1 (MQ=255)

                                                 gCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtcgt                                        >  2:5530195/1‑93 (MQ=255)

                                                 gCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtcgt                                        <  2:5173722/93‑1 (MQ=255)

                                                 gCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATAtcgt                                        >  1:3637110/1‑93 (MQ=255)

                                                    gTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGa                                     <  1:6841707/93‑1 (MQ=255)

                                                     tCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGAt                                    <  1:407707/93‑1 (MQ=255)

                                                      cacaGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGAtt                                   >  1:5635366/1‑93 (MQ=255)

                                                       acaGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtc                                  <  2:4430036/93‑1 (MQ=255)

                                                       acaGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtc                                  >  2:6963822/1‑93 (MQ=255)

                                                        caGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtca                                 <  1:2544345/93‑1 (MQ=255)

                                                        caGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtca                                 <  1:6609715/93‑1 (MQ=255)

                                                         aGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcag                                <  2:4348310/93‑1 (MQ=255)

                                                          gCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGTTtcaga                               >  2:5378500/1‑93 (MQ=255)

                                                          gCTGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcaga                               <  1:525579/93‑1 (MQ=255)

                                                            tGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatc                             >  2:5670953/1‑93 (MQ=255)

                                                            tGAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatc                             <  2:885448/93‑1 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGATGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            <  1:7399019/93‑1 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            >  2:3389854/1‑93 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            <  1:3770019/93‑1 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            >  1:1110906/1‑93 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            >  2:5993703/1‑93 (MQ=255)

                                                             gAAACTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatca                            <  2:4744479/93‑1 (MQ=255)

                                                               aaCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATtcagatcaga                          <  1:2648461/93‑1 (MQ=255)

                                                                aCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTcagatcagat                         <  2:1709905/93‑1 (MQ=255)

                                                                aCTTCTTGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTcagatcagat                         <  2:7294949/93‑1 (MQ=255)

                                                                 cttcttGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCagatcagatc                        <  2:7083445/93‑1 (MQ=255)

                                                                      tGTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGa                   <  1:2754626/93‑1 (MQ=255)

                                                                       gTCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAc                  <  1:3242026/93‑1 (MQ=255)

                                                                        tCACCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAca                 <  2:4605748/93‑1 (MQ=255)

                                                                          aCCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatc               >  2:4297071/1‑93 (MQ=255)

                                                                          aCCCGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatc               <  1:7620581/93‑1 (MQ=255)

                                                                           cccGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatca              <  2:224046/93‑1 (MQ=255)

                                                                            ccGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGAcatcat             >  2:4760943/1‑93 (MQ=255)

                                                                             cGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGCCCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACatcatc            <  1:5429423/93‑1 (MQ=255)

                                                                             cGTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACatcatc            >  2:6270453/1‑93 (MQ=255)

                                                                              gTGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATcc           >  1:1640614/1‑93 (MQ=255)

                                                                               tGGTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCa          >  2:4487619/1‑93 (MQ=255)

                                                                                 gTAAGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCaaa        <  1:2587485/93‑1 (MQ=255)

                                                                                    aGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGa     <  1:5986551/93‑1 (MQ=255)

                                                                                    aGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGa     >  2:4898375/1‑93 (MQ=255)

                                                                                    aGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGa     >  1:5165407/1‑93 (MQ=255)

                                                                                    aGGAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGa     <  2:6419720/93‑1 (MQ=255)

                                                                                     ggAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGAt    >  2:2697269/1‑93 (MQ=255)

                                                                                     ggAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGAt    <  1:923937/93‑1 (MQ=255)

                                                                                      gAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATc   <  1:4243872/93‑1 (MQ=255)

                                                                                      gAGCGAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATc   <  2:6549196/93‑1 (MQ=255)

CACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGCAAGGCTATCTGACCTATGCCGAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCA  >  NC_000913/3210653‑3210832

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