BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,698	A→G	5.1%	H93R (CAC→CGC) ‡	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,180,698	0	A	G	5.1%	415.7 / 8.8	138	H93R (CAC→CGC) ‡	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base A (62/69); new base G (2/5); total (64/74)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 139 total aligned reads displayed.

CCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCGGT  >  NC_000913/3180548‑3180846

cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTG‑CTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCa                                                                                                                                                      >  1:835542/1‑150 (MQ=255)

   aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAca                                                                                                                                                    <  2:1272637/150‑2 (MQ=255)

   aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGAAcc                                                                                                                                                    >  1:1040050/1‑150 (MQ=255)

    ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAcaa                                                                                                                                                   <  1:441055/150‑3 (MQ=255)

     gTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAcaaa                                                                                                                                                  <  1:1627184/150‑4 (MQ=255)

      ttGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGc                                                                                                                                                 <  1:813721/150‑1 (MQ=255)

        gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGCTGACCAAATTCCGCATTCCAAAACACACATGGAGTTTTGTGCGCCAGTCCTGGCTGACGCACCAGCCa                                                                                                                                               >  1:232476/1‑150 (MQ=255)

         aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                                                                              >  2:1708303/1‑150 (MQ=255)

         aaaaGCTGGAAGAAGGCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAt                                                                                                                                              <  1:723966/150‑1 (MQ=255)

          aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtg                                                                                                                                             <  2:1583995/150‑2 (MQ=255)

            aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtgac                                                                                                                                           <  2:12029/150‑4 (MQ=255)

            aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGc                                                                                                                                           <  1:1695381/150‑1 (MQ=255)

              cTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCtt                                                                                                                                         <  1:342144/150‑1 (MQ=255)

               tGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                        <  2:1445972/150‑1 (MQ=255)

                ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTTATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                       >  1:850586/1‑150 (MQ=255)

                ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGCAGTTTTGTGCGCCAGTCATGGATTACGCACCAGCCATCGCTTTa                                                                                                                                       >  1:1125780/1‑150 (MQ=255)

                  aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                                                                     <  2:565879/150‑1 (MQ=255)

                    gaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGCGGAAAAAGTGATCGCAAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                                   <  2:1600789/150‑1 (MQ=255)

                    gaagGCACCGCCGAAATGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTTATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcg                                                                                                                                   <  2:833072/150‑1 (MQ=255)

                     aagGCACCGCCGAACTGCACCTGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAACAACTCCTGGAGTTATGTGCGCCAGTCATGGCTGTCTCACCAGACATCGCTTTATTcgc                                                                                                                                  >  2:887440/1‑150 (MQ=255)

                        tCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTc                                                                                                                               <  1:742993/150‑1 (MQ=255)

                         cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                              <  2:402851/150‑1 (MQ=255)

                         cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCt                                                                                                                              >  2:1151097/1‑150 (MQ=255)

                           ccgccgAACTGCACCATATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAACCATCGCTTTATTCGCGTCTTg                                                                                                                            <  2:596239/150‑1 (MQ=255)

                           ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                            <  2:31202/150‑1 (MQ=255)

                           ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                            <  2:411449/150‑1 (MQ=255)

                           ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                            <  2:325101/150‑1 (MQ=255)

                            cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGTTCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGa                                                                                                                           >  2:501654/1‑150 (MQ=255)

                            cgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGa                                                                                                                           >  2:1651773/1‑150 (MQ=255)

                             gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGACTGACGCACCAGCCATCGCtt                                                                                                                                         <  1:1698737/135‑1 (MQ=255)

                             gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGACTGACGCACCAGCCATCGCtt                                                                                                                                         >  2:1698737/1‑135 (MQ=255)

                              ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                         >  1:1010320/1‑150 (MQ=255)

                              ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                         >  2:1656406/1‑150 (MQ=255)

                              ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                                                                         <  1:679400/150‑1 (MQ=255)

                                gAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTg                                                                                                                       <  1:901200/150‑1 (MQ=255)

                                 taCTGCACCAGATGTGCAGGAAAGTGTTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTgg                                                                                                                      <  1:205325/149‑1 (MQ=255)

                                   cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                                    <  1:179522/150‑1 (MQ=255)

                                   cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                                                                    <  2:976714/150‑1 (MQ=255)

                                    tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGt                                                                                                                   >  2:421062/1‑150 (MQ=255)

                                    tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGt                                                                                                                   >  1:837131/1‑150 (MQ=255)

                                     gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                                  <  2:305259/150‑1 (MQ=255)

                                     gCACCAGATGAGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCAGATGCCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGCCGCACCAGCCATCGCTTTATTCGCGACTTGATCTGGCGTg                                                                                                                  >  2:100551/1‑150 (MQ=255)

                                       aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                                                                >  1:1311907/1‑150 (MQ=255)

                                       aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                                                                <  1:933931/150‑1 (MQ=255)

                                          aGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGCCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATg                                                                                                             >  2:892061/1‑150 (MQ=255)

                                          aGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATg                                                                                                             <  1:547901/150‑1 (MQ=255)

                                           gATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                                                                            >  1:47726/1‑150 (MQ=255)

                                                gCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTg                                                                                                       <  1:272265/150‑1 (MQ=255)

                                                 ccTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTgg                                                                                                      <  1:1708303/150‑1 (MQ=255)

                                                  cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCTAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                     >  2:897896/1‑150 (MQ=255)

                                                  cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                     >  2:987270/1‑150 (MQ=255)

                                                    gAAAGTGGTGGAAAAAGTGAACGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGa                                                                                                   <  2:1190738/150‑1 (MQ=255)

                                                     aaaGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGaa                                                                                                  >  1:343697/1‑150 (MQ=255)

                                                       aGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAAcc                                                                                                <  1:1656406/150‑1 (MQ=255)

                                                            tggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGAGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                                                                           >  1:194634/1‑150 (MQ=255)

                                                            gggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                                                                           <  2:1578808/149‑1 (MQ=255)

                                                             ggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACt                                                                                          <  2:654383/150‑1 (MQ=255)

                                                               aaaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTc                                                                                        <  2:906081/150‑1 (MQ=255)

                                                                aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTTATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                       <  2:850586/150‑1 (MQ=255)

                                                                aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                       <  2:806871/150‑1 (MQ=255)

                                                                     tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaa                                                                                  >  1:923605/1‑150 (MQ=255)

                                                                      gATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaaa                                                                                 <  2:773569/150‑1 (MQ=255)

                                                                           ccAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGGTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACg                                                                            <  1:1155798/150‑1 (MQ=255)

                                                                             aGCGATGAGCTTATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTTTGCGCCATTCATGGCTGACTCACCAGCCATCGCTTTATTCTCGTCTTGATCTGGCGTGGGATGTCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                          <  1:1613359/150‑1 (MQ=255)

                                                                                 aTGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGGGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATa                                                                      >  1:143713/1‑150 (MQ=255)

                                                                                     gCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGcc                                                                  <  2:187030/150‑1 (MQ=255)

                                                                                           gaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTc                                                            >  1:1599092/1‑150 (MQ=255)

                                                                                            aCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCa                                                           >  1:1192690/1‑150 (MQ=255)

                                                                                             ccAAATTCCGCATTCCAACACACAGCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGGGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                          <  1:604626/150‑1 (MQ=255)

                                                                                             ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                          >  2:1469220/1‑150 (MQ=255)

                                                                                             ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                          <  1:501654/150‑1 (MQ=255)

                                                                                               aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta                                                        >  1:629368/1‑150 (MQ=255)

                                                                                                aaTTCCGGATTCCAAAACACACCTGGAGTTTTGTGCGACAGTCAGGTCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGTCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtat                                                       <  1:981263/150‑1 (MQ=255)

                                                                                                aaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtat                                                       >  2:1442594/1‑150 (MQ=255)

                                                                                                 aTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtata                                                      <  1:252170/150‑1 (MQ=255)

                                                                                                 aTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACtata                                                      <  1:743090/150‑1 (MQ=255)

                                                                                                    ccGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGa                                                   >  1:863612/1‑150 (MQ=255)

                                                                                                       cATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggc                                                >  2:1353767/1‑150 (MQ=255)

                                                                                                       cATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggc                                                >  1:510650/1‑150 (MQ=255)

                                                                                                        aTTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggcg                                               >  1:1267664/1‑150 (MQ=255)

                                                                                                            cAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGt                                           <  1:1603528/150‑1 (MQ=255)

                                                                                                             aaaaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGtt                                          >  2:1553664/1‑150 (MQ=255)

                                                                                                                  acacCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAACGCTTTATTATCTTCTTTATCTGGCGTGGGATGTCAATGTTGAACATAAACTTCTGGAAAATAACGCCGATACGCCAACGTAAATATACGAGGCGGAGTTCTTTc                                     >  2:1111775/1‑150 (MQ=255)

                                                                                                                  acacCTGGAGTTGTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTc                                     <  2:908810/150‑1 (MQ=255)

                                                                                                                     ccTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGGCGTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCGTTCAGt                                  >  2:1040383/1‑150 (MQ=255)

                                                                                                                           gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                            >  1:1647728/1‑150 (MQ=255)

                                                                                                                              tggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                         >  2:413913/3‑150 (MQ=255)

                                                                                                                               tgtgCGCCAGTCAGGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggct                        <  2:193374/150‑1 (MQ=255)

                                                                                                                               ggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggct                        >  1:1279397/2‑150 (MQ=255)

                                                                                                                                gtgCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctg                       >  2:742032/1‑150 (MQ=255)

                                                                                                                                  gcgcCAGGCATGGCTGACGTACCAGCCATCGCTTTATTCGCGTCTTGATATGGCGTGGGATGGCACTGTTGATCCTAAACTTCTGGAATATATCGCGGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctggt                     <  2:232476/150‑2 (MQ=255)

                                                                                                                                   cgcCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGaa                    >  1:665347/1‑150 (MQ=255)

                                                                                                                                    taaagcgaTGGCTGGTGCGTCAGCCATCGCTTTTTTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAg                   >  1:113565/8‑150 (MQ=255)

                                                                                                                                     tcAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGa                  >  1:787973/2‑150 (MQ=255)

                                                                                                                                          cATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGc             <  1:1614785/150‑1 (MQ=255)

                                                                                                                                            tGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCtt           >  1:970417/1‑150 (MQ=255)

                                                                                                                                             ggCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTa          >  1:399884/1‑150 (MQ=255)

                                                                                                                                              gCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTaa         >  1:146471/1‑150 (MQ=255)

                                                                                                                                                   cGCACCAGCAATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCg    <  2:1279397/150‑1 (MQ=255)

                                                                                                                                                    gCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCgg   <  1:760023/150‑1 (MQ=255)

CCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACGCCGGT  >  NC_000913/3180548‑3180846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`