BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,302,373	A→C	5.3%	W317G (TGG→GGG)	sbcD ←	exonuclease subunit SbcD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,302,373	0	A	C	5.3%	416.4 / 9.2	126	W317G (TGG→GGG)	sbcD	exonuclease subunit SbcD
Reads supporting (aligned to +/- strand): ref base A (59/60); new base C (5/2); total (64/62)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

AAAGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGT  >  NC_002947/2302224‑2302522

aaaGGCCATGTCATCGTCTTCGGCGATGTTGAGGACGGCGTACTCGGCGCTGATACGAAACAGCCGAACCGCCTTGCATTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAAAACCCGCACCAccagcca                                                                                                                                                        >  2:4178959/1‑150 (MQ=255)

aaaGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTccagcca                                                                                                                                                        >  2:4403274/1‑150 (MQ=255)

aaaGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTccagcca                                                                                                                                                        >  2:1217619/1‑150 (MQ=255)

 aaGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTccagccag                                                                                                                                                       <  1:4084386/150‑1 (MQ=255)

 aaGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTccagccag                                                                                                                                                       <  2:419177/150‑1 (MQ=255)

  aGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAgg                                                                                                                                                      <  2:3725796/150‑1 (MQ=255)

  aGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAgg                                                                                                                                                      <  1:2332197/150‑1 (MQ=255)

      cAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTg                                                                                                                                                  >  2:3487743/1‑150 (MQ=255)

      cAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTg                                                                                                                                                  <  1:66852/150‑1 (MQ=255)

          tCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCgg                                                                                                                                              >  1:1793442/1‑150 (MQ=255)

            aTCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGt‑t                                                                                                                                           >  2:4895329/1‑150 (MQ=255)

            aTCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGt‑t                                                                                                                                           <  1:1907584/150‑1 (MQ=255)

             tCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑Tg                                                                                                                                          <  1:619461/150‑1 (MQ=255)

             tCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCGCCAGCCCTGGCTGGAGGtgc                                                                                                                                           <  2:3371791/150‑3 (MQ=255)

              cGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑Tgg                                                                                                                                         <  1:2298026/150‑1 (MQ=255)

               gTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑Tggg                                                                                                                                        <  1:3668066/150‑1 (MQ=255)

               gTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑Tggg                                                                                                                                        <  1:2807838/150‑1 (MQ=255)

               gTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑Tggg                                                                                                                                        <  1:798677/150‑1 (MQ=255)

                tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGtgcggg                                                                                                                                        <  2:1424292/150‑6 (MQ=255)

                tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑Tgggg                                                                                                                                       >  1:1968730/1‑150 (MQ=255)

                 cTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCCAACCGCCAGCCCTGGCTGGAGGT‑GCGGGt                                                                                                                                      <  1:986547/150‑1 (MQ=255)

                  ttCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTc                                                                                                                                     >  1:1149032/1‑150 (MQ=255)

                  ttCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTc                                                                                                                                     >  2:272910/1‑150 (MQ=255)

                  ttCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGAGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTc                                                                                                                                     <  1:3389793/150‑1 (MQ=255)

                     ggCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTc                                                                                                                                  >  1:4005248/1‑150 (MQ=255)

                     ggCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTc                                                                                                                                  <  1:1374374/150‑1 (MQ=255)

                         ttgttgCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGc                                                                                                                              <  1:1019985/150‑1 (MQ=255)

                           gttgCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAg                                                                                                                            <  1:1217619/150‑1 (MQ=255)

                                 gccggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGAc                                                                                                                      >  2:595029/1‑150 (MQ=255)

                                    ggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCAc                                                                                                                   >  2:1108272/1‑150 (MQ=255)

                                     gcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCGCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACg                                                                                                                  <  2:3416947/150‑1 (MQ=255)

                                         aCTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCa                                                                                                              <  1:1453828/150‑1 (MQ=255)

                                           tCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAgg                                                                                                            >  2:4262938/1‑150 (MQ=255)

                                           tCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAgg                                                                                                            <  2:468340/150‑1 (MQ=255)

                                             ggCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTc                                                                                                          >  1:4984278/1‑150 (MQ=255)

                                               cgcTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCgg                                                                                                        >  2:2317347/1‑150 (MQ=255)

                                               cgcTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCgg                                                                                                        >  2:120031/1‑150 (MQ=255)

                                                gcTGATACGTATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGc                                                                                                       <  2:3266284/150‑1 (MQ=255)

                                                 cTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑GGGGGTCTTCAAGCAGGTCGACCACGGGCAGGGCGGcc                                                                                                      >  1:746803/1‑150 (MQ=255)

                                                      aCGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtg                                                                                                 >  1:4225809/1‑150 (MQ=255)

                                                      aCGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtg                                                                                                 >  2:4647915/1‑150 (MQ=255)

                                                        gAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgct                                                                                               <  2:1700515/150‑1 (MQ=255)

                                                        gAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgct                                                                                               <  2:1700476/150‑1 (MQ=255)

                                                        gAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgct                                                                                               <  1:4733076/150‑1 (MQ=255)

                                                        gAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgct                                                                                               >  1:2218471/1‑150 (MQ=255)

                                                        gAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgct                                                                                               <  1:2106703/150‑1 (MQ=255)

                                                           tCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCa                                                                                            >  1:1027476/1‑150 (MQ=255)

                                                                      gccttgccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGccc                                                                                 <  2:4455063/150‑1 (MQ=255)

                                                                      gccttgccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGccc                                                                                 >  2:3714787/1‑150 (MQ=255)

                                                                         ttgccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGc                                                                              >  2:4149871/1‑150 (MQ=255)

                                                                             cttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTg                                                                          >  1:858810/1‑150 (MQ=255)

                                                                             cttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTg                                                                          >  1:2467731/1‑150 (MQ=255)

                                                                              ttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGc                                                                         >  2:4117105/1‑150 (MQ=255)

                                                                              ttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGc                                                                         >  1:3230685/1‑150 (MQ=255)

                                                                                cAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCg                                                                       >  1:1412683/1‑150 (MQ=255)

                                                                                cAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCg                                                                       >  1:585846/1‑150 (MQ=255)

                                                                                         tttCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGccccacg                                                              >  2:259455/1‑146 (MQ=255)

                                                                                         tttCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCAccc                                                              >  2:3046884/1‑150 (MQ=255)

                                                                                             gATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCaa                                                                                           >  2:3572306/1‑117 (MQ=255)

                                                                                             gATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCaa                                                                                           <  1:3572306/117‑1 (MQ=255)

                                                                                                ctgctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAg                                                       <  1:520094/150‑1 (MQ=255)

                                                                                                  gctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAggg                                                     <  2:1371831/150‑1 (MQ=255)

                                                                                                       cgcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACg                                                >  1:4333035/1‑150 (MQ=255)

                                                                                                          aGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGgcg                                             >  2:604896/1‑150 (MQ=255)

                                                                                                            gtggggTTGCGGCTCGTCGAGCAACACCCGCACCTCCTGCGAGGGCTGGGGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCTCCCGTTGCAGGGCCACGGcgcg                                           >  2:3850454/4‑150 (MQ=255)

                                                                                                             tCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGt                                          <  2:2305071/150‑1 (MQ=255)

                                                                                                             tCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGt                                          >  2:5011483/1‑150 (MQ=255)

                                                                                                             tCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGt                                          <  2:3762487/150‑1 (MQ=255)

                                                                                                             tCGGGGTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGt                                          <  2:1434758/150‑1 (MQ=255)

                                                                                                              cGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                         <  1:1733468/150‑1 (MQ=255)

                                                                                                                  ttGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTAc                                     <  2:816827/150‑1 (MQ=255)

                                                                                                                    gCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACgg                                   <  1:4971501/150‑1 (MQ=255)

                                                                                                                     cGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACggg                                  <  1:4309508/150‑1 (MQ=255)

                                                                                                                      ggCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACgggg                                 >  2:358067/1‑150 (MQ=255)

                                                                                                                         tcgtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGt                              <  2:849785/150‑1 (MQ=255)

                                                                                                                          cgtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                         <  1:2678713/138‑1 (MQ=255)

                                                                                                                          cgtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                         >  2:2678713/1‑138 (MQ=255)

                                                                                                                           gtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTgg                            <  1:4895329/150‑1 (MQ=255)

                                                                                                                           gtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTgg                            >  2:3700041/1‑150 (MQ=255)

                                                                                                                              gAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTc                         >  2:4977672/1‑150 (MQ=255)

                                                                                                                              gAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTc                         >  1:1401789/1‑150 (MQ=255)

                                                                                                                                    cacCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCt                   <  1:1851974/150‑1 (MQ=255)

                                                                                                                                    cacCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCt                   <  2:3925542/150‑1 (MQ=255)

                                                                                                                                      cccGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGa                 >  2:2560651/1‑150 (MQ=255)

                                                                                                                                      cccGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGa                 >  2:2828427/1‑150 (MQ=255)

                                                                                                                                       ccGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                <  1:1839909/150‑1 (MQ=255)

                                                                                                                                        ccaaCCGCCAGCCCTGGCTGGAGGTGCGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  1:897763/4‑150 (MQ=255)

                                                                                                                                        ccaaCCGCCAGCCCTGGCTGGAGGT‑GCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAcc               >  1:2608425/4‑150 (MQ=255)

                                                                                                                                        ccaaCCGCCAGCCCTGGCTGGAGGT‑GCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAcc               >  1:2608301/4‑150 (MQ=255)

                                                                                                                                        cGCACCTCCAGCCAGGGCTGGCGGG‑GGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCTGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAcc               >  1:1402901/1‑150 (MQ=255)

                                                                                                                                         gCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCa              >  1:4788503/1‑150 (MQ=255)

                                                                                                                                         gCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCa              <  2:4225809/150‑1 (MQ=255)

                                                                                                                                          aaCCGCCAGCCCTGGCTGGAGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAg             >  1:2797606/2‑150 (MQ=255)

                                                                                                                                               ccagccagGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGt        <  1:4403274/150‑1 (MQ=255)

                                                                                                                                                cagccagGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGtt       >  2:1276330/1‑150 (MQ=255)

                                                                                                                                                   ccagGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTAcc    >  1:1105102/1‑150 (MQ=255)

                                                                                                                                                   ccagGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTAcc    >  2:3605783/1‑150 (MQ=255)

                                                                                                                                                    cagGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCg   <  1:3870272/150‑1 (MQ=255)

                                                                                                                                                     agGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGt  >  2:4081180/1‑150 (MQ=255)

                                                                                                                                                     agGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGt  >  2:1255788/1‑150 (MQ=255)

AAAGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGT‑TGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGT  >  NC_002947/2302224‑2302522

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