BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,488,073	T→A	6.8%	Y86F (TAC→TTC)	degS ←	serine endoprotease DegS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,488,073	0	T	A	6.8%	1125.3 / 48.6	339	Y86F (TAC→TTC)	degS	serine endoprotease DegS
Reads supporting (aligned to +/- strand): ref base T (155/160); new base A (9/14); total (165/174)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 359 total aligned reads displayed.

TCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATCAGC  >  NC_002947/1487924‑1488222

tCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGgtgt                                                                                                                                                       >  1:7699465/1‑150 (MQ=255)

tCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGgtgt                                                                                                                                                       >  1:5294420/1‑150 (MQ=255)

   gcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCAAGGTGGTGAACa                                                                                                                                                    <  1:5739577/150‑1 (MQ=255)

    cagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAg                                                                                                                                                   <  2:1198870/150‑1 (MQ=255)

    cagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCAAGGTGGTGAACaa                                                                                                                                                   <  1:3856401/150‑2 (MQ=255)

      ggTAGCCCTGTGGGCTCATTATTACCGCCGACCCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCATGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACaaaa                                                                                                                                                 <  1:4057954/150‑4 (MQ=255)

      ggTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGt                                                                                                                                                 >  1:5247333/1‑150 (MQ=255)

      ggTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGGTGGTGAACaaaa                                                                                                                                                 <  1:794739/150‑4 (MQ=255)

      ggTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGGTGGTGAACaaaa                                                                                                                                                 <  1:472534/150‑4 (MQ=255)

        tAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGACCTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTg                                                                                                                                               <  1:3431393/150‑1 (MQ=255)

          gCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGc                                                                                                                                             <  1:224668/150‑1 (MQ=255)

          gCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGc                                                                                                                                             <  1:2599547/150‑1 (MQ=255)

          gCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGc                                                                                                                                             <  2:8425014/150‑1 (MQ=255)

          gCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGc                                                                                                                                             <  1:7217434/150‑1 (MQ=255)

          gCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGc                                                                                                                                             <  2:8823263/150‑1 (MQ=255)

             ctctGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCAc                                                                                                                                          >  2:328491/1‑150 (MQ=255)

              tctGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCAcc                                                                                                                                         >  1:2956079/1‑150 (MQ=255)

                     tCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCgggg                                                                                                                                  >  1:1454402/1‑150 (MQ=255)

                     tCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCgggg                                                                                                                                  >  1:6057989/1‑150 (MQ=255)

                     tCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCcggg                                                                                                                                  >  1:1352501/1‑150 (MQ=255)

                      cATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGggc                                                                                                                                 <  1:1846762/150‑1 (MQ=255)

                          aTTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGcc                                                                                                                             >  2:2656431/1‑150 (MQ=255)

                           ttACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCa                                                                                                                            <  1:321537/150‑1 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCGCCCCGGCGGTGGCCaac                                                                                                                          <  1:6475596/150‑3 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                          <  1:4926955/150‑1 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                          <  2:6693146/150‑1 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                          <  2:8637530/150‑1 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                          <  1:2025833/150‑1 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                          <  1:8030713/150‑1 (MQ=255)

                             aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                          <  2:7897009/150‑1 (MQ=255)

                                gccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTc                                                                                                                       >  1:7544971/1‑150 (MQ=255)

                                 ccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCa                                                                                                                      <  1:5788366/150‑1 (MQ=255)

                                    aaCCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcgg                                                                                                                   >  2:1771607/1‑150 (MQ=255)

                                       ccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCC‑GGGCGGCCAGGGTCAcggcgtc                                                                                                               >  2:4749919/1‑150 (MQ=255)

                                       ccAGGCTTGATTACCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCTGGTCTACGAACAGCGTGTGGGCGCTTTTGTTCACCAACTTGGAGGTGTAGAGGTTGGACACCGCCGGGGAGGCCAGGGTCAcggcgt                                                                                                                >  2:2435251/1‑150 (MQ=255)

                                         aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcg                                                                                                              <  1:5468468/150‑1 (MQ=255)

                                          ggCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTCCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcgg                                                                                                             >  1:3727252/1‑150 (MQ=255)

                                             ttGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcggcgt                                                                                                          >  2:4824750/1‑150 (MQ=255)

                                               gATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTAc                                                                                                        >  2:6135079/1‑150 (MQ=255)

                                                 ttCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGa                                                                                                      <  1:4440743/150‑1 (MQ=255)

                                                 ttCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGa                                                                                                      >  2:6393388/1‑150 (MQ=255)

                                                    ccAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGacac                                                                                                   <  2:2955985/150‑1 (MQ=255)

                                                     cAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACAcc                                                                                                  >  2:2189912/1‑150 (MQ=255)

                                                     cAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACAcc                                                                                                  <  2:4244934/150‑1 (MQ=255)

                                                      aGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCg                                                                                                 <  1:3510003/150‑1 (MQ=255)

                                                            gcTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACAGGGCCtt                                                                                           <  1:3546143/150‑1 (MQ=255)

                                                             cTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTg                                                                                          <  1:8455550/150‑1 (MQ=255)

                                                             cTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTg                                                                                          <  2:731764/150‑1 (MQ=255)

                                                               gTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCa                                                                                        <  2:6940267/150‑1 (MQ=255)

                                                                    ggCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATc                                                                                   <  2:4536156/150‑1 (MQ=255)

                                                                       aGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                                >  1:2702762/1‑150 (MQ=255)

                                                                         gTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGt                                                                              <  2:5959629/150‑1 (MQ=255)

                                                                             tCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTc                                                                          >  1:9074260/1‑150 (MQ=255)

                                                                                ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                                >  1:210124/1‑141 (MQ=255)

                                                                                ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                                <  2:210124/141‑1 (MQ=255)

                                                                                 cGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCtgtg                                                                      <  1:181182/150‑1 (MQ=255)

                                                                                 cGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCtgtg                                                                      <  1:4091497/150‑1 (MQ=255)

                                                                                      aagCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCt                                                                 >  2:1972061/1‑150 (MQ=255)

                                                                                      aagCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCt                                                                 >  1:6874257/1‑150 (MQ=255)

                                                                                        gcggcgGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtg                                                               >  2:2583638/1‑150 (MQ=255)

                                                                                        gcggcgGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCGTTTGTTCACCACCTTGGTGGTGGACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCtgg                                                               >  2:9066029/1‑148 (MQ=255)

                                                                                             ggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCa                                                          >  1:6354028/1‑150 (MQ=255)

                                                                                                  tGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTg                                                     >  2:8173758/1‑150 (MQ=255)

                                                                                                  tGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTg                                                     <  1:8334750/150‑1 (MQ=255)

                                                                                                    cGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGAc                                                   <  2:4828649/150‑1 (MQ=255)

                                                                                                     gggTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACg                                                  >  2:8500453/1‑150 (MQ=255)

                                                                                                         cTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCct                                              <  2:6735632/150‑1 (MQ=255)

                                                                                                             gAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                                          >  2:1573445/1‑150 (MQ=255)

                                                                                                             gAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                                          >  2:3634560/1‑150 (MQ=255)

                                                                                                             gAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                                          >  1:7022819/1‑150 (MQ=255)

                                                                                                                 aGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGc                                      >  1:6217207/1‑150 (MQ=255)

                                                                                                                         ggCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAc                              <  2:5831881/150‑1 (MQ=255)

                                                                                                                          gcgcTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACCCCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAcc                             >  1:1168192/1‑150 (MQ=255)

                                                                                                                            gcTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCa                           >  1:8895258/1‑150 (MQ=255)

                                                                                                                              ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt                         >  2:302364/1‑150 (MQ=255)

                                                                                                                              ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt                         >  2:7267260/1‑150 (MQ=255)

                                                                                                                              ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt                         >  1:5098326/1‑150 (MQ=255)

                                                                                                                                 tGTACACCACCAAGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                         >  2:5193638/1‑115 (MQ=255)

                                                                                                                                 tGTACACCACCAAGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCg                      >  1:4249847/1‑150 (MQ=255)

                                                                                                                                 tGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                                >  1:6696509/1‑92 (MQ=39)

                                                                                                                                  gtacaccacCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCgg                     >  2:7472338/4‑150 (MQ=255)

                                                                                                                                      accaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa                 <  2:5579321/150‑1 (MQ=255)

                                                                                                                                      accaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa                 <  1:1643250/150‑1 (MQ=255)

                                                                                                                                      accaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa                 >  2:7691204/1‑150 (MQ=255)

                                                                                                                                       ccaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTa                                                                             >  2:4344918/1‑89 (MQ=255)

                                                                                                                                        caccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACg               >  2:541368/1‑150 (MQ=255)

                                                                                                                                        caccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACg               >  2:312892/1‑150 (MQ=255)

                                                                                                                                        caccTTGGTGGTGTACAGGATGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGTCCGACCCATTGCGGAAAACg               >  2:3831448/1‑150 (MQ=255)

                                                                                                                                         accTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGc              >  1:8746670/1‑150 (MQ=255)

                                                                                                                                          ccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt             <  2:4025214/150‑1 (MQ=255)

                                                                                                                                          ccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt             <  1:2812509/150‑1 (MQ=255)

                                                                                                                                          ccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt             <  2:1134160/150‑1 (MQ=255)

                                                                                                                                          ccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt             <  2:3728181/150‑1 (MQ=255)

                                                                                                                                          ccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt             <  2:5723161/150‑1 (MQ=255)

                                                                                                                                           cTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTg            <  1:8016128/150‑1 (MQ=255)

                                                                                                                                           cTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTg            <  2:1722874/150‑1 (MQ=255)

                                                                                                                                                 ggtgTTCAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGgatgat      <  2:8840331/150‑1 (MQ=255)

                                                                                                                                                  gtgtACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATc     >  1:598013/1‑150 (MQ=255)

                                                                                                                                                   tgtACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATCa    <  2:1356688/150‑1 (MQ=255)

                                                                                                                                                     tACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATCAGc  <  1:8534500/150‑1 (MQ=255)

TCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGATGATCAGC  >  NC_002947/1487924‑1488222

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