BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,704,880	A→C	6.7%	Y144D (TAC→GAC) ‡	lepB ←	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,704,880	0	A	C	6.7%	401.7 / 18.0	135	Y144D (TAC→GAC) ‡	lepB	signal peptidase I
Reads supporting (aligned to +/- strand): ref base A (65/61); new base C (4/5); total (69/66)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 139 total aligned reads displayed.

GAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCG  >  NC_000913/2704731‑2705029

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

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

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

       cTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGc                                                                                                                                                >  2:1447037/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

             tCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGa                                                                                                                                          >  1:709310/1‑150 (MQ=255)

             tCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGAATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGa                                                                                                                                          <  2:1135022/150‑1 (MQ=255)

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

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

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

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

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

              ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAt                                                                                                                                         <  2:139075/150‑1 (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)

                            gTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTa                                                                                                                           >  1:1870983/1‑150 (MQ=255)

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

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

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

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

                                  ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                     <  1:738676/150‑1 (MQ=255)

                                  ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                     >  1:1261557/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)

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

                                                      cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCgg                                                                                                 <  1:51090/150‑1 (MQ=255)

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

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

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

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

                                                           cGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGAc                                                                                            <  2:709310/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)

                                                                    aTCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGa                                                                                   <  2:1870983/150‑1 (MQ=255)

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

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

                                                                       ccGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCa                                                                                >  1:668426/1‑150 (MQ=255)

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

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

                                                                               aTCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTc                                                                        >  1:1549644/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                         ctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAg                                                              >  1:1312492/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:1960368/1‑150 (MQ=255)

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

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

                                                                                                   cGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAAt                                                    <  2:1610549/150‑1 (MQ=255)

                                                                                                            aGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGc                                           >  1:927125/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:963797/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                               aaaCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTa                                >  1:204657/1‑142 (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                  >  2:1502349/8‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                                             cTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACAAATTaa          >  2:1335399/1‑150 (MQ=255)

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

                                                                                                                                               tGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACa        >  1:1779311/1‑150 (MQ=255)

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

                                                                                                                                                     aaTCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCg  >  1:2162585/1‑150 (MQ=255)

                                                                                                                                                     aaTCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTACTCTACCAGAATAAAATCACCAATTTACAGAGTCg  >  1:1058614/1‑150 (MQ=255)

GAAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCG  >  NC_000913/2704731‑2705029

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: `‑`