BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,871,930	G→C	6.6%	R165R (CGC→CGG)	PP_3420 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,871,930	0	G	C	6.6%	320.1 / 10.1	102	R165R (CGC→CGG)	PP_3420	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base G (47/47); new base C (4/3); total (52/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 104 total aligned reads displayed.

TTGTAGATTTCCACCATGGCCAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAAC  >  NC_002947/3871781‑3872077

ttGTAGATTTCCACCATGGCCAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAA‑ACCGCCTTg                                                                                                                                                     >  1:1408984/1‑150 (MQ=255)

       tttCCACTATGGCCAGCACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGc                                                                                                                                              <  2:1221877/150‑1 (MQ=255)

         tCCACCATGGCCAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCgg                                                                                                                                            >  2:1886825/1‑150 (MQ=255)

           caccaTGGCCAGTACCTGCTCACCCTGGCTGTCGAACATCGGGATGTAGTTCTCGATGAACAGTTAGCGTGGTGCGCGCTGAAACTGCTGCTCCTCACGATCAGCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGaa                                                                                                                                          <  1:921529/150‑1 (MQ=255)

              caTGGCCAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGgc                                                                                                                                       <  2:266822/150‑1 (MQ=255)

               aTGGCCAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGgcg                                                                                                                                      <  2:2196557/150‑1 (MQ=255)

                   ccAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAg                                                                                                                                               <  2:760913/137‑1 (MQ=255)

                   ccAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAg                                                                                                                                               >  1:760913/1‑137 (MQ=255)

                   ccAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTCCGCTCGCGCAAGgcg                                                                                                                                      <  2:2150473/146‑1 (MQ=255)

                   ccAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTCCGCTCGCGCAAGgcg                                                                                                                                      >  1:2150473/1‑146 (MQ=255)

                      gTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCaa                                                                                                                               <  1:283053/150‑1 (MQ=255)

                       tACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAg                                                                                                                              >  1:714483/1‑150 (MQ=255)

                        aCCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGa                                                                                                                             <  1:397704/150‑1 (MQ=255)

                          cTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATc                                                                                                                           >  2:39860/1‑150 (MQ=255)

                          cTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATc                                                                                                                           <  1:1565641/150‑1 (MQ=255)

                              tCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAg                                                                                                                              >  2:1925683/1‑143 (MQ=255)

                              tCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAg                                                                                                                              <  1:1925683/143‑1 (MQ=255)

                                aCCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATc                                                                                                                     >  1:2146008/1‑150 (MQ=255)

                                 cccTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTCCGCTCGCGCAAGgcggt                                                                                                                                    >  2:1377813/1‑132 (MQ=255)

                                 cccTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTCCGCTCGCGCAAGgcggt                                                                                                                                    <  1:1377813/134‑3 (MQ=255)

                                  ccTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAg                                                                                                                   >  1:216299/1‑150 (MQ=255)

                                     ggCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCct                                                                                                                <  2:1980497/150‑1 (MQ=255)

                                        tGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTa                                                                                                             <  2:714483/150‑1 (MQ=255)

                                            gAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCg                                                                                                         >  1:1374858/1‑150 (MQ=255)

                                            gAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCg                                                                                                         >  1:404805/1‑150 (MQ=255)

                                                aGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTgc                                                                                                                                                    <  2:1743447/103‑1 (MQ=255)

                                                 gCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGc                                                                                                    <  1:2354085/150‑1 (MQ=255)

                                                     gATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAAt                                                                                                >  1:2184649/1‑150 (MQ=255)

                                                      aTGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATc                                                                                               >  1:1088844/1‑150 (MQ=255)

                                                           gTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTc                                                                                          >  2:689131/1‑150 (MQ=255)

                                                            ttCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTcc                                                                                         >  1:65811/1‑150 (MQ=255)

                                                               tcGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACg                                                                                      >  1:1210559/1‑150 (MQ=255)

                                                                 gATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCg                                                                                                                          >  2:1202580/1‑112 (MQ=255)

                                                                 gATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCg                                                                                                                          <  1:1202580/112‑1 (MQ=255)

                                                                 gATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGtt                                                                                    >  1:65744/1‑150 (MQ=255)

                                                                  aTGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTtg                                                                                   >  2:849259/1‑150 (MQ=255)

                                                                  aTGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTtg                                                                                   >  1:681867/1‑150 (MQ=255)

                                                                   tGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGACAATCGCCTT‑TTCTATCGGAAGGCGCGCTCAAGATCGACATCCGCCTCTATCCGCTTGCCAATCAGCTCCACGTtgg                                                                                  >  1:2253463/1‑150 (MQ=255)

                                                                   tGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTAC‑CAGCATAAAGAAAAGATCTTGTGATAGCTGGCTGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTtgg                                                                                  >  1:2254540/1‑150 (MQ=255)

                                                                    gAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTtggt                                                                                 <  1:418937/150‑1 (MQ=255)

                                                                     aaCAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCa                                                                                              <  2:1375860/136‑1 (MQ=255)

                                                                     aaCAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCa                                                                                              >  1:1375860/1‑136 (MQ=255)

                                                                     aaCAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTtggtg                                                                                <  2:794735/150‑1 (MQ=255)

                                                                          gTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCa                                                                           <  2:1755396/150‑1 (MQ=255)

                                                                             gCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGAc                                                                        <  2:65744/150‑1 (MQ=255)

                                                                                   tgCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGt                                                                  <  1:39860/150‑1 (MQ=255)

                                                                                    gcgcgcAGAAACTTCTGCGCCTCCCGATCATCCGCA‑ACCTCGTGATAGCTGGCCGAC‑ACCGCCTTGCGCTGGCGGCGGGCGCGCTCATGGCCGCCATCTGCCTCTATCGCGTTGCCTTTCAGGTCCACGTTGGTGGACCAGCCAgtggtg                                                                 <  2:4167/150‑1 (MQ=255)

                                                                                      gcgcTGAAACTTCTGCTCCTCCCGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCCATTCCATGTTGGTGGACCAGACAATGGTGCg                                                               <  1:1880663/150‑1 (MQ=255)

                                                                                        gcTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACTGTTGTGGACCAGACAATGTTGCGGt                                                             >  2:1722020/1‑150 (MQ=255)

                                                                                             aaCTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCgg                                                        <  2:681867/150‑1 (MQ=255)

                                                                                                      tcctcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACa                                               <  1:849259/150‑1 (MQ=255)

                                                                                                        ctcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCa                                                                                                                    <  2:2165785/79‑1 (MQ=255)

                                                                                                         tcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTg                                            <  2:1035659/150‑1 (MQ=255)

                                                                                                         tcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTg                                            >  1:391105/1‑150 (MQ=255)

                                                                                                         tcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTg                                            >  2:207482/1‑150 (MQ=255)

                                                                                                         tcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTg                                            <  2:1374858/150‑1 (MQ=255)

                                                                                                         tcACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTg                                            >  1:2303509/1‑150 (MQ=255)

                                                                                                             gATTATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCa                                        >  2:60185/1‑150 (MQ=255)

                                                                                                             gATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCa                                                                                                                       <  1:1573315/71‑1 (MQ=255)

                                                                                                             gATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCa                                                                                                                       >  2:1573315/1‑71 (MQ=255)

                                                                                                             gATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATcc                                                                                                          <  2:382277/84‑1 (MQ=255)

                                                                                                             gATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATcc                                                                                                          >  1:382277/1‑84 (MQ=255)

                                                                                                                  tcCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAg                                   <  2:2184649/150‑1 (MQ=255)

                                                                                                                  tcCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAg                                   >  2:202540/1‑150 (MQ=255)

                                                                                                                    cGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATa                                                   <  2:2158875/132‑1 (MQ=255)

                                                                                                                    cGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATa                                                   >  1:2158875/1‑132 (MQ=255)

                                                                                                                      cA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCtt                                                                                                      >  1:1643965/1‑79 (MQ=255)

                                                                                                                      cA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCtt                                                                                                      <  2:1643965/79‑1 (MQ=255)

                                                                                                                          ccTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCg                            >  1:164826/1‑150 (MQ=255)

                                                                                                                          ccTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCg                            >  1:1522366/1‑150 (MQ=255)

                                                                                                                          ccTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCg                            <  1:689131/150‑1 (MQ=255)

                                                                                                                           cTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCgg                           >  2:347026/1‑150 (MQ=255)

                                                                                                                            ttGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGc                          >  1:476746/1‑150 (MQ=255)

                                                                                                                              gtgATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAg                        <  2:1522366/150‑1 (MQ=255)

                                                                                                                              gtgATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAg                        <  1:2102975/150‑1 (MQ=255)

                                                                                                                              gtgATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAg                        <  2:476746/150‑1 (MQ=255)

                                                                                                                                gATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGgt                                 <  2:1323306/139‑1 (MQ=255)

                                                                                                                                gATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGgt                                 >  1:1323306/1‑139 (MQ=255)

                                                                                                                                gATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGt                      <  1:202540/150‑1 (MQ=255)

                                                                                                                                      tGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAgg                                  >  2:943944/1‑132 (MQ=255)

                                                                                                                                      tGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAgg                                  <  1:943944/132‑1 (MQ=255)

                                                                                                                                      tGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGa                >  2:2178895/1‑150 (MQ=255)

                                                                                                                                      tGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGa                >  2:779659/1‑150 (MQ=255)

                                                                                                                                           gAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGa           >  1:1043706/1‑150 (MQ=255)

                                                                                                                                            ac‑acCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCAc                                                                                       <  2:1917593/73‑1 (MQ=255)

                                                                                                                                            ac‑acCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCAc                                                                                       >  1:1917593/1‑73 (MQ=255)

                                                                                                                                              gtgcGCCTTCCGCTCGCGCAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGAc               >  1:406442/4‑144 (MQ=255)

                                                                                                                                              gtgcGCCTTCCGCTCGCGCAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGAc               <  2:406442/141‑1 (MQ=255)

                                                                                                                                               acCGCCTTGCGCGAGCGGAAGGCGCGCTCNAGATCGCCATCAGCNTCTANCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGGCATGCTCg        >  1:193018/1‑150 (MQ=255)

                                                                                                                                               acCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAACGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCg        >  1:877972/1‑150 (MQ=255)

                                                                                                                                                ccGCCTTGCGCGAGCGGAAGGCGCGCTCAATATCGCCATCAGCCTCTATCCGCTTTCCTATCGCCACCACGATGGTGGACCATACACTGTTGCTGTCGCGGGCATAGACATTTGCCAACAAGGTGTACCACCGGTACTCGACATGATCGa       <  2:1088844/150‑1 (MQ=255)

                                                                                                                                                 cGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACCTTGGCCAGCAAGGTGTCCGGCAGGGGCGCAACATGATCGAg      >  2:1722003/1‑150 (MQ=255)

                                                                                                                                                 cGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAg      <  1:60185/150‑1 (MQ=255)

                                                                                                                                                 cGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAg      <  1:779659/150‑1 (MQ=255)

                                                                                                                                                 cGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAg      >  1:94825/1‑150 (MQ=255)

                                                                                                                                                 cGCCTTCCGCTCGCGCAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAg      >  1:497143/1‑150 (MQ=255)

                                                                                                                                                  gCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACAGTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAgg     <  2:1311907/150‑1 (MQ=255)

                                                                                                                                                    cTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGaa   >  1:1464008/1‑150 (MQ=255)

                                                                                                                                                     ttGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAAc  <  2:2254540/150‑1 (MQ=255)

                                                                                                                                                     ttGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAAc  <  2:2253463/150‑1 (MQ=255)

TTGTAGATTTCCACCATGGCCAGTACCTGCTCACCCTGGCTGTCGAACAGCGGGATGTAGTTCTCGATGAACAGGTAGCGTGGTGCGCGCTGAAACTTCTGCTCCTCACGATCATCCGCA‑ACCTTGTGATAGCTGGCCGAC‑ACCGCCTTGCGCGAGCGGAAGGCGCGCTCAAGATCGCCATCAGCCTCTATCCGCTTGCCAATCAGCTCCACGTTGGTGGACCAGACAATGGTGCGGTCGCGGGCATAGACATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAAC  >  NC_002947/3871781‑3872077

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