BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,704,879	T→A	6.5%	Y144F (TAC→TTC) ‡	lepB ←	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,704,879	0	T	A	6.5%	438.7 / 13.0	137	Y144F (TAC→TTC) ‡	lepB	signal peptidase I
Reads supporting (aligned to +/- strand): ref base T (66/62); new base A (4/5); total (70/67)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CGAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGT  >  NC_000913/2704730‑2705027

cGAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGt                                                                                                                                                      <  1:640058/150‑1 (MQ=255)

cGAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGt                                                                                                                                                      <  2:1964463/150‑1 (MQ=255)

cGAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGt                                                                                                                                                      <  1:111146/150‑1 (MQ=255)

  aaaTCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTaa                                                                                                                                                    >  2:12812/1‑150 (MQ=255)

    aTCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATc                                                                                                                                                  <  1:730584/150‑1 (MQ=255)

      cGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCaa                                                                                                                                                >  1:1637164/1‑150 (MQ=255)

         tCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCt                                                                                                                                             >  1:1555088/1‑150 (MQ=255)

          cGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAAt                                                                                                                                                   >  1:173432/1‑143 (MQ=255)

           ggTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTGACACGCCTGGCCGGAACTGCATCCCGGGTGTATCGTCAGCTCTTTTGAGACCGGTTCGTAAGTGACTTTATGGGCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgcg                                                                                                                                           <  1:1291512/150‑4 (MQ=255)

           ggTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCttt                                                                                                                                           <  1:1006990/150‑1 (MQ=255)

            gTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgcgc                                                                                                                                          <  1:1075058/150‑5 (MQ=255)

            gTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgcgc                                                                                                                                          <  2:2071172/150‑5 (MQ=255)

            gTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTg                                                                                                                                          >  2:372761/1‑150 (MQ=255)

             ttCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTgg                                                                                                                                         >  2:306008/1‑150 (MQ=255)

              tCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgcgcgg                                                                                                                                        <  1:672240/150‑7 (MQ=255)

              tCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGa                                                                                                                                        <  1:563398/150‑1 (MQ=255)

               cgaCGTTTTAGTAGGTGACCGGCAGCGCGTTTTCAGACGCTTGGGTGTAAATGCATCCCGGTTGAATCGTAAGCTCTTTTGAGTCGGTATCGTAAGTGACTTTATAGGGCGGTAAACTCTCCGTGGTTTTTATGTAATCAAGCTTTGGAt                                                                                                                                       <  2:1339015/148‑1 (MQ=255)

               ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCCCCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                       >  1:1969950/1‑150 (MQ=255)

               ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                       <  2:139075/150‑1 (MQ=255)

               ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                       >  1:1997174/1‑150 (MQ=255)

               ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                       >  1:956302/1‑150 (MQ=255)

               ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                       <  2:164589/150‑1 (MQ=255)

               ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                       >  1:1394239/1‑150 (MQ=255)

                   gTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTc                                                                                                                                   >  1:8609/1‑150 (MQ=255)

                    tttGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTcc                                                                                                                                  >  1:1149035/1‑150 (MQ=255)

                      tGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                                                                                                                >  1:1510082/1‑150 (MQ=255)

                       gAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGa                                                                                                                               >  2:565224/1‑150 (MQ=255)

                            ggTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATAttt                                                                                                                          >  2:1804596/1‑150 (MQ=255)

                                aCCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAg                                                                                                                      <  2:253285/150‑1 (MQ=255)

                                  cGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAc                                                                                                                    >  1:859346/1‑150 (MQ=255)

                                  cGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAc                                                                                                                    >  2:7019/1‑150 (MQ=255)

                                  cGGCAGCGCG‑TTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                   <  1:1223131/150‑1 (MQ=255)

                                   ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGCCTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                   >  1:191498/1‑150 (MQ=255)

                                       gcgcGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGa                                                                                                               >  2:2350253/1‑150 (MQ=255)

                                          cgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGatat                                                                                                            >  2:1796621/1‑150 (MQ=255)

                                                   acGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGtt                                                                                                   <  1:7019/150‑1 (MQ=255)

                                                      ccTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCg                                                                                                >  2:2116377/1‑150 (MQ=255)

                                                       cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTTGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCgg                                                                                               <  1:2281307/150‑1 (MQ=255)

                                                       cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCGTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCgg                                                                                               <  2:1924535/150‑1 (MQ=255)

                                                        tGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGa                                                                                              <  1:1851369/150‑1 (MQ=255)

                                                         ggCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGAt                                                                                             >  2:256858/1‑150 (MQ=255)

                                                           ccGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGa                                                                                           <  1:879927/150‑1 (MQ=255)

                                                           ccGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGa                                                                                           >  2:936439/1‑150 (MQ=255)

                                                            cGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGAc                                                                                          <  2:709310/150‑1 (MQ=255)

                                                            cGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGAc                                                                                          <  1:204680/150‑1 (MQ=255)

                                                              gAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCg                                                                                        <  2:600054/150‑1 (MQ=255)

                                                                 cTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGtt                                                                                     <  1:776505/150‑1 (MQ=255)

                                                                   gCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTc                                                                                   >  2:2352339/1‑150 (MQ=255)

                                                                   gCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATATTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTc                                                                                   >  1:1624000/1‑150 (MQ=255)

                                                                      tCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCCCGATATCGCCGCGTTTCGGATGACCGGTTTCGAt                                                                                >  1:1008495/1‑150 (MQ=255)

                                                                      tCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGAt                                                                                >  1:1056679/1‑150 (MQ=255)

                                                                           gTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCg                                                                           <  2:1997174/150‑1 (MQ=255)

                                                                           gTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCg                                                                           >  2:2206236/1‑150 (MQ=255)

                                                                              gAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGttt                                                                        >  1:1904983/1‑150 (MQ=255)

                                                                                 tCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCt                                                                     <  2:1719665/150‑1 (MQ=255)

                                                                                   gTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTgg                                                                   <  1:565224/150‑1 (MQ=255)

                                                                                   gTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTgg                                                                   <  1:473712/150‑1 (MQ=255)

                                                                                    tCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGt                                                                  <  2:956302/150‑1 (MQ=255)

                                                                                       gCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGa                                                               <  1:125627/150‑1 (MQ=255)

                                                                                       gCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGa                                                               >  1:1800494/1‑150 (MQ=255)

                                                                                        ctctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGAt                                                              <  2:1555088/150‑1 (MQ=255)

                                                                                         tctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATa                                                             <  1:1017520/150‑1 (MQ=255)

                                                                                         tctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATa                                                             <  1:2053190/150‑1 (MQ=255)

                                                                                         tctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATa                                                             <  1:372761/150‑1 (MQ=255)

                                                                                          ctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAg                                                            >  1:1312492/1‑150 (MQ=255)

                                                                                          ctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAg                                                            >  1:883335/1‑150 (MQ=255)

                                                                                             ttGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGAt                                                         >  1:583109/1‑150 (MQ=255)

                                                                                               gagaCCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCt                                                       >  2:382453/1‑150 (MQ=255)

                                                                                                 gaCCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCttt                                                     <  1:1521024/150‑1 (MQ=255)

                                                                                                   ccGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTaa                                                   >  1:2040944/1‑150 (MQ=255)

                                                                                                   ccGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTaa                                                   >  2:575385/1‑150 (MQ=255)

                                                                                                   ccGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTaa                                                   >  2:733321/1‑150 (MQ=255)

                                                                                                                gACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCaaa                                      >  1:1174818/1‑150 (MQ=255)

                                                                                                                  cTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACt                                    >  1:2397505/1‑150 (MQ=255)

                                                                                                                  cTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACt                                    >  2:2460627/1‑150 (MQ=255)

                                                                                                                  cTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACt                                    <  2:2219337/150‑1 (MQ=255)

                                                                                                                    ttATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtc                                  <  2:1149035/150‑1 (MQ=255)

                                                                                                                    ttATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtc                                  >  2:963797/1‑150 (MQ=255)

                                                                                                                     tATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtct                                 >  1:1761155/1‑150 (MQ=255)

                                                                                                                          cccGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTAcc                            <  2:318991/150‑1 (MQ=255)

                                                                                                                             ggTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGa                         >  1:657999/1‑150 (MQ=255)

                                                                                                                                aaaCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTa                              <  2:204657/142‑1 (MQ=255)

                                                                                                                                aaaCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATa                      <  2:792182/150‑1 (MQ=255)

                                                                                                                                 aaCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATaa                     >  1:473318/1‑150 (MQ=255)

                                                                                                                                  aCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATaaa                    >  1:327994/1‑150 (MQ=255)

                                                                                                                                  aCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCATCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGATTTCTGGAAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATaaa                    >  2:2325776/1‑150 (MQ=255)

                                                                                                                                      atccaaagcTTGATTACATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCa                >  1:2481285/8‑150 (MQ=255)

                                                                                                                                      atccaaagcTTGATTACATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCa                >  2:1502349/8‑150 (MQ=255)

                                                                                                                                      aCCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCa                <  1:2361117/150‑1 (MQ=255)

                                                                                                                                        cgcgcgCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCAcc              <  1:1804596/150‑1 (MQ=255)

                                                                                                                                         gcgcgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCa             <  1:694182/150‑1 (MQ=255)

                                                                                                                                         caaagcTTGATTACATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCa             >  2:1345024/5‑150 (MQ=255)

                                                                                                                                            cgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAAtt          <  1:2309421/150‑1 (MQ=255)

                                                                                                                                            cgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAAtt          <  2:1018563/150‑1 (MQ=255)

                                                                                                                                            cgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAAtt          <  2:1904983/150‑1 (MQ=255)

                                                                                                                                            cgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAAtt          <  1:727957/150‑1 (MQ=255)

                                                                                                                                              cTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTaa        >  2:1950239/1‑150 (MQ=255)

                                                                                                                                                tGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACa      <  2:1120548/150‑1 (MQ=255)

                                                                                                                                                   tGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACagag   >  1:320494/1‑150 (MQ=255)

                                                                                                                                                    gTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGt  >  2:1657061/1‑150 (MQ=255)

CGAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGT  >  NC_000913/2704730‑2705027

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