BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,300,416	C→A	8.6%	R314L (CGC→CTC)	sbcD ←	exonuclease subunit SbcD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,300,416	0	C	A	8.6%	352.9 / 18.7	116	R314L (CGC→CTC)	sbcD	exonuclease subunit SbcD
Reads supporting (aligned to +/- strand): ref base C (59/47); new base A (6/4); total (65/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 116 total aligned reads displayed.

CATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTC  >  NC_002947/2300270‑2300558

cATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCgg                                                                                                                                              >  4:213073/1‑149 (MQ=255)

    gTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTgg                                                                                                                                          >  8:262916/1‑149 (MQ=255)

    gTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCGCCAGCCCTGGCTGGAGGtgc                                                                                                                                           <  3:224874/148‑3 (MQ=255)

     tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTggg                                                                                                                                         <  3:205172/149‑1 (MQ=255)

     tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTggg                                                                                                                                         <  5:224166/149‑1 (MQ=255)

     tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTggg                                                                                                                                         <  8:245802/149‑1 (MQ=255)

     tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTggg                                                                                                                                         <  3:213073/149‑1 (MQ=255)

     tCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTggg                                                                                                                                         <  8:34987/149‑1 (MQ=255)

      cTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTgggg                                                                                                                                        >  8:105399/1‑149 (MQ=255)

            cGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCCAACCGCCAGCCCTGGCTGGAGGTGCGGGTGttgc                                                                                                                                  <  6:368696/149‑3 (MQ=38)

                gttgCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTg                                                                                                                                           >  7:3130/1‑136 (MQ=255)

                    cggccggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGtt                                                                                                                                            >  8:139392/1‑131 (MQ=255)

                    cggccggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGtt                                                                                                                                            <  7:139392/131‑1 (MQ=255)

                       ccggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGAc                                                                                                                       >  7:61923/1‑149 (MQ=255)

                        cggcGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCCGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGAc                                                                                                                       >  8:60623/1‑148 (MQ=255)

                           cGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACg                                                                                                                   <  4:282317/149‑1 (MQ=255)

                                          tACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGt                                                                                                    >  5:334496/1‑149 (MQ=255)

                                           aCGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCCGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGtt                                                                                                   >  7:138454/1‑149 (MQ=255)

                                             gAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGGCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgc                                                                                                 >  3:196643/1‑149 (MQ=255)

                                             gAATCAGCCGTACAGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgc                                                                                                 >  5:275288/1‑149 (MQ=255)

                                                tCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTtgctg                                                                                               >  4:133655/1‑148 (MQ=255)

                                                   gCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCaa                                                                                            >  5:204581/1‑148 (MQ=255)

                                                       tACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTc                                                                                                                                   <  3:177397/105‑1 (MQ=255)

                                                       tACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTc                                                                                                                                   >  4:177397/1‑105 (MQ=255)

                                                        aCCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAACTc                                                                                      >  6:64137/1‑149 (MQ=255)

                                                          cGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGc                                                                                    >  3:190083/1‑149 (MQ=255)

                                                               tgccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGc                                                                               >  5:350349/1‑149 (MQ=255)

                                                                gccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGTGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAg                                                                                <  6:275288/147‑1 (MQ=255)

                                                                gccttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCg                                                                              >  8:111506/1‑149 (MQ=255)

                                                                  cttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGccc                                                                                  <  2:297108/143‑1 (MQ=255)

                                                                  cttCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGccc                                                                                  >  1:297108/1‑143 (MQ=255)

                                                                    tCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGc                                                                          >  3:223022/1‑149 (MQ=255)

                                                                       aaggcgGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCgg                                                                       <  2:37585/146‑1 (MQ=255)

                                                                        gcgcGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCggg                                                                      >  6:281020/1‑149 (MQ=255)

                                                                         cgcgGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGc                                                                     >  5:9688/1‑149 (MQ=255)

                                                                          gcgGTTTCGAGCTGCTCGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCGGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGGCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGcc                                                                    <  3:42611/149‑1 (MQ=255)

                                                                              tttCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGACCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCAcc                                                                <  8:61923/149‑1 (MQ=255)

                                                                               ttCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGc                                                                     >  1:216181/1‑143 (MQ=255)

                                                                                 cGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGcc                                                                         >  6:238734/1‑137 (MQ=255)

                                                                                   aTCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTg                                                           >  5:324859/1‑149 (MQ=255)

                                                                                       gctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAg                                                        <  5:137970/148‑1 (MQ=255)

                                                                                       gctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAg                                                        >  6:137970/1‑148 (MQ=255)

                                                                                        ctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAggg                                                      >  1:311736/1‑149 (MQ=255)

                                                                                        ctgGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAgg                                                       <  7:14667/148‑1 (MQ=255)

                                                                                          ggCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCg                                                                              >  1:231027/1‑123 (MQ=255)

                                                                                          ggCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCg                                                                              <  2:231027/123‑1 (MQ=255)

                                                                                             gcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCgg                                                                                                         >  7:61764/1‑93 (MQ=255)

                                                                                             gcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCgg                                                                                                         <  8:61764/93‑1 (MQ=255)

                                                                                             gcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCgg                                                                       <  3:316081/127‑1 (MQ=255)

                                                                                             gcAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCgg                                                                       >  4:316081/1‑127 (MQ=255)

                                                                                                gatCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGgcg                                              >  5:172165/3‑149 (MQ=255)

                                                                                                 gTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACgggtg                                             <  7:111506/149‑4 (MQ=255)

                                                                                                  tCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGcgcg                                            >  6:74494/1‑149 (MQ=255)

                                                                                                    gggTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                          <  7:123597/149‑1 (MQ=255)

                                                                                                     ggTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTgg                                         >  1:95732/1‑149 (MQ=255)

                                                                                                      gTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGgc                                               <  4:344914/142‑1 (MQ=255)

                                                                                                      gTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGgc                                               >  3:344914/1‑142 (MQ=255)

                                                                                                       ttGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGcgcg                                            >  3:3825/1‑144 (MQ=255)

                                                                                                          cGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACgg                                    <  2:311736/149‑1 (MQ=255)

                                                                                                          cGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAACTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACgg                                    <  5:64137/149‑1 (MQ=255)

                                                                                                            gCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACgggg                                  >  4:37252/1‑149 (MQ=255)

                                                                                                             cTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGc                                 >  7:340526/1‑149 (MQ=255)

                                                                                                             cTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGc                                 <  5:74494/149‑1 (MQ=255)

                                                                                                               cgtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCAc                                                                                                                    <  2:140880/64‑1 (MQ=255)

                                                                                                               cgtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCAc                                                                                                                    >  1:140880/1‑64 (MQ=255)

                                                                                                                gtcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTg                              <  2:16279/149‑1 (MQ=255)

                                                                                                                 tcgAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGcgcg                                            >  6:127010/1‑134 (MQ=255)

                                                                                                                   gAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGATGCAGGGCCACGCCCCGTACCACGCGCCGTGGcc                           <  1:82013/149‑2 (MQ=255)

                                                                                                                    aGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTc                          <  4:190083/149‑1 (MQ=255)

                                                                                                                      cAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGa                        >  1:208008/1‑149 (MQ=255)

                                                                                                                      cAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGa                        >  2:265783/1‑149 (MQ=255)

                                                                                                                       aaCACCCGCACCTCCAGCCAGGGCGGGCGGTTGGGGGCTTCAAGCAGGTCGACCACGGGCAGGGCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGAc                       <  7:60623/149‑1 (MQ=255)

                                                                                                                        acacCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACg                      >  7:272178/1‑149 (MQ=255)

                                                                                                                         cacCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGc                     >  1:124380/1‑149 (MQ=255)

                                                                                                                           cccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTg                   <  3:133655/149‑1 (MQ=255)

                                                                                                                           cccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTg                   <  8:272178/149‑1 (MQ=255)

                                                                                                                            cccaaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                          >  3:238602/5‑125 (MQ=255)

                                                                                                                            cccaaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                          <  4:238602/121‑1 (MQ=255)

                                                                                                                            cccaaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGt                               >  4:79889/5‑136 (MQ=255)

                                                                                                                            cccaaCCGCCAGCCCTGGCTGGAGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGt                               <  3:79889/132‑1 (MQ=255)

                                                                                                                            cccaaCCGCCAGCCCTGGCTGGAGGTGCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGa                  >  2:82406/5‑149 (MQ=255)

                                                                                                                            ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGcc                                                                         >  1:37231/1‑94 (MQ=255)

                                                                                                                            ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGcc                                                                         <  2:37231/94‑1 (MQ=255)

                                                                                                                            ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTgg                             <  4:94493/138‑1 (MQ=255)

                                                                                                                            ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTgg                             >  3:94493/1‑138 (MQ=255)

                                                                                                                            ccGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGa                  >  5:365903/1‑149 (MQ=255)

                                                                                                                             ccaaCCGCCAGCCCTGGCTGGAGGTGCGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                 >  4:228032/4‑149 (MQ=255)

                                                                                                                             cGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                 <  8:340526/149‑1 (MQ=255)

                                                                                                                             cGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGAc                 >  3:211304/1‑149 (MQ=255)

                                                                                                                                aCCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAg              <  5:281020/149‑1 (MQ=255)

                                                                                                                                  cgccagccctGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCCCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGc            >  8:78078/3‑149 (MQ=255)

                                                                                                                                  cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTgg                             >  4:204016/1‑132 (MQ=255)

                                                                                                                                  cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGc            <  5:173333/149‑1 (MQ=255)

                                                                                                                                  cTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGc            >  2:10404/1‑149 (MQ=255)

                                                                                                                                    ccagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCg          >  5:69060/1‑149 (MQ=255)

                                                                                                                                     cagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGcc           <  6:204581/147‑1 (MQ=255)

                                                                                                                                     cagccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGt         >  4:220849/1‑149 (MQ=255)

                                                                                                                                       gccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTa       >  3:20219/1‑149 (MQ=255)

                                                                                                                                         cagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTAcc     <  6:9688/149‑1 (MQ=255)

                                                                                                                                             gCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTc  >  3:238525/1‑148 (MQ=255)

CATCGTCTTCGGCGTTGTTGCGGCCGGCGTACTCGGCGCTGATACGAATCAGCCGTACCGCCTTGCCTTCAAGCGCGGTTTCGATCTGCTGGCGCAGGTCGGGTTGCGGCTCGTCGAGCAACACCCGCACCTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTC  >  NC_002947/2300270‑2300558

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