BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,361,724	C→G	6.0%	R718P (CGG→CCG)	PP_2076 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,361,724	0	C	G	6.0%	345.3 / 7.4	113	R718P (CGG→CCG)	PP_2076	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (52/54); new base G (2/5); total (54/59)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 124 total aligned reads displayed.

GCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCCA  >  NC_002947/2361579‑2361868

gctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTAcc                                                                                                                                                   <  1:230408/146‑1 (MQ=255)

 ctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                               <  1:117670/149‑1 (MQ=255)

 ctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                               >  2:261513/1‑149 (MQ=255)

  tgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                               >  5:117058/1‑148 (MQ=255)

     ttttGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                               >  8:138977/1‑145 (MQ=255)

     ttttGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                               <  7:138977/145‑1 (MQ=255)

         gTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCGGTACGGCGCGaacagt                                                                                                                                        <  7:141864/148‑3 (MQ=255)

            aTGCTCAATCGAGTCAATGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGt                                                                                                                                    >  7:223715/1‑149 (MQ=255)

            aTGCTCAATCGAGTCAATGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGt                                                                                                                                    <  8:223715/149‑1 (MQ=255)

                tCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCca                                                                                                                                 >  1:16146/1‑148 (MQ=255)

                  aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTAcc                                                                                                                                                   <  6:90277/128‑1 (MQ=255)

                  aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTAcc                                                                                                                                                   >  5:90277/1‑128 (MQ=255)

                  aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCcaca                                                                                                                               >  5:154969/1‑148 (MQ=255)

                  aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTAAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCCCCCTGTTCGCGCCGTACCGCGCGACCACCTGTTCcaca                                                                                                                               >  6:71409/1‑148 (MQ=255)

                      gAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACt                                                                                                                          <  2:16146/149‑1 (MQ=255)

                       aGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTg                                                                                                                         <  6:117058/149‑1 (MQ=255)

                           aGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTc                                                                                                                     >  4:346111/1‑149 (MQ=255)

                            gTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTc                                                                                                                     >  6:323348/1‑148 (MQ=255)

                             tGAAATCCCACCACGCGAGTACCAACCCACCGCTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTACTCGTCCGGTCACACTGTTCGCGACCAACCCCGCGACAACCTCGTCCACAAAACTAGCTcc                                                                                                                    >  6:362666/1‑148 (MQ=17)

                              gAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCaa                                                                                                                  >  8:10041/1‑149 (MQ=255)

                                  tCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTc                                                                                                              <  3:346111/149‑1 (MQ=255)

                                    ccaccaCGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCg                                                                                                             >  7:179894/1‑148 (MQ=255)

                                     caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAAc                                                                                                                 >  3:330039/1‑143 (MQ=255)

                                     caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAAc                                                                                                                 <  4:330039/143‑1 (MQ=255)

                                     caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                           <  5:323348/149‑1 (MQ=255)

                                     caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                           <  7:9314/149‑1 (MQ=255)

                                      accacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAg                                                                                                          >  2:98197/1‑149 (MQ=255)

                                       ccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGc                                                                                                         <  7:112948/149‑1 (MQ=255)

                                        cacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGcc                                                                                                        >  4:107309/1‑149 (MQ=255)

                                             gAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCccgccg                                                                                                   >  6:175344/1‑149 (MQ=255)

                                                  ccAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTcc                                                                                                                                  <  3:264949/113‑1 (MQ=255)

                                                  ccAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTcc                                                                                                                                  >  4:264949/1‑113 (MQ=255)

                                                   cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTc                                                                                                                     >  6:40970/1‑125 (MQ=255)

                                                   cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTc                                                                                                                     <  5:40970/125‑1 (MQ=255)

                                                   cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAAc                                                                                                                 >  6:361958/1‑129 (MQ=255)

                                                   cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCccgc                                                                                                     >  6:319090/1‑141 (MQ=255)

                                                   cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCccgc                                                                                                     <  5:319090/141‑1 (MQ=255)

                                                   cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCg                                                                                             >  4:3824/1‑149 (MQ=255)

                                                    aaCCCCCCGGTGCCAGCTAAACTGGCGCCCCATACTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCg                                                                                             >  6:143314/1‑148 (MQ=255)

                                                     aCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                           <  1:98197/149‑1 (MQ=255)

                                                     aCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                           <  7:10041/149‑1 (MQ=255)

                                                      ccccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAAACCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTACAACTTCGAACCCCGCCGCAGCCGTac                                                                                          >  8:194755/1‑148 (MQ=255)

                                                         cccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCa                                                                                      >  7:5084/1‑149 (MQ=255)

                                                          ccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCAt                                                                                     <  3:3824/149‑1 (MQ=255)

                                                            ggTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATcg                                                                                   <  1:261226/149‑1 (MQ=40)

                                                                ccAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGgc                                                                               <  3:107309/149‑1 (MQ=38)

                                                                 cAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTc                                                                                                                                   >  3:303752/1‑97 (MQ=37)

                                                                       aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCa                                                                          <  7:194755/147‑1 (MQ=37)

                                                                       aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCATCTTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCAACTGGTCCACATAACTGGCTCCACCTTCGAAGCACGCCGAAGCCGTAG‑TCATcgccgctggccgt                                                                        >  1:123890/1‑142 (MQ=11)

                                                                          tGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGt                                                                     >  7:312062/1‑149 (MQ=39)

                                                                               ccccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGgccgc                                                                            >  3:177091/1‑137 (MQ=37)

                                                                               ccccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGgccgc                                                                            <  4:177091/137‑1 (MQ=37)

                                                                                cccATTGTTTCTAAGGCCCTAGAGCGTCCGCTGCTACTCGTCCGGTCACACTGTTCGGGCCGTACCGCGCGACCACTTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCATCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGt                                                                <  5:362666/148‑1 (MQ=14)

                                                                                cccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                           <  4:306118/106‑1 (MQ=40)

                                                                                cccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTc                                                               <  3:297091/149‑1 (MQ=39)

                                                                                cccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTc                                                               <  6:205012/149‑1 (MQ=39)

                                                                                cccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTc                                                               >  4:296395/1‑149 (MQ=39)

                                                                                  cATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGCCCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCt                                                              >  2:204049/1‑148 (MQ=255)

                                                                                   aTAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCtt                                                             >  5:287045/1‑148 (MQ=39)

                                                                                    tAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCgg                                                                                <  8:113562/128‑1 (MQ=38)

                                                                                    tAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCgg                                                                                >  7:113562/1‑128 (MQ=38)

                                                                                          cAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCcgc                                                      <  3:296395/148‑1 (MQ=39)

                                                                                          cAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCcgc                                                      <  8:179894/148‑1 (MQ=39)

                                                                                           aaaCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCcgcg                                                     >  8:59798/1‑148 (MQ=39)

                                                                                            aaagCCCTAGAGCGGCCGCTGCTACTAGTACGGTCACACTGTTCTCGCCGTACCGCGCGACCACCTGGTCCACAAAACTGGCTCCAACTTCGAAGCCCGCCCCAGCCGTAG‑TAATCGCCGCCGCCAGGAGTTc                                                                   >  5:231742/4‑133 (MQ=14)

                                                                                            aaCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTc                                                                   <  6:231742/133‑1 (MQ=40)

                                                                                                  cTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTg                                             >  4:208440/1‑149 (MQ=39)

                                                                                                     gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGt                                            <  8:5084/147‑1 (MQ=38)

                                                                                                       gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAg                                                                         >  6:309529/1‑116 (MQ=35)

                                                                                                       gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAg                                                                         <  5:309529/116‑1 (MQ=35)

                                                                                                       gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCa                                                           <  1:29578/130‑1 (MQ=38)

                                                                                                       gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCa                                                           >  2:29578/1‑130 (MQ=38)

                                                                                                         ggCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGc                                       <  1:261513/149‑1 (MQ=255)

                                                                                                             gctgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                           >  8:314814/1‑93 (MQ=37)

                                                                                                             gctgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                           <  7:314814/93‑1 (MQ=37)

                                                                                                              ctgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCccg                                                                                                      >  1:97052/1‑81 (MQ=38)

                                                                                                              ctgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCccg                                                                                                      <  2:97052/81‑1 (MQ=38)

                                                                                                                 ctcctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcgc                              >  1:132196/1‑149 (MQ=39)

                                                                                                                 ctcctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcgc                              <  8:312062/149‑1 (MQ=39)

                                                                                                                  tcctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcgcg                             >  3:189610/1‑149 (MQ=39)

                                                                                                                    ctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCt                           >  4:246354/1‑149 (MQ=39)

                                                                                                                       gTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCtt                          <  5:175344/147‑1 (MQ=39)

                                                                                                                          cGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                      <  3:229560/148‑1 (MQ=39)

                                                                                                                          cGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                      <  5:71409/148‑1 (MQ=39)

                                                                                                                           ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcg                               >  5:84455/1‑138 (MQ=40)

                                                                                                                           ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcg                               <  6:84455/138‑1 (MQ=40)

                                                                                                                           ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAgg                    <  1:204049/149‑1 (MQ=39)

                                                                                                                               acacTGTTCGCGCCGTACCGCGCGACCACCTTGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCTCAGCCGTAG‑TCATCGCGGACTCCAGGAGTTCGGTCTTCAACCGAGCATCACTGTCATGCAGTTCGCGAGCTTTGAAGGACa                 >  4:199316/1‑148 (MQ=14)

                                                                                                                                cacTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATc               <  7:59798/149‑1 (MQ=39)

                                                                                                                                 ggtggtCGCGCGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCAt                <  2:89128/142‑1 (MQ=255)

                                                                                                                                 ggtggtCGCGCGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCAt                >  1:89128/6‑147 (MQ=255)

                                                                                                                                 acTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAgg                    >  4:26847/1‑143 (MQ=39)

                                                                                                                                 acTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAgg                    <  3:26847/143‑1 (MQ=39)

                                                                                                                                 acTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCa                 >  5:36342/1‑146 (MQ=39)

                                                                                                                                        gcgcGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGACAGTTCGCGCGCTTTGAAGGCCATCAGCATccc       <  1:135643/149‑1 (MQ=255)

                                                                                                                                           ccGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTc    <  5:143314/149‑1 (MQ=39)

                                                                                                                                           ccGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTc    <  3:208440/149‑1 (MQ=39)

                                                                                                                                            cGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCTAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATc         >  7:267444/1‑143 (MQ=255)

                                                                                                                                            agTAACGCGCGACCACCTGGTACACATAACTGGCTCCTACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCTAGGAGTTCGGTCTTCAGCCGCGCATCTCTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCAtg         <  8:267444/142‑2 (MQ=14)

                                                                                                                                             gTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCca  <  4:189610/149‑1 (MQ=39)

GCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCCA  >  NC_002947/2361579‑2361868

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