BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,057,686	G→C	5.7%	H139Q (CAC→CAG)	PP_3572 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,057,686	0	G	C	5.7%	829.0 / 29.7	244	H139Q (CAC→CAG)	PP_3572	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (123/107); new base C (9/5); total (132/112)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 252 total aligned reads displayed.

ATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTG  >  NC_002947/4057537‑4057835

aTCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAgg                                                                                                                                                       <  1:2949210/150‑1 (MQ=255)

    ggTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGgc                                                                                                                                                   >  2:5507540/1‑150 (MQ=255)

          tGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCgg                                                                                                                                             <  1:1415202/150‑1 (MQ=255)

           gTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGc                                                                                                                                            <  2:4089786/150‑1 (MQ=255)

             cGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCt                                                                                                                                          <  2:844537/150‑1 (MQ=255)

              ggCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCAATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTg                                                                                                                                         >  2:4678278/1‑150 (MQ=255)

               gCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGt                                                                                                                                        >  1:1018936/1‑150 (MQ=255)

               gCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGt                                                                                                                                        >  1:1018909/1‑150 (MQ=255)

                cccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGtt                                                                                                                                       <  2:1524449/150‑1 (MQ=255)

                 ccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGGCCGACCGCCTgctg                                                                                                                                      <  1:3450461/150‑4 (MQ=255)

                 ccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGGCCGACCGCCTgctg                                                                                                                                      <  1:637704/150‑4 (MQ=255)

                  cAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCa                                                                                                                                     >  2:3758913/1‑150 (MQ=255)

                    aTCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaa                                                                                                                                   <  2:814775/150‑1 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  <  1:2711177/150‑1 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  <  1:3209187/150‑1 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  <  1:3147680/150‑1 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  >  1:4976794/1‑150 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  <  1:2797367/150‑1 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  <  2:1096008/150‑1 (MQ=255)

                     tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                  <  1:169559/150‑1 (MQ=255)

                       gACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGc                                                                                                                                <  2:2534340/150‑1 (MQ=255)

                       gACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGc                                                                                                                                >  1:4060608/1‑150 (MQ=255)

                        aCGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTgg                                                                                                                                                    >  1:1468835/1‑129 (MQ=255)

                             cAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCAc                                                                                                                          <  1:5507540/150‑1 (MQ=255)

                                gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                                       >  2:4778539/1‑150 (MQ=255)

                                gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                                       >  2:361314/1‑150 (MQ=255)

                                 tCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATg                                                                                                                      >  1:742110/1‑150 (MQ=255)

                                  ccaccaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGgccgaccgc                                                                                                                                            <  2:3186143/127‑6 (MQ=255)

                                    accaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCCGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCt                                                                                                                   >  2:364607/1‑150 (MQ=255)

                                    accaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCt                                                                                                                   >  2:4475660/1‑150 (MQ=255)

                                      caGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTc                                                                                                                 <  1:2136324/150‑1 (MQ=255)

                                      caGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTc                                                                                                                 >  2:4923173/1‑150 (MQ=255)

                                           gCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCa                                                                                                            >  2:1596606/1‑150 (MQ=255)

                                               ggTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTg                                                                                                        <  1:4607271/150‑1 (MQ=255)

                                                   aGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTc                                                                                                    >  2:3374247/1‑150 (MQ=255)

                                                     cTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGa                                                                                                  <  1:3645729/150‑1 (MQ=255)

                                                     cTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGa                                                                                                  <  2:274370/150‑1 (MQ=255)

                                                      tCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                                 <  2:3294244/150‑1 (MQ=255)

                                                      tCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                                 <  1:982959/150‑1 (MQ=255)

                                                      tCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                                 <  1:4085711/150‑1 (MQ=255)

                                                        gCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGAt                                                                                                              >  2:570327/1‑135 (MQ=255)

                                                        gCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGAt                                                                                                              <  1:570327/135‑1 (MQ=255)

                                                         ccAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAc                                                                                              >  2:4131211/1‑150 (MQ=255)

                                                         ccAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAc                                                                                              >  2:735909/1‑150 (MQ=255)

                                                          cAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                                             <  2:1187850/150‑1 (MQ=255)

                                                          cAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                                             >  1:1868359/1‑150 (MQ=255)

                                                          cAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                                             >  2:3980101/1‑150 (MQ=255)

                                                           aGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcag                                                                                            <  1:4369900/150‑1 (MQ=255)

                                                           aGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcag                                                                                            >  1:4816522/1‑150 (MQ=255)

                                                             cAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcagca                                                                                          <  2:2509492/150‑1 (MQ=255)

                                                              aaTCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcagcag                                                                                         >  2:2246732/1‑150 (MQ=255)

                                                                  ttCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGATCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTg                                                                                     <  2:2000050/150‑1 (MQ=255)

                                                                   tCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                                    >  2:1723707/1‑150 (MQ=255)

                                                                   tCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAGGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                                    >  1:5568474/1‑150 (MQ=255)

                                                                   tCAGCGCCGCCGAGTGCACGGGCCCGGATTCGCCACCCTCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGATTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                                    <  1:3554295/150‑1 (MQ=255)

                                                                     aGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGtt                                                                                      >  1:238087/1‑146 (MQ=255)

                                                                             cgAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                          <  1:1121329/150‑1 (MQ=255)

                                                                             cgAGTGCACGGGCCCGGCTTCACCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                          <  1:5150078/150‑1 (MQ=255)

                                                                              gAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCcgcg                                                                         >  2:5427884/1‑150 (MQ=255)

                                                                                    aCGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCGGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGt                                                                   >  1:185045/1‑150 (MQ=255)

                                                                                       ggCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCg                                                                >  2:5077869/1‑150 (MQ=255)

                                                                                        gCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGc                                                               <  1:2729986/150‑1 (MQ=255)

                                                                                         cccGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGcc                                                              <  2:3435593/150‑1 (MQ=255)

                                                                                          ccGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCg                                                             >  1:1066298/1‑150 (MQ=255)

                                                                                           cGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCgg                                                            >  2:4209401/1‑150 (MQ=255)

                                                                                             gCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGcc                                                          >  2:3508419/1‑150 (MQ=255)

                                                                                             gCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGcc                                                          >  1:2040440/1‑150 (MQ=255)

                                                                                                   ccACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCa                                                    >  1:4289020/1‑150 (MQ=255)

                                                                                                     aCCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATg                                                  >  2:5125234/1‑150 (MQ=255)

                                                                                                     aCCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATg                                                  <  2:4857070/150‑1 (MQ=255)

                                                                                                       ccGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGgt                                                >  1:4837726/1‑150 (MQ=255)

                                                                                                         gCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGgtgt                                              <  2:2944410/150‑1 (MQ=255)

                                                                                                              cATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGcc                                         <  2:5451408/150‑1 (MQ=255)

                                                                                                               aTCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGccc                                        >  2:1130244/1‑150 (MQ=255)

                                                                                                                cgcccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCa                                       >  1:4541417/4‑150 (MQ=255)

                                                                                                                  gccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGt                                     >  2:2660119/1‑150 (MQ=255)

                                                                                                                       gTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTc                                >  1:2548976/1‑150 (MQ=255)

                                                                                                                        tGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCGGTTCAAAAGCCAGCACCATGGCTTCGATCACCGGGGTCGAGGACAGCAGGGTGCCAGCCGCCGCGCATTGTTCGCCGGCCCGTGCATGGTGTGTGCCCAGTGCTTcc                               >  2:1323457/1‑150 (MQ=255)

                                                                                                                         gCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTccc                              <  2:3354486/150‑1 (MQ=255)

                                                                                                                         gCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTccc                              <  1:352159/150‑1 (MQ=255)

                                                                                                                             ggCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTcccgccc                          >  1:4412835/1‑150 (MQ=255)

                                                                                                                                tGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGt                       <  2:875837/150‑1 (MQ=255)

                                                                                                                                tGCCAGCAGGCGGTCGGCCAGGGGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGt                       >  1:2232435/1‑150 (MQ=255)

                                                                                                                                    agcagGCGGTCGGCCAGGTGGCCGCCGGCCGGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAAc                   <  2:1484541/150‑1 (MQ=255)

                                                                                                                                       agGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAgg                <  2:1002866/150‑1 (MQ=255)

                                                                                                                                        ggCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAggg               <  1:4923173/150‑1 (MQ=255)

                                                                                                                                           ggTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCa            >  1:1800728/1‑150 (MQ=255)

                                                                                                                                              cGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCt         >  1:3001186/1‑150 (MQ=255)

                                                                                                                                              cGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCt         >  1:3676756/1‑150 (MQ=255)

                                                                                                                                                gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                          <  2:3814308/83‑1 (MQ=255)

                                                                                                                                                gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                          >  1:3814308/1‑83 (MQ=255)

                                                                                                                                                gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg       >  2:71332/1‑150 (MQ=255)

                                                                                                                                                gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg       >  1:1976115/1‑150 (MQ=255)

                                                                                                                                                gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg       >  2:3671954/1‑150 (MQ=255)

                                                                                                                                                gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg       >  1:4883729/1‑150 (MQ=255)

                                                                                                                                                 ccAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg       >  1:3356291/1‑149 (MQ=255)

                                                                                                                                                 ccAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGc      >  2:2150516/1‑150 (MQ=255)

                                                                                                                                                 ccAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGc      >  1:3825815/1‑150 (MQ=255)

                                                                                                                                                    ggTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGCCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGCTGAACAGGGCCACCTGGCCGt   <  2:1908654/150‑1 (MQ=255)

                                                                                                                                                     gTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTg  <  1:2926177/150‑1 (MQ=255)

ATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTG  >  NC_002947/4057537‑4057835

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