BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,199	G→C	5.1%	I263M (ATC→ATG)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,199	0	G	C	5.1%	655.0 / 18.1	195	I263M (ATC→ATG)	dnaG	DNA primase
Reads supporting (aligned to +/- strand): ref base G (91/94); new base C (2/8); total (93/102)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 207 total aligned reads displayed.

CGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGCCCG  >  NC_002947/474051‑474346

cGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACgatga                                                                                                                                                     >  1:4489192/1‑150 (MQ=255)

  ccGTCAAAACAGACCAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGAtc                                                                                                                                                   <  1:289318/150‑1 (MQ=255)

    gTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGAtc‑tc                                                                                                                                                <  2:4323012/150‑1 (MQ=255)

       aaaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtc                                                                                                                                             <  2:1872941/150‑1 (MQ=255)

       aaaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtc                                                                                                                                             <  2:452765/150‑1 (MQ=255)

        aaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtcg                                                                                                                                            >  1:5067700/1‑150 (MQ=255)

         aacagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGa                                                                                                                                           <  1:4689536/150‑1 (MQ=255)

         aacagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATATAGCCCTCGACGACGATGATC‑TCGTCGa                                                                                                                                           >  1:4939582/1‑150 (MQ=255)

          acagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAg                                                                                                                                          >  1:4229598/1‑150 (MQ=255)

          acagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                                                                          <  1:1384053/150‑1 (MQ=255)

          acagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                                                                          <  2:836745/150‑1 (MQ=255)

            agaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGt                                                                                                                                        >  1:4241103/1‑150 (MQ=255)

              aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGACCTCGACGA‑GATCATCGTCGTCGAgggct                                                                                                                                      <  2:2181820/150‑4 (MQ=255)

              aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGACCTCGACGA‑GATCATCGTCGTCGAgggct                                                                                                                                      <  2:1532781/150‑4 (MQ=255)

                cagCACGCTTGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCg                                                                                                                                    >  2:1885613/1‑150 (MQ=255)

                 cgCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCgg                                                                                                                                   <  2:3344235/149‑1 (MQ=255)

                 agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCTCTCGACGA‑GATCATCGTCGTCGAgggctaca                                                                                                                                   <  1:3743794/150‑7 (MQ=255)

                 agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCgg                                                                                                                                   <  2:3253168/150‑1 (MQ=255)

                  gCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGt                                                                                                                                  >  1:322178/1‑150 (MQ=255)

                  gCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGt                                                                                                                                  >  2:3184644/1‑150 (MQ=255)

                   cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGtg                                                                                                                                 >  2:4884159/1‑149 (MQ=255)

                     cGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTtgt                                                                                                                               <  2:4105182/150‑1 (MQ=255)

                         gggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCtt                                                                                                                           <  2:148421/150‑1 (MQ=255)

                         gggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCtt                                                                                                                           <  2:750461/150‑1 (MQ=255)

                              ccacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcg                                                                                                                      >  1:5367778/1‑150 (MQ=255)

                               cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                                     >  1:2471214/1‑150 (MQ=255)

                               cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                                     <  1:2790573/150‑1 (MQ=255)

                               cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                                     >  2:1691897/1‑150 (MQ=255)

                                acGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGcc                                                                                                                    >  1:2546915/1‑150 (MQ=255)

                                  gcgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTc                                                                                                                  <  1:5357293/150‑1 (MQ=255)

                                  gcgGAACAGGCGCTGGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTc                                                                                                                  <  2:2669925/150‑1 (MQ=255)

                                        cAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAg                                                                                                            <  2:2965538/150‑1 (MQ=255)

                                        cAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAg                                                                                                            <  1:5022037/150‑1 (MQ=255)

                                           gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                         >  2:4606428/1‑150 (MQ=255)

                                           gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                         >  2:1963020/1‑150 (MQ=255)

                                              cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                                                                      <  1:4959051/150‑1 (MQ=255)

                                              cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                                                                      <  1:1185449/150‑1 (MQ=255)

                                                tGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACa                                                                                                    <  2:3811675/150‑1 (MQ=255)

                                                 gAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAg                                                                                                   <  1:3879477/150‑1 (MQ=255)

                                                        tctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGc                                                                                            <  1:5301755/150‑1 (MQ=255)

                                                        tctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGc                                                                                            >  2:2300367/1‑150 (MQ=255)

                                                        tctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGCACAGCCCGTACAGCTCCTGc                                                                                            >  1:4122246/1‑150 (MQ=255)

                                                          ttCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGccc                                                                                          <  2:4838000/150‑1 (MQ=255)

                                                          ttCGCTGGGGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGccc                                                                                          <  2:2905666/150‑1 (MQ=255)

                                                               tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGGCAGGGCAATGAGGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCGTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                     >  2:5159410/1‑150 (MQ=255)

                                                               tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                     >  1:4770396/1‑150 (MQ=255)

                                                               tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                     >  2:2590119/1‑150 (MQ=255)

                                                               tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                     >  1:4770589/1‑150 (MQ=255)

                                                                     cGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGaaaa                                                                               >  1:3515383/1‑150 (MQ=255)

                                                                          ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                          <  1:3063867/150‑1 (MQ=255)

                                                                          ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                          <  1:2680571/150‑1 (MQ=255)

                                                                          ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                          <  2:1668956/150‑1 (MQ=255)

                                                                           cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgggg                                                                         >  1:1812678/1‑150 (MQ=255)

                                                                            aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                        >  2:832652/1‑150 (MQ=255)

                                                                            aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                        <  1:1316729/150‑1 (MQ=255)

                                                                            aaGGGTGGCCACGGCATAGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                        >  2:522614/1‑150 (MQ=255)

                                                                               ggTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTc                                                                     >  2:1108573/1‑150 (MQ=255)

                                                                                gTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTcc                                                                    >  1:799926/1‑150 (MQ=255)

                                                                                 tGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCg                                                                   >  2:3926293/1‑150 (MQ=255)

                                                                                     cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                               >  2:183117/1‑150 (MQ=255)

                                                                                     cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                               <  1:3184644/150‑1 (MQ=255)

                                                                                      aCGTCATTGCCCTGGCCCAGCAAGGCCTGCGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAg                                                              >  1:2310035/1‑150 (MQ=38)

                                                                                      aCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCg                                                                   <  2:3548362/145‑1 (MQ=255)

                                                                                      aCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAg                                                              >  2:4690279/1‑150 (MQ=255)

                                                                                         gCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTc                                                           <  2:3039296/150‑1 (MQ=255)

                                                                                                gcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACt                                                    >  2:1711170/1‑150 (MQ=255)

                                                                                                 cAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACtt                                                   >  1:176466/1‑150 (MQ=255)

                                                                                                      cTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCt                                              >  1:1472239/1‑150 (MQ=255)

                                                                                                           tGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTgtcg                                         <  2:4442346/150‑1 (MQ=255)

                                                                                                               ccAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGc                                     >  2:1284403/1‑150 (MQ=255)

                                                                                                                cAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGcc                                    >  1:3283266/1‑150 (MQ=255)

                                                                                                                 aGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCg                                   <  1:3565802/150‑1 (MQ=255)

                                                                                                                   ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAg                                 <  2:3515383/150‑1 (MQ=255)

                                                                                                                   ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAg                                 >  1:1558995/1‑150 (MQ=255)

                                                                                                                   ggCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAg                                 >  2:1758316/1‑150 (MQ=255)

                                                                                                                         gACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                           >  2:2827288/1‑150 (MQ=255)

                                                                                                                         gACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                           <  2:4995867/150‑1 (MQ=255)

                                                                                                                            gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGcc                        >  1:3545752/1‑150 (MQ=255)

                                                                                                                             tCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCa                       >  1:2981237/1‑150 (MQ=255)

                                                                                                                                aTGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCg                    >  2:3630340/1‑150 (MQ=255)

                                                                                                                                  gTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGaa                  >  1:3931059/1‑150 (MQ=255)

                                                                                                                                   tAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAg                 <  2:2513613/150‑1 (MQ=255)

                                                                                                                                     gCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGc               <  2:4939582/150‑1 (MQ=255)

                                                                                                                                      cccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCg              >  1:948437/1‑150 (MQ=255)

                                                                                                                                      cccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCg              <  2:111365/150‑1 (MQ=255)

                                                                                                                                       ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGaa                  <  1:4424877/145‑1 (MQ=255)

                                                                                                                                       ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga             >  2:592491/1‑150 (MQ=255)

                                                                                                                                       ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga             >  2:2895117/1‑150 (MQ=255)

                                                                                                                                       ccTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga             >  1:1986472/1‑150 (MQ=255)

                                                                                                                                       ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                           >  1:1597964/1‑88 (MQ=38)

                                                                                                                                        cTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgat            <  1:2396076/150‑1 (MQ=255)

                                                                                                                                         tCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg           >  2:1452833/1‑150 (MQ=255)

                                                                                                                                         tCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg           >  1:646035/1‑150 (MQ=255)

                                                                                                                                         tCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg           >  2:2396597/1‑150 (MQ=255)

                                                                                                                                               acgaTGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATgcgc     <  2:4411949/150‑1 (MQ=255)

                                                                                                                                               acgaTGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTGTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATgcgc     <  1:4755038/150‑1 (MQ=255)

                                                                                                                                                 gatgatC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGccc   >  1:4977452/1‑150 (MQ=255)

                                                                                                                                                 gatgatC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGccc   >  1:3800291/1‑150 (MQ=255)

                                                                                                                                                 gatgatC‑TCGTCGAGGTTGCGGTTGCGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTggg                                                >  1:2263414/1‑105 (MQ=255)

                                                                                                                                                  atgatC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCGCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGCCcg  >  2:3543334/1‑150 (MQ=255)

CGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGCCCG  >  NC_002947/474051‑474346

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