BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,793,037	A→G	Q161R (CAG→CGG)	rfaF →	ADP‑heptose‑‑lps heptosyltransferase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,793,037	0	A	G	100.0%	884.0 / NA	272	Q161R (CAG→CGG)	rfaF	ADP‑heptose‑‑lps heptosyltransferase II
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (135/137); total (135/137)

Only 100 of 283 total aligned reads displayed.

TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126

tggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                         >  2:398415/1‑93 (MQ=255)

   tggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGa                                                                                      <  2:5464280/93‑1 (MQ=255)

    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     <  1:4098332/93‑1 (MQ=255)

    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     <  1:469104/93‑1 (MQ=255)

    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     >  1:1517603/1‑93 (MQ=255)

    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     <  1:3053967/93‑1 (MQ=255)

    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     <  1:2262039/93‑1 (MQ=255)

     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                    <  2:1241077/93‑1 (MQ=255)

     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                    >  2:1277829/1‑93 (MQ=255)

       aCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGa                                                                                  <  1:604788/93‑1 (MQ=255)

         gCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAg                                                                                >  1:2164035/1‑93 (MQ=255)

          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               <  2:3360520/93‑1 (MQ=255)

          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               <  1:826356/93‑1 (MQ=255)

          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               >  2:919441/1‑93 (MQ=255)

           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              <  2:1428049/93‑1 (MQ=255)

           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              <  2:3488345/93‑1 (MQ=255)

           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              >  2:3662867/1‑93 (MQ=255)

           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              >  1:1634681/1‑93 (MQ=255)

           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              >  2:4197211/1‑93 (MQ=255)

            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             <  2:3246291/93‑1 (MQ=255)

            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             <  1:1314254/93‑1 (MQ=255)

             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            >  1:122003/1‑93 (MQ=255)

             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            >  1:5777173/1‑93 (MQ=255)

             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            <  1:5260609/93‑1 (MQ=255)

              ataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaa                                                                           <  1:5582428/93‑1 (MQ=255)

               taGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaa                                                                          <  2:3870329/93‑1 (MQ=255)

               taGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaa                                                                          >  1:6174944/1‑93 (MQ=255)

                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                        >  2:5280918/1‑93 (MQ=255)

                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                        >  2:152834/1‑93 (MQ=255)

                   gcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATc                                                                      >  2:4306788/1‑93 (MQ=255)

                    cTGGCCTATGACAAAGGCATTATGCGCACAGCATAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCa                                                                     <  1:3981165/93‑1 (MQ=255)

                    cTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCa                                                                     >  1:3273902/1‑93 (MQ=255)

                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                   <  2:2062375/93‑1 (MQ=255)

                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                   <  2:1114941/93‑1 (MQ=255)

                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCaaa                                                                   >  1:3466002/1‑91 (MQ=255)

                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                  <  1:3642044/93‑1 (MQ=255)

                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                  <  2:5564605/93‑1 (MQ=255)

                       gCCTATGACAAAGACATTATGCGCACAGCACAAGATCTGCCTCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                  <  1:2811535/93‑1 (MQ=255)

                         cTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAc                                                                >  1:5102554/1‑93 (MQ=255)

                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              <  1:1077567/93‑1 (MQ=255)

                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              <  2:440014/93‑1 (MQ=255)

                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              >  2:4647675/1‑93 (MQ=255)

                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              >  1:1486513/1‑93 (MQ=255)

                            tGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTg                                                             <  1:2961885/93‑1 (MQ=255)

                             gACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCAGTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGt                                                            >  1:4440059/1‑93 (MQ=255)

                              aCAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTa                                                           >  2:891199/1‑93 (MQ=255)

                               cAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTaa                                                          <  2:768004/93‑1 (MQ=255)

                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         >  1:6485/1‑93 (MQ=255)

                                 aaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATc                                                        >  1:2450353/1‑93 (MQ=255)

                                  aGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCa                                                       >  1:4372702/1‑93 (MQ=255)

                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                      >  2:1475535/1‑93 (MQ=255)

                                    gCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCTCAGTGGCGGGTGAGCGAAGGTGAAAAGTCATATACCTGTAATCAAt                                                     <  2:3043738/93‑1 (MQ=255)

                                     cATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtt                                                    <  1:427429/93‑1 (MQ=255)

                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 >  2:3713897/1‑93 (MQ=255)

                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 >  2:2447836/1‑93 (MQ=255)

                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 >  2:170361/1‑93 (MQ=255)

                                         aTGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCg                                                <  1:5031223/93‑1 (MQ=255)

                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  1:5178062/93‑1 (MQ=255)

                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  2:2819686/93‑1 (MQ=255)

                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                            <  1:5518684/93‑1 (MQ=255)

                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                            <  1:4181236/93‑1 (MQ=255)

                                              cacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtc                                           <  2:5030474/93‑1 (MQ=255)

                                                    acaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGa                                     >  1:3374718/1‑93 (MQ=255)

                                                    acaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGa                                     <  1:2295401/93‑1 (MQ=255)

                                                      aaGATCTGCCGCAGCCATTATTATGGCCTCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   <  2:2978765/93‑1 (MQ=255)

                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                  >  2:3754310/1‑93 (MQ=255)

                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                 <  2:471061/93‑1 (MQ=255)

                                                         aTCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTc                                <  2:2217922/93‑1 (MQ=255)

                                                          tCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAGTTTTCGCTTTCATCAGAACGTcc                               >  2:1851560/1‑93 (MQ=255)

                                                          tCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTcc                               >  1:211676/1‑93 (MQ=255)

                                                          tCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTcc                               <  1:3789541/93‑1 (MQ=255)

                                                           cTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCg                              >  1:2587240/1‑93 (MQ=255)

                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGCACGTCCga                             >  1:5876509/1‑93 (MQ=255)

                                                             gCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                            >  2:1915383/1‑93 (MQ=255)

                                                             gCCGCAGCCATTGTTATGGACGCAGTTGCGGGGGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                            >  2:5913763/1‑93 (MQ=255)

                                                               cGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatga                          <  2:1043295/93‑1 (MQ=255)

                                                               cGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatga                          <  1:321896/93‑1 (MQ=255)

                                                               cGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatga                          >  1:1740063/1‑93 (MQ=255)

                                                                gCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGAGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatgat                         >  1:4077689/1‑93 (MQ=255)

                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgg                      >  1:2034206/1‑93 (MQ=255)

                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgg                      >  2:2287472/1‑93 (MQ=255)

                                                                      aTTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttt                   >  1:2702691/1‑93 (MQ=255)

                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  <  1:3529842/93‑1 (MQ=255)

                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  <  2:1603214/93‑1 (MQ=255)

                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  <  2:6084741/93‑1 (MQ=255)

                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  >  2:3733182/1‑93 (MQ=255)

                                                                           tataGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGcc              >  2:6282641/1‑93 (MQ=255)

                                                                           tATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGcc              <  2:4096917/93‑1 (MQ=255)

                                                                            aTGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGccc             <  1:1590341/93‑1 (MQ=255)

                                                                            aTGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGccc             <  2:4879504/93‑1 (MQ=255)

                                                                            aTGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGccc             <  1:5013437/93‑1 (MQ=255)

                                                                             tGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCg            >  2:1680526/1‑93 (MQ=255)

                                                                               gCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          <  1:919202/93‑1 (MQ=255)

                                                                                 cGCAGTGGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTg        >  2:5056412/1‑93 (MQ=255)

                                                                                  gCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGc       >  2:3141121/1‑93 (MQ=255)

                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    >  2:5770531/1‑93 (MQ=255)

                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    >  1:3151029/1‑93 (MQ=255)

                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    >  2:840789/1‑93 (MQ=255)

                                                                                      ttGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAg   <  2:2379339/93‑1 (MQ=255)

                                                                                       tGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGt  <  1:5587285/93‑1 (MQ=255)

TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126

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