BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,346,311	C→G	7.6%	Q261H (CAG→CAC)	PP_2941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,346,311	0	C	G	7.6%	457.3 / 20.3	142	Q261H (CAG→CAC)	PP_2941	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (81/49); new base G (5/6); total (86/56)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 143 total aligned reads displayed.

GCCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACT  >  NC_002947/3346162‑3346459

gCCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAc                                                                                                                                                      <  1:2620868/150‑1 (MQ=255)

   aTTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGa                                                                                                                                                   >  1:5293773/1‑150 (MQ=255)

          aagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                            <  1:1872209/150‑1 (MQ=255)

             ggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTa                                                                                                                                         >  1:2459436/1‑150 (MQ=255)

              gaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAt                                                                                                                                        <  1:3254202/150‑1 (MQ=255)

              gaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAt                                                                                                                                        <  2:3588038/150‑1 (MQ=255)

              gaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAt                                                                                                                                        <  2:1721858/150‑1 (MQ=255)

               aaggTTATGGCTGAACAGTACGTGGAAAGTCAGACAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                       <  2:3294845/150‑1 (MQ=255)

               aaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                       >  1:4071154/1‑150 (MQ=255)

                aggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtca                                                                                                                                      >  1:5433343/1‑150 (MQ=255)

                          tGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGcc                                                                                                                            <  2:2622768/150‑1 (MQ=255)

                            aaCAGTACTTGGAATGTCAGAAAAGCGGAGCTTGTGCTGCAAGGTGATTTCGCCGAGGTCTGATCGGACCCACCGAAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGATTCACTGAAGCGGCTCTATCAGCTCAGCGCCtt                                                                                                                          <  1:252688/150‑1 (MQ=255)

                            aaCAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCtt                                                                                                                          <  1:932549/150‑1 (MQ=255)

                                gTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGAtt                                                                                                                      <  2:3133773/150‑1 (MQ=255)

                                gTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGAtt                                                                                                                      >  2:1009575/1‑150 (MQ=255)

                                     tGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtca                                                                                                                                      <  2:3984920/129‑1 (MQ=255)

                                     tGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtca                                                                                                                                      >  1:3984920/1‑129 (MQ=255)

                                        aaaGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCCCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtcagct                                                                                                                                   <  1:4969619/129‑1 (MQ=255)

                                        aaaGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCCCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtcagct                                                                                                                                   >  2:4969619/1‑129 (MQ=255)

                                          aGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTg                                                                                                            >  2:3473937/1‑150 (MQ=255)

                                           gTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGCGTTCGCCGAGGTCTGATCGGACCCCCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                           <  1:104409/150‑1 (MQ=255)

                                            tCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGc                                                                                                          >  2:515995/1‑150 (MQ=255)

                                            tCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGc                                                                                                          <  1:878685/150‑1 (MQ=255)

                                             cAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCa                                                                                                         >  1:5406164/1‑150 (MQ=255)

                                             cAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCa                                                                                                         >  1:15448/1‑150 (MQ=255)

                                              aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCaa                                                                                                        >  2:4910762/1‑150 (MQ=255)

                                                   aGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCt                                                                                                   <  1:400029/150‑1 (MQ=255)

                                                    gCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtg                                                                                                  <  2:5335805/150‑1 (MQ=255)

                                                       gAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctct                                                                                                                                          <  1:3456989/107‑1 (MQ=255)

                                                       gAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctct                                                                                                                                          >  2:3456989/1‑107 (MQ=255)

                                                         gCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                            <  2:1542329/103‑1 (MQ=255)

                                                         gCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                            >  1:4288237/1‑103 (MQ=255)

                                                           ttGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCAc                                                                                           >  2:3340239/1‑150 (MQ=255)

                                                             gtgCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTc                                                                                         >  2:2359014/1‑150 (MQ=255)

                                                                ctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGAt                                                                                      <  1:2343276/150‑1 (MQ=255)

                                                                 tgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATa                                                                                     <  1:220475/150‑1 (MQ=255)

                                                                   ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAcc                                                                                   >  1:2813471/1‑150 (MQ=255)

                                                                   ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAcc                                                                                   >  2:550779/1‑150 (MQ=255)

                                                                     aGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCAc                                                                                 >  1:4211649/1‑150 (MQ=255)

                                                                      ggTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACt                                                                                >  1:1564420/1‑150 (MQ=255)

                                                                       gTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTc                                                                               <  2:2459436/150‑1 (MQ=255)

                                                                            ttCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAgg                                                                          >  2:3726809/1‑150 (MQ=255)

                                                                             tCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGc                                                                         <  2:4081100/150‑1 (MQ=255)

                                                                               gCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCt                                                                       >  2:2414177/1‑150 (MQ=255)

                                                                                 cGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCg                                                                     <  1:3580367/150‑1 (MQ=255)

                                                                                   aGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggc                                                                   <  2:3012307/150‑1 (MQ=255)

                                                                                    ggTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggct                                                                  <  1:5073237/150‑1 (MQ=255)

                                                                                      tCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctgg                                                                <  1:795141/150‑1 (MQ=255)

                                                                                      tCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctgg                                                                >  1:1109941/1‑150 (MQ=255)

                                                                                       cTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggc                                                               <  1:2127706/150‑1 (MQ=255)

                                                                                        tGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGACTCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATCCCACTCGAAGGCCACggctggct                                                              >  1:3090041/1‑150 (MQ=255)

                                                                                        tGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggct                                                              >  1:2291114/1‑150 (MQ=255)

                                                                                        tGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggct                                                              >  1:2291096/1‑150 (MQ=255)

                                                                                        tGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggct                                                              >  1:2291115/1‑150 (MQ=255)

                                                                                        tGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACAATGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggct                                                              >  2:1362613/1‑150 (MQ=255)

                                                                                        tGATCGGACCAACCGCAACCACGGGAAGGTCGACCACCAGCGAATCTGTTCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGACTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggct                                                              >  2:1183061/1‑150 (MQ=255)

                                                                                          aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCg                                                                                                                   <  2:4927270/95‑1 (MQ=255)

                                                                                          aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCg                                                                                                                   >  1:4927270/1‑95 (MQ=255)

                                                                                          aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCg                                                            >  2:4656961/1‑150 (MQ=255)

                                                                                          aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCg                                                            >  1:3763945/1‑150 (MQ=255)

                                                                                          aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGAAGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCGCACTCGATACCACTCGAAGGCCTCGGCTGGCTCg                                                            >  2:835508/1‑150 (MQ=255)

                                                                                               aCCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAATCGGCTCTATCAGCTCCGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAg                                                       >  1:2428122/1‑150 (MQ=255)

                                                                                               aCCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAg                                                       >  2:3413174/1‑150 (MQ=255)

                                                                                                 ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGGACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGt                                                     >  1:2097680/1‑150 (MQ=255)

                                                                                                 ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGt                                                     >  2:1847573/1‑150 (MQ=255)

                                                                                                 ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGt                                                     >  1:2177857/1‑150 (MQ=255)

                                                                                                    ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGa                                                  <  1:2041951/150‑1 (MQ=255)

                                                                                                       cAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAgcg                                               >  1:1777897/1‑150 (MQ=255)

                                                                                                       cAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAgcg                                               >  1:69126/1‑150 (MQ=255)

                                                                                                           cACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCAt                                           <  1:1347463/150‑1 (MQ=255)

                                                                                                            aCGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATc                                          >  2:2691413/1‑150 (MQ=255)

                                                                                                             ttcgCAGGTGGTCGACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCt                                         >  2:2140188/4‑150 (MQ=255)

                                                                                                             cGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCt                                         <  2:745606/150‑1 (MQ=255)

                                                                                                                    gTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGc                                  >  1:2690855/1‑150 (MQ=255)

                                                                                                                    gTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGc                                  <  2:3222065/150‑1 (MQ=255)

                                                                                                                        accaccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                              <  2:4071154/150‑1 (MQ=255)

                                                                                                                          caccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctgg                                                                >  2:2623185/1‑114 (MQ=255)

                                                                                                                          caccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctgg                                                                <  1:2623185/114‑1 (MQ=255)

                                                                                                                            ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGc                          <  2:4699984/150‑1 (MQ=255)

                                                                                                                            ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGc                          <  1:2140188/150‑1 (MQ=255)

                                                                                                                            ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGc                          >  1:3305471/1‑150 (MQ=255)

                                                                                                                             cTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGcc                         >  1:3720202/1‑150 (MQ=255)

                                                                                                                              tGCGCATCTGTGCAGCCGCTTCACTGAAGCGGCTCTCTCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                        >  2:2192484/1‑150 (MQ=255)

                                                                                                                              tGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                        >  1:4709500/1‑150 (MQ=255)

                                                                                                                              tGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                        <  1:1716930/150‑1 (MQ=255)

                                                                                                                              tGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGAAACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                        >  1:4709501/1‑150 (MQ=255)

                                                                                                                                  aaTCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCg                    >  1:5109610/1‑150 (MQ=255)

                                                                                                                                    tCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTc                  >  2:1714609/1‑150 (MQ=255)

                                                                                                                                        tgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCgg              >  1:721232/1‑150 (MQ=255)

                                                                                                                                            gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                           <  2:5420161/85‑1 (MQ=255)

                                                                                                                                            gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                           >  1:5420161/1‑85 (MQ=255)

                                                                                                                                            gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTg          >  2:3029498/1‑150 (MQ=255)

                                                                                                                                            gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCGTGTTCGCCCGGCGTCGCGGGTTg          >  2:4702069/1‑150 (MQ=255)

                                                                                                                                            gCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTg          >  1:3752617/1‑150 (MQ=255)

                                                                                                                                            gCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTg          >  2:2019114/1‑150 (MQ=255)

                                                                                                                                             ccGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGc         >  1:2916298/1‑150 (MQ=255)

                                                                                                                                               gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCAt       >  1:526587/1‑150 (MQ=255)

                                                                                                                                                 ttCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATcc     >  1:3018578/1‑150 (MQ=255)

                                                                                                                                                  tCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCa    >  1:4999451/1‑150 (MQ=255)

                                                                                                                                                   cACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCAc   <  2:4488203/150‑1 (MQ=255)

GCCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACT  >  NC_002947/3346162‑3346459

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