BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,161,324	G→A	V143M (GTG→ATG)	baeS →	BaeS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,161,324	0	G	A	99.6%	870.6 / NA	228	V143M (GTG→ATG)	baeS	BaeS
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (95/132); minor base C (1/0); total (96/132)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

AACAACAAAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGCC  >  NC_000913/2161234‑2161413

aacaacAAAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATg                                                                                         >  1:3171369/1‑93 (MQ=255)

aacaacAAAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATg                                                                                         <  1:1396341/93‑1 (MQ=255)

  caacaaAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCgatgat                                                                                       <  1:5092922/93‑1 (MQ=255)

  caacaaAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCgatgat                                                                                       >  1:896914/1‑93 (MQ=255)

  caacaaAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCgatgat                                                                                       <  2:1924795/93‑1 (MQ=255)

   aacaaAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATc                                                                                      <  2:622819/93‑1 (MQ=255)

    acaaAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCg                                                                                     >  2:3756888/1‑93 (MQ=255)

      aaaGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGcc                                                                                   <  1:2524261/93‑1 (MQ=255)

        aGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTc                                                                                 >  2:6323183/1‑93 (MQ=255)

         gTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTcc                                                                                >  2:2518175/1‑93 (MQ=255)

         gTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTcc                                                                                <  1:624385/93‑1 (MQ=255)

           gCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTcccc                                                                              >  2:2737125/1‑93 (MQ=255)

            cTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTccccc                                                                             <  2:1553155/93‑1 (MQ=255)

             tggttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCg                                                                            >  2:1101087/1‑93 (MQ=255)

              ggttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGt                                                                           >  1:5487181/1‑93 (MQ=255)

              ggttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGt                                                                           <  2:5324171/93‑1 (MQ=255)

              ggttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGt                                                                           >  1:4589543/1‑93 (MQ=255)

              ggttggtCCACGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGt                                                                           <  2:5241440/93‑1 (MQ=255)

               gttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGtt                                                                          >  1:5807565/1‑93 (MQ=255)

                ttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTg                                                                         >  1:1013277/1‑93 (MQ=255)

                ttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTg                                                                         <  1:5990074/93‑1 (MQ=255)

                ttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTg                                                                         <  2:5124294/93‑1 (MQ=255)

                ttggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTg                                                                         <  1:1900818/93‑1 (MQ=255)

                  ggtCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAg                                                                       >  2:2078884/1‑93 (MQ=255)

                   gtCCGCGAGCGCCGATTCCACCTGACTGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTAAGc                                                                      >  2:4407342/1‑93 (MQ=255)

                    tCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCg                                                                     <  2:3080097/93‑1 (MQ=255)

                    tCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCg                                                                     <  1:6117389/93‑1 (MQ=255)

                      cgcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACTGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGt                                                                   <  2:4518257/93‑1 (MQ=255)

                      cgcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGt                                                                   <  1:2617049/93‑1 (MQ=255)

                       gcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGtt                                                                  >  2:3983098/1‑93 (MQ=255)

                       gcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGtt                                                                  <  2:2992091/93‑1 (MQ=255)

                       gcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGtt                                                                  <  1:2442514/93‑1 (MQ=255)

                       gcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGtt                                                                  >  1:3113113/1‑93 (MQ=255)

                       gcgAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGtt                                                                  >  1:4371455/1‑93 (MQ=255)

                         gAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTaa                                                                <  2:5371693/93‑1 (MQ=255)

                         gAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTaa                                                                >  2:2926234/1‑93 (MQ=255)

                         gAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTaa                                                                >  1:2695387/1‑93 (MQ=255)

                          aGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAAc                                                               >  1:945827/1‑93 (MQ=255)

                          aGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAAc                                                               >  2:2580104/1‑93 (MQ=255)

                           gcgcCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAAcg                                                              >  2:1314986/1‑93 (MQ=255)

                           gcgcCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAAcg                                                              <  1:1691076/93‑1 (MQ=255)

                              ccGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACgcgc                                                           <  1:5679548/93‑1 (MQ=255)

                              ccGATTCCACCTGACGGTACACGGCGACCCATTCTAGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGCTGAGCGGTTAACgcgc                                                           <  1:4003437/93‑1 (MQ=255)

                                gATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCaa                                                         >  2:506072/1‑93 (MQ=255)

                                  ttCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATa                                                       >  1:6379915/1‑93 (MQ=255)

                                    ccACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACt                                                     <  2:4352591/93‑1 (MQ=255)

                                    ccACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACt                                                     <  2:5912237/93‑1 (MQ=255)

                                    ccACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACt                                                     <  1:4253898/93‑1 (MQ=255)

                                      aCCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGa                                                   >  2:1533202/1‑93 (MQ=255)

                                        cTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGata                                                 <  2:1698342/93‑1 (MQ=255)

                                        cTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGata                                                 <  1:2047727/93‑1 (MQ=255)

                                          gACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATc                                               >  2:4876305/1‑93 (MQ=255)

                                            cGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCaa                                             <  1:1922747/93‑1 (MQ=255)

                                            cGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCaa                                             <  1:3892707/93‑1 (MQ=255)

                                             ggTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAAt                                            <  1:4808387/93‑1 (MQ=255)

                                             ggTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAAt                                            <  1:6015418/93‑1 (MQ=255)

                                              gTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAAtt                                           >  1:236567/1‑93 (MQ=255)

                                                acacGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTc                                         <  2:4928349/93‑1 (MQ=255)

                                                  acGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGa                                       >  1:3345789/1‑93 (MQ=255)

                                                   cGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGAt                                      <  1:553120/93‑1 (MQ=255)

                                                    ggCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATa                                     >  2:5894992/1‑93 (MQ=255)

                                                    ggCGACCCATTCTGGTCAACGGGGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATa                                     >  1:3909454/1‑93 (MQ=255)

                                                     gCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATaa                                    >  2:4598156/1‑93 (MQ=255)

                                                     gCGACCAATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATaa                                    <  1:2923564/93‑1 (MQ=255)

                                                       gACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAaac                                  <  1:371229/93‑1 (MQ=255)

                                                         cccATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAaacaa                                <  1:789634/93‑1 (MQ=255)

                                                          ccATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAaacaac                               <  1:4567343/93‑1 (MQ=255)

                                                          ccATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAaacaac                               <  2:4261719/93‑1 (MQ=255)

                                                          ccATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAaacaac                               <  1:5313463/93‑1 (MQ=255)

                                                           cATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAacaaca                              <  2:751055/93‑1 (MQ=255)

                                                           cATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAacaaca                              <  2:47475/93‑1 (MQ=255)

                                                           cATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAacaaca                              >  1:4944570/1‑93 (MQ=255)

                                                           cATTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAacaaca                              >  1:3261108/1‑93 (MQ=255)

                                                            aTTCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAg                             >  2:5681879/1‑93 (MQ=255)

                                                             ttCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGc                            <  2:6356443/93‑1 (MQ=255)

                                                             ttCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGc                            <  2:5570341/93‑1 (MQ=255)

                                                              tCTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCg                           <  2:5170759/93‑1 (MQ=255)

                                                               cTGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCgg                          <  2:2317112/93‑1 (MQ=255)

                                                                tGGTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTAGATAAACAACAGCGGc                         >  1:4973100/1‑93 (MQ=255)

                                                                 ggTCAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCa                        >  1:279270/1‑93 (MQ=255)

                                                                    cAACGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAAc                     <  1:4180446/93‑1 (MQ=255)

                                                                     aaCGGTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAAcc                    >  2:1610774/1‑93 (MQ=255)

                                                                        ggTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGc                 >  1:1965427/1‑93 (MQ=255)

                                                                        ggTGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGc                 <  2:1062009/93‑1 (MQ=255)

                                                                          tGCGGAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTg               <  1:147717/93‑1 (MQ=255)

                                                                             ggAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGtt            >  2:891460/1‑93 (MQ=255)

                                                                             ggAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGtt            >  1:617877/1‑93 (MQ=255)

                                                                              gAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTg           <  2:3336399/93‑1 (MQ=255)

                                                                              gAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGGTg           >  1:2286679/1‑93 (MQ=255)

                                                                              gAAGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCACAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTg           <  1:1910364/93‑1 (MQ=255)

                                                                               aaGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGttta          <  1:1102760/93‑3 (MQ=255)

                                                                               aaGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGa          <  2:1981514/93‑1 (MQ=255)

                                                                                aGTTGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGAt         >  1:1147689/1‑93 (MQ=255)

                                                                                   tGGCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGt      <  2:1185566/93‑1 (MQ=255)

                                                                                    ggCGCGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTc     >  1:3506041/1‑93 (MQ=255)

                                                                                     gcgcGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCTATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCg    <  2:3433772/93‑1 (MQ=255)

                                                                                     gcgcGATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCg    >  2:3924262/1‑93 (MQ=255)

                                                                                      cgcgATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGc   <  2:6208242/93‑1 (MQ=255)

                                                                                      cgcgATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGc   <  2:644407/93‑1 (MQ=255)

                                                                                       gcgATGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGcc  <  1:5893096/93‑1 (MQ=255)

AACAACAAAGTGCTGGTTGGTCCGCGAGCGCCGATTCCACCTGACGGTACACGGCGACCCATTCTGGTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGCC  >  NC_000913/2161234‑2161413

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