BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,361,719	C→G	6.6%	G720R (GGC→CGC)	PP_2076 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,361,719	0	C	G	6.6%	377.0 / 11.2	120	G720R (GGC→CGC)	PP_2076	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (57/55); new base G (2/6); total (59/61)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 124 total aligned reads displayed.

CCACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTC  >  NC_002947/2361571‑2361866

ccACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACAAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTCCTCCTCGTCCGGTCACACTGTTCGCGcc                                                                                                                                                      >  2:261226/1‑149 (MQ=255)

   ctctGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTa                                                                                                                                                   >  8:71993/1‑149 (MQ=255)

        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                                                                                                                                             <  7:138977/145‑1 (MQ=255)

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

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

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

                    aTGCTCAATCGAGTCAATGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGt                                                                                                                                  >  7:223715/1‑149 (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)

                                            ccaccaCGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgcg                                                                                                                                            >  8:304328/1‑115 (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)

                                                     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)

                                                            aaCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCg                                                                                           >  4:120583/1‑148 (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)

                                                                               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)

                                                                                        cccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                         >  3:306118/1‑106 (MQ=40)

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

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

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

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

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

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

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

                                                                                              gTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCAc                                                        <  3:120583/148‑1 (MQ=39)

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

                                                                                                  cAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCcgc                                                    <  3:296395/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)

                                                                                                                 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)

                                                                                                                          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                    <  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)

                                                                                                                                     tCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGcc                <  3:246354/149‑1 (MQ=39)

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

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

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

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

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

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

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

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

                                                                                                                                                gcgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATccc     >  2:332770/1‑149 (MQ=39)

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

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

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

CCACTCTGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAG‑TCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTC  >  NC_002947/2361571‑2361866

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