BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,704,880	A→C	5.2%	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	5.2%	351.4 / 7.4	109	Y144D (TAC→GAC) ‡	lepB	signal peptidase I
Reads supporting (aligned to +/- strand): ref base A (55/48); new base C (3/3); total (58/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 110 total aligned reads displayed.

ATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCG  >  NC_000913/2704734‑2705029

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

  cGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGGCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCaa                                                                                                                                                  <  2:1115188/150‑1 (MQ=255)

   gCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAg                                                                                                                                                 >  2:908184/1‑150 (MQ=255)

    cTCGGTTCCACGTTTGAGTAGGTGACCTGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTTACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGc                                                                                                                                                <  2:254546/150‑1 (MQ=255)

       ggTTCCACGTTTGAGTAGGTGAGCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATGCCGGTTGAATCGTCAGCTCTTTTGAGACAGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATGAAGCttt                                                                                                                                             <  2:1073284/150‑1 (MQ=255)

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

         ttgcACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATACCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTTACTTTATCGCCCGGTAAACCCACCGCGGGCTTGATGTAATCAAGCTTTgg                                                                                                                                           <  2:289154/147‑1 (MQ=255)

         ttCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGAGCTTGATTACATCAAGcgcgcg                                                                                                                                           <  1:681047/150‑6 (MQ=255)

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

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

            cACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGAta                                                                                                                                        >  2:300908/1‑149 (MQ=255)

            cACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATc                                                                                                                                        >  2:840963/1‑150 (MQ=255)

              cGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCCTCCCGGTTGAATCGTCCGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCCCGCGCGGGATGTAATCAAGCTTTGGATCtt                                                                                                                                      >  1:345231/1‑150 (MQ=255)

                 ttGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCg                                                                                                                                   >  1:1013166/1‑150 (MQ=255)

                 ttGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCg                                                                                                                                   <  2:1229819/150‑1 (MQ=255)

                 ttGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCg                                                                                                                                   <  1:609872/150‑1 (MQ=255)

                 ttGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCg                                                                                                                                   >  1:1043415/1‑150 (MQ=255)

                  tGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                                                                                                                  <  1:801777/150‑1 (MQ=255)

                  agAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                                                                                                                  >  2:35776/2‑150 (MQ=255)

                    aGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGat                                                                                                                                <  1:338615/150‑1 (MQ=255)

                         gTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTa                                                                                                                           <  1:164796/150‑1 (MQ=255)

                           gACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTaaa                                                                                                                         >  2:1040274/1‑150 (MQ=255)

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

                               ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGTATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                     >  2:599679/1‑150 (MQ=255)

                               ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                     >  1:646960/1‑150 (MQ=255)

                                   gcgcGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCAACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCAcca                                                                                                                 >  2:544755/1‑148 (MQ=255)

                                      cgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGatat                                                                                                              <  1:35776/150‑1 (MQ=255)

                                        ttttCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcg                                                                                                            >  1:357048/1‑150 (MQ=255)

                                          ttCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgcc                                                                                                          >  2:1251921/1‑150 (MQ=255)

                                            cacaCGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgccgc                                                                                                        >  2:218829/1‑150 (MQ=255)

                                                 gCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTc                                                                                                   <  2:1352197/150‑1 (MQ=255)

                                                 gCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTc                                                                                                   <  1:1040274/150‑1 (MQ=255)

                                                  ccTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCg                                                                                                  <  1:337707/150‑1 (MQ=255)

                                                   cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCgg                                                                                                 >  1:803884/1‑150 (MQ=255)

                                                    tGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATTTTTAAAGACCACGATATCGCCGCGTTTCGGa                                                                                                >  1:1392126/1‑150 (MQ=255)

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

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

                                                       ccGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCCTAAGTGACTTTAACGCCCGGTAAACCAACCGCGCGCTTGATGTAATCAAGCTTTGCATCTTCCGCATATTTAAAGACCACGATATCGCCCCGATTCGGATGa                                                                                             >  2:491545/1‑150 (MQ=255)

                                                              accATCCACGTTTAATCATCAGATCATATTAGACAAGATAGAAAATGACTTTATAGAAAGGAAAACAAAACCCACGCTTAATGTAATAAAGCTTTGGATCTTACGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGttt                                                                                      >  2:712514/3‑150 (MQ=255)

                                                               gCATCCCGGTTGAATCGTCAGCTCTTTTGAGCCCGGATCGTAAGTGACTTTATCGCCCGGTACACCCCCCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCCCGATATCGCCGCGTTTCGGATGCCCGGGTTc                                                                                     >  2:1344777/1‑150 (MQ=255)

                                                               gCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTc                                                                                     <  2:1534905/150‑1 (MQ=255)

                                                                cATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCg                                                                                    >  2:325746/1‑150 (MQ=255)

                                                                  tCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCa                                                                                                                                                   <  2:503978/85‑1 (MQ=255)

                                                                  tCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCa                                                                                                                                                   >  1:503978/1‑85 (MQ=255)

                                                                      ggTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGc                                                                              <  1:840963/150‑1 (MQ=255)

                                                                          gAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTc                                                                                     <  1:669351/139‑1 (MQ=255)

                                                                          gAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTc                                                                                     >  2:669351/1‑139 (MQ=255)

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

                                                                              cGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTg                                                                      >  2:61967/1‑150 (MQ=255)

                                                                               gTCAGCTGTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTgg                                                                     <  1:490604/150‑1 (MQ=255)

                                                                               gTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTgg                                                                     >  2:1018955/1‑150 (MQ=255)

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

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

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

                                                                                 cAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTa                                                                   <  1:1349187/150‑1 (MQ=255)

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

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

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

                                                                                            agaCCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCtt                                                        >  1:1478909/1‑150 (MQ=255)

                                                                                             gaCCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCttt                                                       >  1:1121380/1‑150 (MQ=255)

                                                                                             gaCCGGATCGTAAGTGACTTTATCCCCCGGTAAACCAACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCttt                                                       >  2:530238/1‑150 (MQ=255)

                                                                                                cGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAAt                                                    >  1:289767/1‑150 (MQ=255)

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

                                                                                                 ggATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATg                                                   >  1:1505501/1‑150 (MQ=255)

                                                                                                 ggATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATg                                                   <  2:646960/150‑1 (MQ=255)

                                                                                                 ggATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCCCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTATTCTCGTAGATAGGATCTTTAATg                                                   >  1:1505504/1‑150 (MQ=255)

                                                                                                  gATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGc                                                  >  2:594426/1‑150 (MQ=255)

                                                                                                      gTAAGTGACTTTATCGCCCGGTAAACCCACAGCGCGCTTGATGTAATCAAGCTTTCGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATa                                              >  2:754119/1‑150 (MQ=255)

                                                                                                        aaGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAg                                            <  1:1372345/150‑1 (MQ=255)

                                                                                                        aaGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAg                                            >  2:434503/1‑150 (MQ=255)

                                                                                                         aGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGc                                           <  1:625840/150‑1 (MQ=255)

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

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

                                                                                                                  aTCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtctc                                  <  1:115147/150‑1 (MQ=255)

                                                                                                                  aTCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtctc                                  <  2:510473/150‑1 (MQ=255)

                                                                                                                  aTCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtctc                                  >  2:958048/1‑150 (MQ=255)

                                                                                                                   tCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTctct                                 >  1:1520619/1‑150 (MQ=255)

                                                                                                                   tCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTctct                                 >  1:1521729/1‑150 (MQ=255)

                                                                                                                    cGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTa                                <  2:721961/150‑1 (MQ=255)

                                                                                                                          gTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGaa                          <  1:61967/150‑1 (MQ=255)

                                                                                                                               cccACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATaa                                             >  2:1323857/1‑126 (MQ=255)

                                                                                                                               cccACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATaa                                             <  1:1323857/126‑1 (MQ=255)

                                                                                                                               cccACCGCGCGCTTCATGTAATCAAGCTTAGGATCTAACGGATATTTAAAGACAACAATATCGCCGCGTTTCGGATGACCTGTTTCGAACAGCGTTTTAAGTTAGAAAGGATCTTTAATGCCATAAGCAAACTTCACTACCAGAATaaaa                     >  2:391889/1‑150 (MQ=255)

                                                                                                                                ccACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAAt                    >  2:340648/1‑150 (MQ=255)

                                                                                                                                 cACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATTTTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATc                   <  2:1392126/150‑1 (MQ=255)

                                                                                                                                   tccaaagcTTGATTACATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGAGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCAc                 >  2:887937/7‑150 (MQ=255)

                                                                                                                                   ccGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCAc                 >  1:426444/1‑150 (MQ=255)

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

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

                                                                                                                                     caaagcTTGATTACATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCa                             >  2:745493/5‑136 (MQ=255)

                                                                                                                                     caaagcTTGATTACATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCa                             <  1:745493/132‑1 (MQ=255)

                                                                                                                                        cgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAAtt            >  1:499237/1‑150 (MQ=255)

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

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

                                                                                                                                             gATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACag       >  1:647127/1‑150 (MQ=255)

                                                                                                                                              aTGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACaga      >  1:1225975/1‑150 (MQ=255)

                                                                                                                                              aTGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACaga      >  2:1044434/1‑150 (MQ=255)

                                                                                                                                                 tAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTc   >  2:1455521/1‑150 (MQ=255)

                                                                                                                                                 tAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTc   >  2:687039/1‑150 (MQ=255)

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

ATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCG  >  NC_000913/2704734‑2705029

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