BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,305,400	A→C	5.7%	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,305,400	0	A	C	5.7%	860.2 / 25.6	252	W317G (TGG→GGG)	sbcD	exonuclease subunit SbcD
Reads supporting (aligned to +/- strand): ref base A (133/104); new base C (13/2); total (146/106)
Fisher's exact test for biased strand distribution p-value = 2.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

AAAGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGT  >  NC_002947/2305251‑2305549

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

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

   ggCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGATTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAggg                                                                                                                                                    >  2:4791085/1‑150 (MQ=255)

        ggTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGc                                                                                                                                               <  1:3350477/150‑1 (MQ=255)

            aTCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGtt                                                                                                                                           >  2:5182004/1‑150 (MQ=255)

             tCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCCTGGCTGGCGGTTg                                                                                                                                          >  1:2804511/1‑150 (MQ=255)

             tCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTg                                                                                                                                          <  2:5260653/150‑1 (MQ=255)

             tCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTg                                                                                                                                          <  1:4506664/150‑1 (MQ=255)

              cGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTgg                                                                                                                                         <  1:3920207/150‑1 (MQ=255)

              cGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTgg                                                                                                                                         <  2:1677026/150‑1 (MQ=255)

               gTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTggg                                                                                                                                        <  1:1033496/150‑1 (MQ=255)

                  ttCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTc                                                                                                                                     <  2:4941738/150‑1 (MQ=255)

                      gCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCa                                                                                                                                 >  2:1248642/1‑150 (MQ=255)

                       cGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCaa                                                                                                                                <  1:1440739/150‑1 (MQ=255)

                         ttgttgCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGc                                                                                                                              >  1:3154625/1‑150 (MQ=255)

                           gttgCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCGAGCAg                                                                                                                            >  2:610432/1‑150 (MQ=255)

                           gttgCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAg                                                                                                                            <  1:2485997/150‑1 (MQ=255)

                              gCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTc                                                                                                                         <  1:5423191/150‑1 (MQ=255)

                                  ccggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAGGCAGGTCGAcc                                                                                                                     >  1:2732724/1‑150 (MQ=255)

                                   cggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCa                                                                                                                    >  1:2041850/1‑150 (MQ=255)

                                     gcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACg                                                                                                                  >  1:3613671/1‑150 (MQ=255)

                                      cGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGTGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACgg                                                                                                                 <  2:833312/150‑1 (MQ=255)

                                      cGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACgg                                                                                                                 <  1:4939694/150‑1 (MQ=255)

                                        tACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGc                                                                                                               <  2:4433373/150‑1 (MQ=255)

                                          cTCGGCGCTGATACGATTCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAg                                                                                                             >  2:3183457/1‑150 (MQ=255)

                                          cTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAg                                                                                                             <  2:4578222/150‑1 (MQ=255)

                                           tCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAgg                                                                                                            >  1:728932/1‑150 (MQ=255)

                                             ggCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTc                                                                                                          >  1:3430839/1‑150 (MQ=255)

                                             ggCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTc                                                                                                          >  1:1529276/1‑150 (MQ=255)

                                              gcgcTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCg                                                                                                         >  1:304319/1‑150 (MQ=255)

                                              gcgcTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCg                                                                                                         >  2:1837954/1‑150 (MQ=255)

                                                gcTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGc                                                                                                       >  2:4271838/1‑150 (MQ=255)

                                                  tGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCa                                                                                                     >  2:1306335/1‑150 (MQ=255)

                                                    aTACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGNTCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGGTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGt                                                                                                   >  1:418234/1‑150 (MQ=255)

                                                    aTACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGt                                                                                                   <  2:4398031/150‑1 (MQ=255)

                                                     tACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGtt                                                                                                  >  2:858954/1‑150 (MQ=255)

                                                     tACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGtt                                                                                                  >  2:5503511/1‑150 (MQ=255)

                                                          aTCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgctgc                                                                                             <  2:3482056/150‑1 (MQ=255)

                                                          aTCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgctgc                                                                                             >  2:4510482/1‑150 (MQ=255)

                                                             aGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAAc                                                                                          <  1:518576/150‑1 (MQ=255)

                                                              gCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACa                                                                                         >  1:2922795/1‑150 (MQ=255)

                                                              gCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACa                                                                                         <  2:3745958/150‑1 (MQ=255)

                                                                 gTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCt                                                                                      >  1:2524292/1‑150 (MQ=255)

                                                                     cGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGc                                                                                                                              <  2:2836873/106‑1 (MQ=255)

                                                                       ccttgccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCa                                                                                >  2:243444/1‑150 (MQ=255)

                                                                             cttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTg                                                                          >  2:3915724/1‑150 (MQ=255)

                                                                             cttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTg                                                                          >  2:749914/1‑150 (MQ=255)

                                                                              ttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGt                                                                                                                                      <  1:1090604/89‑1 (MQ=255)

                                                                                 aaGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCgg                                                                      <  2:1982009/150‑1 (MQ=255)

                                                                                   gcgcGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCGACAGCTCGCCCAGCGGTGCCGGGc                                                                    >  1:2017738/1‑150 (MQ=255)

                                                                                   gcgcGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGc                                                                    <  1:1970483/150‑1 (MQ=255)

                                                                                       ggTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCAc                                                                <  2:3390283/150‑1 (MQ=255)

                                                                                         tttCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCAccc                                                              >  1:2701151/1‑150 (MQ=255)

                                                                                         tttCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCAccc                                                              >  1:1571260/1‑150 (MQ=255)

                                                                                           tCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGt                                                            >  1:4240347/1‑150 (MQ=255)

                                                                                             gATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTg                                                          <  1:5520860/150‑1 (MQ=255)

                                                                                              aTCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGc                                                         >  1:4398062/1‑150 (MQ=255)

                                                                                                 tgctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAgg                                                      <  2:1189319/150‑1 (MQ=255)

                                                                                                      gcgcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCAc                                                 >  1:1175676/1‑150 (MQ=255)

                                                                                                       cgcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACg                                                <  1:791400/150‑1 (MQ=255)

                                                                                                             tCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGt                                          <  2:2732590/150‑1 (MQ=255)

                                                                                                              cGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                         <  2:1615636/150‑1 (MQ=255)

                                                                                                                  ttGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTAc                                     >  2:3721264/1‑150 (MQ=255)

                                                                                                                     cGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACggg                                  <  1:3883169/150‑1 (MQ=255)

                                                                                                                     cGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACggg                                  <  1:2532127/150‑1 (MQ=255)

                                                                                                                       gCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGc                                >  1:4670952/1‑150 (MQ=255)

                                                                                                                         tcgtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGt                              <  2:1099776/150‑1 (MQ=255)

                                                                                                                            tcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGCCCCGTACCACGCGCCGTGGc                           <  2:5028026/150‑1 (MQ=255)

                                                                                                                             cgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGCCCCGTACCACGCGCCGTGGcc                          <  1:4266293/150‑2 (MQ=255)

                                                                                                                              gAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTc                         <  1:2026455/150‑1 (MQ=255)

                                                                                                                              gAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTc                         >  1:951319/1‑150 (MQ=255)

                                                                                                                              gAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTc                         >  1:2016784/1‑150 (MQ=255)

                                                                                                                                gCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGa                       >  2:4861497/1‑150 (MQ=255)

                                                                                                                                    cacCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCt                   >  2:270261/1‑150 (MQ=255)

                                                                                                                                     acCCGCACCTCCAGCCCGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTg                  >  1:1759033/1‑150 (MQ=255)

                                                                                                                                     acCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTg                  >  2:907416/1‑150 (MQ=255)

                                                                                                                                       cccaaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  1:5274958/5‑150 (MQ=255)

                                                                                                                                       cccaaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  1:2796597/5‑150 (MQ=255)

                                                                                                                                       cccaaCCGCCAGCCCTGGCTGGAGGTGCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  1:1897471/5‑150 (MQ=255)

                                                                                                                                       ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  1:3804653/1‑150 (MQ=255)

                                                                                                                                       ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  2:1923024/1‑150 (MQ=255)

                                                                                                                                       ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  2:2313785/1‑150 (MQ=255)

                                                                                                                                       ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                >  1:3998094/1‑150 (MQ=255)

                                                                                                                                        ccaaCCGCCAGCCCTGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAcc               >  1:2307524/4‑150 (MQ=255)

                                                                                                                                         gCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCa              <  1:3054079/150‑1 (MQ=255)

                                                                                                                                         caaCCGCCAGCCCTGGCTGGAGGTGCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTCCGGGGCGTGGCTCGACGCTGACCa              >  2:3085690/3‑150 (MQ=255)

                                                                                                                                         caaCCGCCAGCCCTGGCTGGAGGTGCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCa              >  1:4695787/3‑150 (MQ=255)

                                                                                                                                          aaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAg             >  2:2243839/2‑150 (MQ=255)

                                                                                                                                             cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCACCAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGcc          <  2:2017652/150‑1 (MQ=255)

                                                                                                                                             cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGcc          >  2:4822851/1‑150 (MQ=255)

                                                                                                                                             cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGcc          >  1:3132858/1‑150 (MQ=255)

                                                                                                                                             cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGcc          >  1:3530858/1‑150 (MQ=255)

                                                                                                                                              tCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCg         >  1:5601297/1‑150 (MQ=255)

                                                                                                                                                cagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGtt       >  2:2078451/1‑150 (MQ=255)

                                                                                                                                                cagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGtt       >  2:2678028/1‑150 (MQ=255)

                                                                                                                                                cagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGtt       >  1:2203927/1‑150 (MQ=255)

                                                                                                                                                cagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGtt       >  1:4864112/1‑150 (MQ=255)

                                                                                                                                                   ccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTAcc    <  2:3617484/150‑1 (MQ=255)

                                                                                                                                                     agGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGt  <  2:627037/150‑1 (MQ=255)

                                                                                                                                                     agGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGt  >  2:4084536/1‑150 (MQ=255)

AAAGGCCAGGTCATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGT  >  NC_002947/2305251‑2305549

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