BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,057,678	T→G	5.8%	D142A (GAC→GCC)	PP_3572 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,057,678	0	T	G	5.8%	823.3 / 29.1	240	D142A (GAC→GCC)	PP_3572	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (111/114); new base G (8/6); total (119/121)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 240 total aligned reads displayed.

CACGCCGATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCA  >  NC_002947/4057530‑4057825

 aCGCCGATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCg                                                                                                                                                   <  2:385197/150‑1 (MQ=255)

 aCGCCGATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATAGCCGCGTGCATGGCTGCCAGCAGGCGGTCg                                                                                                                                                   <  2:744483/150‑1 (MQ=255)

   gCCGATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGc                                                                                                                                                 <  2:4887218/150‑1 (MQ=255)

      gATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAg                                                                                                                                              <  1:859728/150‑1 (MQ=255)

      gATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGGCGGCCac                                                                                                                                              <  1:1437996/150‑2 (MQ=255)

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

              cGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGgccgc                                                                                                                                      <  1:4551308/150‑1 (MQ=255)

              cGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGgccgc                                                                                                                                      >  2:3298877/1‑150 (MQ=255)

               gCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGgccgac                                                                                                                                     >  2:637704/1‑148 (MQ=255)

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

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

                     ggCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTg                                                                                                                               >  2:5436775/1‑150 (MQ=255)

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

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

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

                       cccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGtt                                                                                                                             >  1:2725701/1‑150 (MQ=255)

                        ccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTc                                                                                                                            >  1:5116705/1‑150 (MQ=255)

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

                          aaTCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaa                                                                                                                          >  2:254176/1‑150 (MQ=255)

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

                            tCGGCGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                        <  2:4116080/150‑1 (MQ=255)

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

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

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

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

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

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

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

                               aCGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGcc                                                                                                                     >  1:2793751/1‑150 (MQ=255)

                               aCGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGcc                                                                                                                     >  1:604208/1‑150 (MQ=255)

                                   gcAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCa                                                                                                                 >  1:1908654/1‑150 (MQ=255)

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

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

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

                                          caccaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGc                                                                                                          >  1:4960549/1‑150 (MQ=255)

                                           accaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCGGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCt                                                                                                         >  1:336169/1‑150 (MQ=255)

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

                                                 ggCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATc                                                                                                   <  1:660016/150‑1 (MQ=255)

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

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

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

                                                       gTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTgg                                                                                             <  1:4844585/150‑1 (MQ=255)

                                                          aGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTc                                                                                          <  1:4074234/150‑1 (MQ=255)

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

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

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

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

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

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

                                                               gCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGa                                                                                     <  2:4162370/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                          tCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                          <  1:2370911/150‑1 (MQ=255)

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

                                                                           cAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGa                                                                                     >  2:1906704/1‑138 (MQ=255)

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

                                                                                  gccgAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAgccgcc                                                                  <  2:3577131/150‑1 (MQ=255)

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

                                                                                          cACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTg                                                          >  2:4294518/1‑150 (MQ=255)

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

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

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

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

                                                                                                    gCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTACAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGcccccc                                                <  2:5116705/150‑5 (MQ=255)

                                                                                                       tCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCAtt                                             <  1:5404172/150‑1 (MQ=255)

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

                                                                                                          ccACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAACCGCCGCGCATTGTTCGCCGGCCATTGCa                                          >  2:1745811/1‑150 (MQ=255)

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

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

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

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

                                                                                                                    ccATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGc                                <  1:1754376/150‑1 (MQ=255)

                                                                                                                      aTCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGccc                              <  1:2688843/150‑1 (MQ=255)

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

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

                                                                                                                          ccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTg                          <  1:3467905/150‑1 (MQ=255)

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

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

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

                                                                                                                                  aTGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTcccgc                  >  2:2663011/1‑150 (MQ=255)

                                                                                                                                   tGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTcccgcc                 >  1:2063258/1‑150 (MQ=255)

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

                                                                                                                                      ttttgaacagGCCGGCGGCCACCTGGCCGACCGCCTGCTCAAAAGCCAGCACCCTGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCgg              >  2:5394287/8‑150 (MQ=38)

                                                                                                                                      cTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGACGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCgg              >  2:3669336/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                          cagcagGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                                   <  2:4616342/77‑1 (MQ=255)

                                                                                                                                           agcagGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAAc         <  1:5549302/150‑1 (MQ=255)

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

                                                                                                                                                gCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTgg                                                                                             <  2:2990427/61‑1 (MQ=255)

CACGCCGATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCA  >  NC_002947/4057530‑4057825

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