BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,638,545	A→T	5.1%	Q168L (CAG→CTG)	lepB →	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,638,545	0	A	T	5.1%	967.0 / 27.6	304	Q168L (CAG→CTG)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand): ref base A (162/120); new base T (3/13); total (171/133)
Fisher's exact test for biased strand distribution p-value = 3.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 314 total aligned reads displayed.

GGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTG  >  NC_002947/1638396‑1638694

ggTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCa                                                                                                                                                       <  2:3985180/150‑1 (MQ=255)

ggTCACCGACCTCGATGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCa                                                                                                                                                       >  2:2502576/1‑150 (MQ=255)

  tCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGc                                                                                                                                                     >  2:5286617/1‑150 (MQ=255)

    aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                   >  1:5868090/1‑150 (MQ=255)

    aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                   >  2:3765150/1‑150 (MQ=255)

    aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                   >  2:6050935/1‑150 (MQ=255)

    aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                   >  1:1617756/1‑150 (MQ=255)

         ggTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGGTGc                                                                                                                                              >  1:7228627/1‑150 (MQ=255)

          gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGcc                                                                                                                                             >  1:1448977/1‑150 (MQ=255)

          gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTgac                                                                                                                                             <  2:3199777/150‑3 (MQ=255)

          gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTgac                                                                                                                                             <  2:6262581/150‑3 (MQ=255)

          gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGACAACGGGCTGGCCGTTgac                                                                                                                                             <  2:3608804/150‑3 (MQ=255)

           tCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCg                                                                                                                                            <  2:2563400/150‑1 (MQ=255)

            cGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGa                                                                                                                                           <  2:6416047/150‑1 (MQ=255)

             ggTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaa                                                                                                                                          >  1:629209/1‑150 (MQ=255)

             ggTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCTACGGGCTGGCCGTTGACGaa                                                                                                                                          <  2:5398879/150‑1 (MQ=255)

             ggTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaa                                                                                                                                          <  1:4617030/150‑1 (MQ=255)

                   ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                              <  2:2025681/140‑1 (MQ=255)

                   ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtg                                                                                                                                    >  1:4063557/1‑150 (MQ=255)

                     gCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggt                                                                                                                                  <  2:1579804/150‑1 (MQ=255)

                       aaCGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCACCGGCCAGCCCGTTGCCGAACAACtggtgg                                                                                                                                <  2:4604822/150‑1 (MQ=255)

                       aaCGCGGGGATGTAATGGGGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtgg                                                                                                                                <  1:850732/150‑1 (MQ=255)

                        aCGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTggg                                                                                                                               <  2:3636701/150‑1 (MQ=255)

                          gcgGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCa                                                                                                                             >  2:4058760/1‑150 (MQ=255)

                                aTGTAATGGCGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGc                                                                                                                       >  2:6766748/1‑150 (MQ=255)

                                   tAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                                    <  1:1413527/150‑1 (MQ=255)

                                    aaTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCCGCGCGCCggg                                                                                                                   >  1:1029200/1‑150 (MQ=255)

                                    aaTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCATGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCggg                                                                                                                   >  1:524718/1‑150 (MQ=255)

                                       ggTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAc                                                                                                                >  1:620904/1‑150 (MQ=255)

                                           ttCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                            <  1:5891157/150‑1 (MQ=255)

                                           ttCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                            >  2:3403857/1‑150 (MQ=255)

                                             ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                                          <  1:4495901/150‑1 (MQ=255)

                                             ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                                          <  2:284932/150‑1 (MQ=255)

                                             ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                                          >  1:436165/1‑150 (MQ=255)

                                             ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                                          <  2:7056856/150‑1 (MQ=255)

                                             ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                            <  2:4101761/148‑1 (MQ=255)

                                               gCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                                        <  1:108256/150‑1 (MQ=255)

                                                  aCCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTg                                                                                                     <  2:1534441/150‑1 (MQ=255)

                                                       aGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAg                                                                                                <  1:5517117/150‑1 (MQ=255)

                                                            cccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                           <  2:599162/150‑1 (MQ=255)

                                                            cccGAATGTCAACGACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                           <  2:5027195/150‑1 (MQ=255)

                                                             ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCCGTGCTCAGCTGTAc                                                                                          >  1:3561185/1‑150 (MQ=255)

                                                             ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTAc                                                                                          >  2:326064/1‑150 (MQ=255)

                                                             ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTAc                                                                                          >  2:4784379/1‑150 (MQ=255)

                                                             ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACAGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTAc                                                                                          >  1:6724822/1‑150 (MQ=255)

                                                              cGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACAGGTGGGCAGCGAGCCGGGTACCCAGGGCAGTGCTCAGCAGTACa                                                                                         >  1:5052910/1‑150 (MQ=255)

                                                                aaTGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAg                                                                                       <  1:7190753/150‑1 (MQ=255)

                                                                    tCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGaaa                                                                                   >  1:7329992/1‑150 (MQ=255)

                                                                      aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAcc                                                                                                               >  2:6223580/1‑120 (MQ=255)

                                                                      aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                                 >  1:294206/1‑150 (MQ=255)

                                                                      aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                                 <  1:1339232/150‑1 (MQ=255)

                                                                      aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                                 <  1:2063731/150‑1 (MQ=255)

                                                                      aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                                 >  2:350132/1‑150 (MQ=255)

                                                                            aTCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAACCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgc                                                                           >  1:1088509/1‑150 (MQ=255)

                                                                            aTCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAACCCGTTGCCGAACAACTGGTGGGCCGCGAGCCGGGGACCCTGGGCAGTGCTCAGCTGTACAATGAAAATCACGgc                                                                           >  1:3718977/1‑150 (MQ=255)

                                                                                 gCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGCACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGc                                                                      >  2:6953594/1‑150 (MQ=255)

                                                                                 gCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGc                                                                      >  1:1916566/1‑150 (MQ=255)

                                                                                 gCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGc                                                                      <  2:4647802/150‑1 (MQ=255)

                                                                                   gtgtGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCCGCCCGTTGCCGAACAACTGGGGGGCAGCGAGCCGGGTACCCGGGGCAGTGCTCAGCTGTACAAGGAAAAGCTAGGCGCGGCCg                                                                    >  1:7020962/1‑150 (MQ=255)

                                                                                     gtgGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                            <  1:2697341/108‑1 (MQ=255)

                                                                                     gtgGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                            >  2:2697778/1‑108 (MQ=255)

                                                                                        gTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAacac                                                               >  2:7429901/1‑150 (MQ=255)

                                                                                        gTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAacac                                                               <  2:5489663/150‑1 (MQ=255)

                                                                                            gCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGa                                                           <  1:1572878/150‑1 (MQ=255)

                                                                                            gCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGa                                                           >  2:1086955/1‑150 (MQ=255)

                                                                                                 ccGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                                      <  1:6584496/150‑1 (MQ=255)

                                                                                                 ccGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                                      >  2:1357078/1‑150 (MQ=255)

                                                                                                 ccGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                                      >  1:5463540/1‑150 (MQ=255)

                                                                                                           tGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGa                                            <  2:1840881/150‑1 (MQ=255)

                                                                                                           tGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGa                                            <  1:3359238/150‑1 (MQ=255)

                                                                                                           tGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGa                                            <  1:3301123/150‑1 (MQ=255)

                                                                                                              tCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCCGCGAGCCGGGTACCCTGGGCAGTGCTCCGCTGTACCAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGcc                                         >  2:1236217/1‑150 (MQ=255)

                                                                                                              tCCGCTACACCAACGACAAGCGGATGTTCGTCAACGGCCAGCCCGTTGCCGAAAAACTGGTGGGCAGCGACCCGGGAACCCTGGGCAGTGCACAGCTGTACAAGGAAAACCACGGCGCGGCCGAACACCTCATCCGCAAGGAAATGAccc                                         >  2:5114078/1‑150 (MQ=255)

                                                                                                               ccGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCg                                        >  1:2846762/1‑150 (MQ=255)

                                                                                                                cGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGt                                       <  1:6754274/150‑1 (MQ=255)

                                                                                                                        cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCAt                               >  2:3054754/1‑150 (MQ=255)

                                                                                                                        cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCAt                               >  1:3071309/1‑150 (MQ=255)

                                                                                                                        cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCGCAGCGGTACAAGGAAAAGCTCGGCGCGGCCGCACACCTGATCCGCAAGGAAATGAGCCGTTATCGCag                               >  1:5182986/1‑149 (MQ=255)

                                                                                                                               aaGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGccgccg                        <  1:4448315/150‑1 (MQ=255)

                                                                                                                                 gCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGa                      >  1:2328919/1‑150 (MQ=255)

                                                                                                                                   ggCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGAcc                    >  2:1024485/1‑150 (MQ=255)

                                                                                                                                      tGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagc                 >  2:72049/1‑150 (MQ=255)

                                                                                                                                      tGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagc                 >  2:4381516/1‑150 (MQ=255)

                                                                                                                                        ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagcag               >  2:335251/1‑150 (MQ=255)

                                                                                                                                        ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagcag               >  2:2220809/1‑150 (MQ=255)

                                                                                                                                        ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagcag               >  1:6917979/1‑150 (MQ=255)

                                                                                                                                        ttCGGCAACGGGCTGGCCGTTGACGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagcag               >  1:3441419/1‑150 (MQ=255)

                                                                                                                                         tCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGt              <  1:5763972/150‑1 (MQ=255)

                                                                                                                                         tCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGt              >  1:2379172/1‑150 (MQ=255)

                                                                                                                                         tCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGt              <  2:587113/150‑1 (MQ=255)

                                                                                                                                             cAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAc          >  2:5719832/1‑150 (MQ=255)

                                                                                                                                              aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAcc         >  1:6031408/1‑150 (MQ=255)

                                                                                                                                              aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAcc         >  1:6482968/1‑150 (MQ=255)

                                                                                                                                              aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAcc         >  1:4284750/1‑150 (MQ=255)

                                                                                                                                              aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAcc         >  2:872754/1‑150 (MQ=255)

                                                                                                                                                    cAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCt   >  1:7048406/1‑150 (MQ=255)

                                                                                                                                                    cAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCt   >  2:7279768/1‑150 (MQ=255)

                                                                                                                                                    cAGCCCGTTGCCGAACAACTGGCGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGACATGAGCCGGTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCt   >  2:7016902/1‑150 (MQ=255)

                                                                                                                                                     aGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTg  >  2:6520159/1‑150 (MQ=255)

                                                                                                                                                     aGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTg  >  1:4867044/1‑150 (MQ=255)

GGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTG  >  NC_002947/1638396‑1638694

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