BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	321,039	G→A	D70N (GAC→AAC)	ykgE →	putative dehydrogenase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	321,039	0	G	A	100.0%	790.5 / NA	211	D70N (GAC→AAC)	ykgE	putative dehydrogenase subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (119/92); total (119/92)

Only 100 of 219 total aligned reads displayed.

GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCA  >  NC_000913/320949‑321128

   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                      >  1:1130293/1‑93 (MQ=255)

   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                      >  2:2513159/1‑93 (MQ=255)

   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                      >  1:1869320/1‑93 (MQ=255)

     ctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATcc                                                                                    <  1:5105439/93‑1 (MQ=255)

      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATccc                                                                                   >  1:4388934/1‑93 (MQ=255)

      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATccc                                                                                   <  2:3182450/93‑1 (MQ=255)

      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATccc                                                                                   >  1:2265623/1‑93 (MQ=255)

         ggTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatt                                                                                <  1:1963368/93‑1 (MQ=255)

         ggTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatt                                                                                <  2:3359278/93‑1 (MQ=255)

          gTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatta                                                                               <  2:3769788/93‑1 (MQ=255)

          gTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatta                                                                               >  2:5673682/1‑93 (MQ=255)

           tCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCattat                                                                              <  2:2789503/93‑1 (MQ=255)

            cAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCattatt                                                                             <  2:5542157/93‑1 (MQ=255)

             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            >  2:4410926/1‑93 (MQ=255)

             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            <  1:4403687/93‑1 (MQ=255)

             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            <  2:2480930/93‑1 (MQ=255)

             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                            <  1:6084703/93‑1 (MQ=255)

              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                           >  1:1545445/1‑93 (MQ=255)

              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                           <  2:5790610/93‑1 (MQ=255)

              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                           <  1:622506/93‑1 (MQ=255)

              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                           >  2:669917/1‑93 (MQ=255)

               ccTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCa                                                                          >  2:4648855/1‑93 (MQ=255)

               ccTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCa                                                                          <  1:6131721/93‑1 (MQ=255)

                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        >  2:3521446/1‑93 (MQ=255)

                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        >  1:951787/1‑93 (MQ=255)

                  gCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCg                                                                       >  2:4142862/1‑93 (MQ=255)

                   cGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCgg                                                                      >  1:4909832/1‑93 (MQ=255)

                    gATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggc                                                                     >  1:3718415/1‑93 (MQ=255)

                     aTCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggcg                                                                    >  2:457721/1‑92 (MQ=255)

                      tCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctg                                                                   >  2:1122746/1‑93 (MQ=255)

                       cAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctgg                                                                  >  2:3402559/1‑93 (MQ=255)

                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                 <  1:4455999/93‑1 (MQ=255)

                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                 <  2:1038453/93‑1 (MQ=255)

                         aTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggct                                                                >  1:5445226/1‑93 (MQ=255)

                          tAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctc                                                               <  1:234783/93‑1 (MQ=255)

                          tAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctc                                                               <  1:1932213/93‑1 (MQ=255)

                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              >  2:5299768/1‑93 (MQ=255)

                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              <  2:1407385/93‑1 (MQ=255)

                              ggTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGc                                                           >  2:4306280/1‑93 (MQ=255)

                              ggTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGc                                                           >  1:4818782/1‑93 (MQ=255)

                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCa                                                          <  2:5087834/93‑1 (MQ=255)

                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCa                                                          >  1:3007669/1‑93 (MQ=255)

                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCa                                                          >  1:1199025/1‑93 (MQ=255)

                                ttATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCAc                                                         >  1:852917/1‑93 (MQ=255)

                                    atCAAAGAGGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     >  2:6199767/1‑93 (MQ=255)

                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     <  1:5844447/93‑1 (MQ=255)

                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     <  1:3213464/93‑1 (MQ=255)

                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     >  2:4149135/1‑93 (MQ=255)

                                      cAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGc                                                   >  1:285563/1‑93 (MQ=255)

                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                  >  2:1705047/1‑93 (MQ=255)

                                         agaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGt                                                >  1:443276/1‑93 (MQ=255)

                                         agaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGt                                                <  1:2013586/93‑1 (MQ=255)

                                         agaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGt                                                >  2:6020195/1‑93 (MQ=255)

                                          gaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTa                                               >  2:4359815/1‑93 (MQ=255)

                                          gaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTa                                               <  2:2363964/93‑1 (MQ=255)

                                            agCGATTCCAGGGATGAAAAATCTGATAGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaa                                             >  2:714842/1‑93 (MQ=255)

                                             gCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaa                                            >  2:1354821/1‑93 (MQ=255)

                                              cGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaaa                                           <  1:3534586/93‑1 (MQ=255)

                                                aTTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGt                                         >  2:502938/1‑93 (MQ=255)

                                                aTTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGt                                         >  2:1465397/1‑93 (MQ=255)

                                                 ttCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGtt                                        >  2:767068/1‑93 (MQ=255)

                                                 ttCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGtt                                        <  2:5493741/93‑1 (MQ=255)

                                                   ccAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAc                                      >  1:5703633/1‑93 (MQ=255)

                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  2:1644177/93‑1 (MQ=255)

                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     >  1:2655470/1‑93 (MQ=255)

                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     >  1:1393372/1‑93 (MQ=255)

                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     >  1:4507734/1‑93 (MQ=255)

                                                    cAGCGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  1:1894270/93‑1 (MQ=255)

                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                    <  2:828836/93‑1 (MQ=255)

                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                    <  2:3319527/93‑1 (MQ=255)

                                                     aGGGATGAAAAATCTGATAGCCGCACTGGAGGATGACAACGATCCCATTATTTCACCGGCTGGCTCATGCACCTATGCCGTAAAAAGATAccc                                    >  2:403619/1‑93 (MQ=255)

                                                       ggATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGa                                  <  1:1069718/93‑1 (MQ=255)

                                                        gATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGAc                                 <  2:5747008/93‑1 (MQ=255)

                                                        gATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGAc                                 >  1:4611406/1‑93 (MQ=255)

                                                         aTGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACg                                >  1:5279053/1‑93 (MQ=255)

                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                               >  1:4324711/1‑93 (MQ=255)

                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                               >  1:4358024/1‑93 (MQ=255)

                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                               >  2:2776935/1‑93 (MQ=255)

                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTa                              >  1:4523485/1‑93 (MQ=255)

                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                            <  2:191933/93‑1 (MQ=255)

                                                              aaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCt                           <  1:3496528/93‑1 (MQ=255)

                                                               aaTCTGATCGCCGCACTGGAGGATAACAACGATCCCACTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTg                          >  1:5620751/1‑93 (MQ=255)

                                                                 tCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGc                        <  2:88069/93‑1 (MQ=255)

                                                                   tGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCgg                      >  2:2585537/1‑93 (MQ=255)

                                                                   tGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCgg                      >  2:2063971/1‑93 (MQ=255)

                                                                      tCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATg                   <  1:5315644/93‑1 (MQ=255)

                                                                      tCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATg                   <  1:123071/93‑1 (MQ=255)

                                                                      tCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATg                   <  1:5998027/93‑1 (MQ=255)

                                                                        gccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGaa                 >  1:4007071/1‑93 (MQ=255)

                                                                        gccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGaa                 <  1:917032/93‑1 (MQ=255)

                                                                          cgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAcc               >  1:5736202/1‑93 (MQ=255)

                                                                            cACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTg             <  2:5339299/93‑1 (MQ=255)

                                                                            cACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTg             <  1:3811118/93‑1 (MQ=255)

                                                                             aCTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGa            >  2:3661109/1‑93 (MQ=255)

                                                                               tGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAAt          >  1:4256239/1‑93 (MQ=255)

                                                                                ggaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATg         >  1:5510150/1‑93 (MQ=255)

                                                                                  aggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATggg       <  2:2322752/93‑1 (MQ=255)

                                                                                  aggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATggg       <  2:2942490/93‑1 (MQ=255)

                                                                                   ggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGAGGATGAACCTGAATGGGc      <  1:575690/93‑1 (MQ=255)

                                                                                      tAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATc   >  2:3891167/1‑93 (MQ=255)

GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCA  >  NC_000913/320949‑321128

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